Genomic Analysis Research Articles

CFPG: Creating a Common Fungal Pathogenic Genes Database through Data Mining

Fu ngal pathogenicity is one of the most vigorously tackled ecological and medicinal issues facing many scientists. Comparative genomics is an extremely important methodology and tool used to understand fungal pathogenicity, and it allows the development of early diagnostic tools for fungal-inflicted diseases across different host organisms. However, comparative genomics depends heavily on readily available fungal pathogenic gene databases to enable downstream genomics study and the development of new diagnosis and detection methods. Here, we have developed the Common Fungal Pathogenic Genes Database through comparative genomic analysis using 86 publicly available fungal genomic data against fungal pathogenicity-related databases, such as Pathogenic-Host Interaction Database (PHI-base), Carbohydrate-Active enZymes Database (CAZy), and Database of Fungal Virulence Factory (DFVF).

ITRAQ-based proteomics revealed the acaricidal mechanism of methyl isoeugenol against Dermatophagoides farinae

ABSTRACT Dermatophagoides farinae is a major source of allergens that are the cause of allergic diseases in humans. Methyl isoeugenol (MIE) is a naturally occurring compound in plants and its strong acaricidal toxicity has been established; however, the mechanism behind its acaricidal properties remains unclear. In this study, to investigate the action mechanism of MIE on D. farinae the isobaric tags for relative and absolute quantification (iTRAQ) proteomics method is used. Protein samples were extracted from D. farinae before and after 24 h of MIE treatment for iTRAQ technique analysis. Bioinformatics analysis revealed that most differentially expressed proteins were located in the cytoplasm, impacting cytoplasmic function. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed significant enrichment of the differential proteins in the cellular autophagy signalling pathway. The autophagy-related proteins Autophagy-Related Protein 3 (ATG3) and Cathepsin L (CTSL) were up-regulated, disrupting the autophagic process. We found that the main action mechanism of MIE against D. farinae might be cellular autophagy. The transcriptional validation of six differential proteins demonstrated consistent expression trends with the proteomics data. This study reveals the acaricidal mechanism of MIE for the first time and lays the foundation for further research on D. farinae control.

Functional genomic analysis of the isolated potential probiotic Lactobacillus delbrueckii subsp. indicus TY-11 and its comparison with other Lactobacillus delbrueckii strains.

This study aimed to conduct functional genomic analysis to uncover the probiotic potential of Lactobacillus delbrueckii subsp. indicus TY-11 isolated from native yogurt in Bangladesh. We also performed transmission electron microscopic (TEM) analysis, comparative genomic study as well as phylogenetic tree construction with 332 core genes from 262 genomes. In our current investigation, we revealed a number of common and unique excellences of the probiotic Lactobacillus delbrueckii subsp. indicus TY-11 that are likely to be important to illustrate its intestinal residence and probiotic roles. This is the first study to explore the molecular insights into intestinal residence and probiotic roles, including antimicrobial activities and antibiotic sensitivity, of a representative strain (TY-11) of Lactobacillus delbrueckii subsp. indicus.

Comprehensive comparative analysis of the periodontal pathogen Porphyromonas gingivalis: exploring the pan-genome, the reconstruction of the gene regulatory network and genome-scale metabolic network.

Porphyromonas gingivalis is a nonmotile, obligate anaerobic, Gram-negative bacterium known for its association with periodontal disease and its involvement in systemic diseases such as atherosclerosis, cardiovascular disease, colon cancer and Alzheimer's disease. This bacterium produces several virulence factors, including capsules, fimbriae, lipopolysaccharides, proteolytic enzymes and hemagglutinins. A comparative genomic analysis revealed the open pangenome of P. gingivalis and identified complete type IV secretion systems (T4SS) in strain KCOM2805 and almost complete type VI secretion systems (T6SS) in strains KCOM2798 and ATCC49417, which is a new discovery as previous studies did not find the proteins involved in secretion systems IV and VI. Conservation of some virulence factors between different strains was observed, regardless of their genetic diversity and origin. In addition, we performed for the first time a reconstruction analysis of the gene regulatory network (GRN), identifying transcription factors and proteins involved in the regulatory mechanisms of bacterial pathogenesis. In particular, QseB regulates the expression of hemagglutinin and arginine deaminase, while Rex may suppress the release of gingipain through interactions with PorV and the formatum/nitrate transporter. Our study highlights the central role of conserved virulence factors and regulatory pathways, particularly QseB and Rex, in P. gingivalis and provides insights into potential therapeutic targets.

Raman cell sorting for single-cell research

Cells constitute the fundamental units of living organisms. Investigating individual differences at the single-cell level facilitates an understanding of cell differentiation, development, gene expression, and cellular characteristics, unveiling the underlying laws governing life activities in depth. In recent years, the integration of single-cell manipulation and recognition technologies into detection and sorting systems has emerged as a powerful tool for advancing single-cell research. Raman cell sorting technology has garnered attention owing to its non-labeling, non-destructive detection features and the capability to analyze samples containing water. In addition, this technology can provide live cells for subsequent genomics analysis and gene sequencing. This paper emphasizes the importance of single-cell research, describes the single-cell research methods that currently exist, including single-cell manipulation and single-cell identification techniques, and highlights the advantages of Raman spectroscopy in the field of single-cell analysis by comparing it with the fluorescence-activated cell sorting (FACS) technique. It describes various existing Raman cell sorting techniques and introduces their respective advantages and disadvantages. The above techniques were compared and analyzed, considering a variety of factors. The current bottlenecks include weak single-cell spontaneous Raman signals and the requirement for a prolonged total cell exposure time, significantly constraining Raman cell sorting technology’s detection speed, efficiency, and throughput. This paper provides an overview of current methods for enhancing weak spontaneous Raman signals and their associated advantages and disadvantages. Finally, the paper outlines the detailed information related to the Raman cell sorting technology mentioned in this paper and discusses the development trends and direction of Raman cell sorting.

Occurrence Autosomal Mosaicism in Karan Fries Cattle and Murrah Buffalo with Poor Reproduction

Background: The present study was conducted at the Livestock Genome Analysis Laboratory (LGAL), ICAR-National Dairy Research Institute (NDRI), Karnal, Haryana, India, to explore, chromosomal anomalies in reproductively inefficient animals. Chromosomal abnormalities are deviations in normal genetic architecture and lead to disorders in bearer individuals. These can be both in number and structure of autosomes or sex chromosomes, usually inflict adverse effects on anatomy, survival and reproductive performance in domestic animals. Methods: Present report is on a female Murrah buffalo (Bubalus bubalis) and a Karan Fries cattle which were identified during a routine cytogenetic reproductively inefficient females at cattle yard, National Dairy Research Institute, Karnal. Metaphase chromosome preparations were prepared using standard lymphocyte culture method. Slides were stained with Giemsa and treated for R-banding using fluorochrome-photolysis-Giemsa technique (RB-FPG technique). Karyotypes were prepared and abnormal chromosome was identified. Result: Routine examination using Giemsa-stained slides revealed two cell lines with 60 and 61 chromosomal constitutions in cattle and 3 cell lines in buffalo. Among 182 metaphase plates examined were examined for Karan fries which revealed two types of chromosome configuration, the frequency of cells carrying 60 and 61 chromosomes was 83.5% and 16.5%, respectively. R banding technique distinguished autosome 16 responsible for mosaicism. Cytogenetic evaluation of 407 metaphase spreads of Murrah buffalo revealed three types of chromosome configurations, viz. 49, 50 and 51 all with XX and the overall frequency was 16.7, 76.7 and 6.6 per cent, respectively. The R- bandings confirmed autosome 11 both in monosomy (2n=49) and trisomy (2n=51), respectively. This unusual chromosomal number might have arisen due to non-disjunction during early stage of zygotic development of the buffalo and cattle.

Detection and evaluation of susceptibility to antibiotics in non-hydrogen sulfide-producing antibiotic-resistant soil microbe: Pseudomonas guariconensis.

Antimicrobial resistance in bacteria is a global threat that can make antibacterial treatments ineffective. One well-known method of antibiotic resistance and a common defensive mechanism in many harmful bacteria is the synthesis of endogenous hydrogen sulfide (H2S) in bacteria. In this study, soil bacteria were screened using the lead acetate agar test and the triple sugar iron test to determine that they were non-endogenous H2S producers. This was further validated by full genome analysis of the identified organism against the gene sequences of H2S-producing genes. Antibacterial resistance of the bacteria was phenotypically analyzed using the Kirby-Bauer disk diffusion method. Then, the effect of exogenous H2S on the antibiotic-resistant bacteria was checked in sodium sulfide, leading to antibiotic re-sensitization.

New Solutions in Single-Cell Protein Production from Methane: Construction of Glycogen-Deficient Mutants of Methylococcus capsulatus MIR

The biotechnology of converting methane to single-cell protein (SCP) implies using fast-growing thermotolerant aerobic methanotrophic bacteria. Among the latter, members of the genus Methylococcus received significant research attention and are used in operating commercial plants. Methylococcus capsulatus MIR is a recently discovered member of this genus with the potential to be used for the purpose of SCP production. Like other Methylococcus species, this bacterium stores carbon and energy in the form of glycogen, particularly when grown under nitrogen-limiting conditions. The genome of strain MIR encodes two glycogen synthases, GlgA1 and GlgA2, which are only moderately related to each other. To obtain glycogen-free cell biomass of this methanotroph, glycogen synthase mutants, ΔglgA1, ΔglgA2, and ΔglgA1ΔglgA2, were constructed. The mutant lacking both glycogen synthases exhibited a glycogen-deficient phenotype, whereas the intracellular glycogen content was not reduced in strains defective in either GlgA1 or GlgA2, thus suggesting functional redundancy of these enzymes. Inactivation of the glk gene encoding glucokinase also resulted in a sharp decrease in glycogen content and accumulation of free glucose in cells. Wild-type strain MIR and the mutant strain ΔglgA1ΔglgA2 were also grown in a bioreactor operated in batch and continuous modes. Cell biomass of ΔglgA1ΔglgA2 mutant obtained during batch cultivation displayed high protein content (71% of dry cell weight (DCW) compared to 54% DCW in wild-type strain) as well as a strong reduction in glycogen content (10.8 mg/g DCW compared to 187.5 mg/g DCW in wild-type strain). The difference in protein and glycogen contents in biomass of these strains produced during continuous cultivation was less pronounced, yet biomass characteristics relevant to SCP production were slightly better for ΔglgA1ΔglgA2 mutant. Genome analysis revealed the presence of glgA1-like genes in all methanotrophs of the Gammaproteobacteria and Verrucomicrobia, while only a very few methanotrophic representatives of the Alphaproteobacteria possessed these determinants of glycogen biosynthesis. The glgA2-like genes were present only in genomes of gammaproteobacterial methanotrophs with predominantly halo- and thermotolerant phenotypes. The role of glycogen in terms of energy reserve is discussed.

Genomic identification and evolutionary analysis of chemosensory receptor gene families in two Phthorimaea pest species: insights into chemical ecology and host adaptation

BackgroundInsects rely on sophisticated sensitive chemosensory systems to sense their complex chemical environment. This sensory process involves a combination of odorant receptors (ORs), gustatory receptors (GRs) and ionotropic receptors (IRs) in the chemosensory system. This study focused on the identification and characterization of these three types of chemosensory receptor genes in two closely related Phthorimaea pest species, Phthorimaea operculella (potato tuber moth) and Phthorimaea absoluta (tomato leaf miner).ResultsBased on manual annotation of the genome, we identified a total of 349 chemoreceptor genes from the genome of P. operculella, including 93 OR, 206 GR and 50 IR genes, while for P. absoluta, we identified 72 OR, 122 GR and 46 IR genes. Through phylogenetic analysis, we observed minimal differences in the number and types of ORs and IRs between the potato tuber moth and tomato leaf miner. In addition, we found that compared with those of tomato leaf miners, the gustatory receptor branch of P. operculella has undergone a large expansion, which may be related to P. absoluta having a narrower host range than P. operculella. Through analysis of differentially expressed genes (DEGs) of male and female antennae, we uncovered 45 DEGs (including 32ORs, 9 GRs, and 4 IRs).ConclusionsOur research provides a foundation for exploring the chemical ecology of these two pests and offers new insights into the dietary differentiation of lepidopteran insects, while simultaneously providing molecular targets for developing environmentally friendly pest control methods based on insect chemoreception.

Abstract A009: Integrative drug screening and multi-omic characterization of patient-derived bladder cancer organoids reveals novel molecular correlates of chemotherapy response

Abstract Introduction: Predicting response to therapy for each patient’s tumor is critical to improving long-term outcomes for muscle-invasive bladder cancer (MIBC). Our objective is to establish ex vivo bladder cancer patient-derived organoid (PDO) models that are representative of patients’ tumors and determine potential efficacy of standard-of-care and curated experimental therapies. Methods: Tumor material was prospectively collected from consented patients with bladder cancer to generate short-term PDO models, which were screened against a panel of clinically relevant drugs in ex vivo 3D culture. Drugs were tested at the maximum plasma concentration (Cmax) in human trials, so as to provide physiologic relevance. When material was sufficient, cells were additionally tested using full dose response. The number of drugs screened was dependent upon amount of tissue available, with up to 34 drugs screened at Cmax and 9 drugs screened in dose response format. Multiomic profiling was utilized to validate the PDO models, establish the molecular characteristics of each tumor, assess gene expression (GEX) patterns between paired primary tissue and PDOs, and identify potential biomarkers of drug response using Pearson correlation coefficients and Kruskal-Wallis tests with Dunn’s post-hoc pairwise comparison testing. Results: 106 tumors were collected from 97 patients, with 65 samples yielding sufficient material for complete multiomic molecular characterization and PDO screening with 6 to 32 drugs/combinations. RNA sequencing and GSEA enriched pathways across organoids were compared to other published datasets, noting high levels of correlation between our organoid models and in particular the Cancer Cell Line Encyclopedia urethral lines. Moreover, analysis of subtype and mutation/copy number data further indicate that our organoid dataset is representative of the full spectrum of disease. RNA sequencing and subtyping of cultured organoids compared to patient tissue indicate that organoid models are representative of patient tissue and do not undergo subtype shifts in our short-term 3D cultures. Our drug screening results highlight the large spectrum of response to chemotherapies that in fact, clinically benefit only a small proportion of patients in the neoadjuvant setting. Utilizing an integrative approach, novel correlations between ex vivo drug responses and genomic alterations, GEX, and protein expression were identified, including a multiomic signature of gemcitabine response. Conclusions: Rapid organoid development, characterization, and drug screening allows for the prediction of therapeutic response in ~10 days following sample collection. Use of this technique on tissue provided during disease work-up may further guide selection of effective therapeutic agents in patients with bladder cancer, beyond the typical unranked lists derived from genomic analyses. This work would help maximize the overall benefits of standard of care therapies and could be used to identify alternative therapeutics for patients who progress on standard therapies. Citation Format: Nathan M. Merrill, Samuel D. Kaffenberger, Liwei Bao, Nathalie Vandecan, Laura E. Goo, Athena Apfel, Xu Cheng, Zhaoping Qin, Chia-Jen Liu, Armand Bankhead, Yin Wang, Varun Kathawate, Lila Tudrick, Habib Serhan, Zackariah Farah, Chad Ellimoottil, Khaled S. Hafez, Lindsey A. Herrel, Jeffrey S. Montgomery, Todd M. Morgan, Simpa S. Salami, Alon Z. Weizer, Peter J. Ulintz, Mark L. Day, Matthew B. Soellner, Phillip L. Palmbos, Sofia D. Merajver, Aaron M. Udager. Integrative drug screening and multi-omic characterization of patient-derived bladder cancer organoids reveals novel molecular correlates of chemotherapy response [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2024 May 17-20; Charlotte, NC. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(10_Suppl):Abstract nr A009.

Framework for exploring the sensory repertoire of the human gut microbiota.

Bacteria sense changes in their environment and transduce signals to adjust their cellular functions accordingly. For this purpose, bacteria employ various sensors feeding into multiple signal transduction pathways. Signal recognition by bacterial sensors is studied mainly in a few model organisms, but advances in genome sequencing and analysis offer new ways of exploring the sensory repertoire of many understudied organisms. The human gut is a natural target of this line of study: it is a nutrient-rich and dynamic environment and is home to thousands of bacterial species whose activities impact human health. Many gut commensals are also poorly studied compared to model organisms and are mainly known through their genome sequences. To begin exploring the signals human gut commensals sense and respond to, we have designed a framework that enables the identification of sensory domains, prediction of signals that they recognize, and experimental verification of these predictions. We validate this framework's functionality by systematically identifying amino acid sensors in selected bacterial genomes and metagenomes, characterizing their amino acid binding properties, and demonstrating their signal transduction potential.IMPORTANCESignal transduction is a central process governing how bacteria sense and respond to their environment. The human gut is a complex environment with many living organisms and fluctuating streams of nutrients. One gut inhabitant, Escherichia coli, is a model organism for studying signal transduction. However, E. coli is not representative of most gut microbes, and signaling pathways in the thousands of other organisms comprising the human gut microbiota remain poorly understood. This work provides a foundation for how to explore signals recognized by these organisms.

Genomics and transcriptomics reveal β-carotene synthesis mechanism in Dunaliella salina

Dunaliella salina is by far the most salt-tolerant organism and contains many active substances, including β-carotene, glycerol, proteins, and vitamins, using in the production of dried biomass or cell extracts for the biofuels, pharmaceutical formulations, food additives, and fine chemicals, especially β-carotene. We report a high-quality genome sequence of D. Salina FACHB435, which has a 472 Mb genome size, with a contig N50 of 458 Kb. A total of 30,752 protein-coding genes were predicted. The annotation results evaluated by BUSCO was shown that completeness was 91.0% and replication was 53.1%. The fragments were 6.3% and the deletions were 2.6%. Phylogenomic and comparative genomic analyses revealed that A. thaliana diverged from Volvocales about 448 million years ago, then Volvocales C. eustigma, D. salina, and other species diverged about 250 million years ago. High light could promote the accumulation of β-carotene in D. salina at a 13 d stage of culture. The enrichment of DEGs in KEGG, it notes that the predicted up-regulated genes of carotenoid metabolic pathway include DsCrtB, DsPDS, DsZ-ISO, DsZDS, DsCRTISO, DsLUT5, DsCrtL-B, and DsCCD8, while the predicted down-regulated genes include DsCrtF, and DsLUT1. The four genes that were both up-regulated and down-regulated were DsZEP, DsCrtR-b, DsCruA/P and DsCrtZ 4. The research results can provide scientific basis for the industrialization practice of D. salina.

The whole genome sequence of Polish vaccine strain Mycobacterium bovis BCG Moreau.

The whole genome sequence obtained is the only genomic sequence of the strain that has been used for vaccine production in Poland since 1955. Sequencing of different BCG lots showed that the strain was stable over a period of 59 years. The comprehensive genomic analysis performed not only enriches knowledge about the microevolution and attenuation of the BCG vaccine substrains but also enables the utilization of identified markers as a reference point in the genetic control and identity tests of the stability of the vaccine strain in the future.

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.

Comparative 3D genome analysis between neural retina and retinal pigment epithelium reveals differential cis-regulatory interactions at retinal disease loci

BackgroundVision depends on the interplay between photoreceptor cells of the neural retina and the underlying retinal pigment epithelium (RPE). Most genes involved in inherited retinal diseases display specific spatiotemporal expression within these interconnected retinal components through the local recruitment of cis-regulatory elements (CREs) in 3D nuclear space.ResultsTo understand the role of differential chromatin architecture in establishing tissue-specific expression at inherited retinal disease loci, we mapped genome-wide chromatin interactions using in situ Hi-C and H3K4me3 HiChIP on neural retina and RPE/choroid from human adult donor eyes. We observed chromatin looping between active promoters and 32,425 and 8060 candidate CREs in the neural retina and RPE/choroid, respectively. A comparative 3D genome analysis between these two retinal tissues revealed that 56% of 290 known inherited retinal disease genes were marked by differential chromatin interactions. One of these was ABCA4, which is implicated in the most common autosomal recessive inherited retinal disease. We zoomed in on retina- and RPE-specific cis-regulatory interactions at the ABCA4 locus using high-resolution UMI-4C. Integration with bulk and single-cell epigenomic datasets and in vivo enhancer assays in zebrafish revealed tissue-specific CREs interacting with ABCA4.ConclusionsThrough comparative 3D genome mapping, based on genome-wide, promoter-centric, and locus-specific assays of human neural retina and RPE, we have shown that gene regulation at key inherited retinal disease loci is likely mediated by tissue-specific chromatin interactions. These findings do not only provide insight into tissue-specific regulatory landscapes at retinal disease loci, but also delineate the search space for non-coding genomic variation underlying unsolved inherited retinal diseases.Graphical

Upgraded durian genome reveals the role of chromosome reshuffling during ancestral karyotype evolution, lignin biosynthesis regulation, and stress tolerance.

Durian (Durio zibethinus) is a tropical fruit that has a unique flavor and aroma. It occupies a significant phylogenetic position within the Malvaceae family. Extant core-eudicot plants are reported to share seven ancestral karyotypes that have undergone reshuffling, resulting in an abundant genomic diversity. However, the ancestral karyotypes of the Malvaceae family, as well as the evolution trajectory leading to the 28 chromosomes in durian, remain poorly understood. Here, we report the high-quality assembly of the durian genome with comprehensive comparative genomic analyses. By analyzing the collinear blocks between cacao and durian, we inferred 11 Malvaceae ancestral karyotypes. These blocks were present in a single-copy form in cacao and mainly in triplicates in durian, possibly resulting from a recent whole genome triplication (WGT) event that led to hexaploidization of the durian genome around 20 (17-24) million years ago. A large proportion of the duplicated genes in durian, such as those involved in the lignin biosynthesis module for phenylpropane biosynthesis, are derived directly from whole genome duplication, which makes it an important force in reshaping its genomic architecture. Transcriptome studies have revealed that genes involved in feruloyl-CoA formations were highly preferentially expressed in fruit peels, indicating that the thorns produced on durian fruit may comprise guaiacyl and syringyl lignins. Among all the analyzed transcription factors (TFs), members of the heat shock factor family (HSF) were the most significantly upregulated under heat stress. All subfamilies of genes encoding heat shock proteins (HSPs) in the durian genome appear to have undergone expansion. The potential interactions between HSF Dzi05.397 and HSPs were examined and experimentally verified. Our study provides a high-quality durian genome and reveals the reshuffling mechanism of ancestral Malvaceae chromosomes to produce the durian genome. We also provide insights into the mechanism underlying lignin biosynthesis and heat stress tolerance.

The chromosome-level genome and functional database accelerate research about biosynthesis of secondary metabolites in Rosa roxburghii

Rosa roxburghii Tratt, a valuable plant in China with long history, is famous for its fruit. It possesses various secondary metabolites, such as L-ascorbic acid (vitamin C), alkaloids and poly saccharides, which make it a high nutritional and medicinal value. Here we characterized the chromosome-level genome sequence of R. roxburghii, comprising seven pseudo-chromosomes with a total size of 531 Mb and a heterozygosity of 0.25%. We also annotated 45,226 coding gene loci after masking repeat elements. Orthologs for 90.1% of the Complete Single-Copy BUSCOs were found in the R. roxburghii annotation. By aligning with protein sequences from public platform, we annotated 85.89% genes from R. roxburghii. Comparative genomic analysis revealed that R. roxburghii diverged from Rosa chinensis approximately 5.58 to 13.17 million years ago, and no whole-genome duplication event occurred after the divergence from eudicots. To fully utilize this genomic resource, we constructed a genomic database RroFGD with various analysis tools. Otherwise, 69 enzyme genes involved in L-ascorbate biosynthesis were identified and a key enzyme in the biosynthesis of vitamin C, GDH (L-Gal-1-dehydrogenase), is used as an example to introduce the functions of the database. This genome and database will facilitate the future investigations into gene function and molecular breeding in R. roxburghii.

P32 Melanin confers protection against ultraviolet radiation-mediated mitochondrial DNA damage in skin of colour

Abstract Introduction and aims Melanin is a photoprotective pigment that effectively absorbs and scatters solar ultraviolet radiation (UVR) affording protection against nuclear DNA photodamage. In lightly pigmented skin, UVR-mediated damage of cellular DNA is not limited to nuclear DNA; mitochondrial DNA (mtDNA) is also susceptible to UVR-induced damage leading to an increased frequency of deletions and mutations. However, the protective role of melanin against UVR-induced mtDNA damage is unknown. Therefore, we sought to determine whether melanin affords protection against UVR-induced mtDNA damage. Methods Using Genome Analysis Toolkit best practices, mtDNA was extracted from photoexposed forearm and photoprotected buttock of Black African (n = 3) and White Northern European (n = 3) volunteers and sequenced to identify somatic mutations. Results Comparison of photoexposed forearm with photoprotected buttock identified a total of 21 somatic mutations, including single-nucleotide transitions, deletions and additions. Of these, three somatic mutations were observed in Black skin, whereas 18 were identified in White skin (P < 0.01). Analysis of single-nucleotide transitions revealed a total of 13 variations and of these, two were identified in Black skin and 11 in White skin. Interestingly, of the somatic mutations observed in White skin, there was enrichment of the cytosine to thymine (C>T) transitions in all volunteers. In contrast, photoexposed Black skin had a thymine to guanine transition at position 414 (T414G) identified for two of the three Black volunteers. Conclusions Here we identified that Black skin exhibited fewer mtDNA mutations compared with White skin; enrichment of C>T transitions in White skin is consistent with a previously reported UVR-signature of mtDNA damage, whereas the T414G transition identified in Black photoexposed skin has previously been suggested as a mitochondrial biomarker of photoageing. Our study suggests that melanin confers protection against UVR-mediated mtDNA damage and that the mtDNA damage signature differs between skin types.

Genome-wide identification and expression analysis of Ubiquitin-specific protease gene family in maize (Zea mays L.)

BackgroundUbiquitin-specific proteases (UBPs) are a large family of deubiquitinating enzymes (DUBs). They are widespread in plants and are critical for plant growth, development, and response to external stresses. However, there are few studies on the functional characteristics of the UBP gene family in the important staple crop, maize (Zea mays L.).ResultsIn this study, we performed a bioinformatic analysis of the entire maize genome and identified 45 UBP genes. Phylogenetic analysis indicated that 45 ZmUBP genes can be divided into 15 subfamilies. Analysis of evolutionary patterns and divergence levels indicated that ZmUBP genes were present before the isolation of dicotyledons, were highly conserved and subjected to purifying selection during evolution. Most ZmUBP genes exhibited different expression levels in different tissues and developmental stages. Based on transcriptome data and promoter element analysis, we selected eight ZmUBP genes whose promoters contained a large number of plant hormones and stress response elements and were up-regulated under different abiotic stresses for RT-qPCR analysis, results showed that these genes responded to abiotic stresses and phytohormones to varying degrees, indicating that they play important roles in plant growth and stress response.ConclusionsIn this study, the structure, location and evolutionary relationship of maize UBP gene family members were analyzed for the first time, and the ZmUBP genes that may be involved in stress response and plant growth were identified by combining promoter element analysis, transcriptome data and RT-qPCR analysis. This study informs research on the involvement of maize deubiquitination in stress response.

Whole genome and transcriptome integrated analyses guide clinical care of pediatric poor prognosis cancers

The role for routine whole genome and transcriptome analysis (WGTA) for poor prognosis pediatric cancers remains undetermined. Here, we characterize somatic mutations, structural rearrangements, copy number variants, gene expression, immuno-profiles and germline cancer predisposition variants in children and adolescents with relapsed, refractory or poor prognosis malignancies who underwent somatic WGTA and matched germline sequencing. Seventy-nine participants with a median age at enrollment of 8.8 y (range 6 months to 21.2 y) are included. Germline pathogenic/likely pathogenic variants are identified in 12% of participants, of which 60% were not known prior. Therapeutically actionable variants are identified by targeted gene report and whole genome in 32% and 62% of participants, respectively, and increase to 96% after integrating transcriptome analyses. Thirty-two molecularly informed therapies are pursued in 28 participants with 54% achieving a clinical benefit rate; objective response or stable disease ≥6 months. Integrated WGTA identifies therapeutically actionable variants in almost all tumors and are directly translatable to clinical care of children with poor prognosis cancers.