Expression Of Various Genes Research Articles

Variability in non-tumor areas of colorectal cancer patients as revealed by endoscopic intestinal step biopsies

A comprehensive endoscopic small and large intestinal untargeted step biopsy procedure was conducted to compare gene expression between the normal intestinal mucosa of healthy individuals and that of patients with colorectal tumors. From 78 participants (healthy individuals [n = 17], patients with colorectal conventional adenomas [n = 6], patients with Tis–T1 colorectal cancer [n = 41], patients with T2–4 colorectal cancer [n = 14]), biopsies of normal mucosa of the terminal ileum, right-sided colon (cecum and ascending colon), and left-sided colorectum (descending colon, sigmoid colon, and rectum) were obtained using a lower gastrointestinal endoscope. RNA was extracted from all samples, and total transcriptome sequencing was performed. Transcriptome data from 388 samples was analyzed. DNA was also extracted from tumor biopsy tissues and analyzed for whole-exome sequencing. In healthy individuals, gene expression differed significantly among the terminal ileum, right-sided colon, and left-sided colorectum, presumably linked to embryological factors. There were differences in gene expression in the normal mucosa in colorectal cancer patients, compared to healthy controls. Patients with tumors, especially T2–4 colorectal cancer, showed considerable variation in gene expression in non-tumor tissues, even in the terminal ileum distant from the tumor site. Based on endoscopic biopsies, the results imply cancer-predisposing conditions in seemingly normal tissues. The present study points to the importance of small intestine and cancer-predisposing conditions in the colon of colorectal cancer patients, with possible implications for developing novel immunotherapy and other therapeutic modalities.

Simultaneous single-nucleus RNA sequencing and single-nucleus ATAC sequencing of neuroblastoma cell lines

Neuroblastoma is the most common extracranial solid tumor in children, and a leading cause of childhood cancer deaths. All neuroblastomas arise from neural crest-derived sympathetic neuronal progenitors, but numerous mutations, the most common of which is MYCN amplification, give rise to these lesions. Epigenetic aberrations also play a role in oncogenesis and tumor progression. To better understand biologic diversity of neuroblastomas, we performed joint single-nucleus ATAC sequencing and single-nucleus RNA sequencing on six neuroblastoma cell lines, three of which are MYCN amplified. After standard filtering for high-quality nuclei, we obtained chromatin accessibility and transcript abundance data from 41,733 neuroblastoma tumor cells. Preliminary analysis reveals significant diversity in chromatin landscape and gene expression across neuroblastoma cell lines. This dataset is a valuable resource for studying the transcriptional and epigenetic mechanisms of this deadly childhood disease.

Comparative transcriptomic insights into the evolutionary origin of the tetrapod double cone.

The tetrapod double cone is a pair of tightly associated cones called the "principal" and the "accessory" member. It is found in amphibians, reptiles, and birds, as well as monotreme and marsupial mammals but is absent in fish and eutherian mammals. To explore the potential evolutionary origins of the double cone, we analyzed single-cell and -nucleus transcriptomic atlases of photoreceptors from six vertebrate species: zebrafish, chicken, lizard, opossum, ground squirrel, and human. Computational analyses separated the principal and accessory members in chicken and lizard, identifying molecular signatures distinguishing either member from single cones and rods in the same species. Comparative transcriptomic analyses suggest that both the principal and accessory originated from ancestral red cones. Furthermore, the gene expression variation among cone subtypes mirrors their spectral order (red → green → blue → UV), suggesting a constraint in their order of emergence during evolution. Finally, we find that rods are equally dissimilar to all cone types, suggesting that they emerged before the spectral diversification of cones.

Insights into molecular mechanisms of phytochemicals in quorum sensing modulation for bacterial biofilm control.

Biofilm formation is a common mechanism by which bacteria undergo phenotypic changes to adapt to environmental stressors. The formation of biofilms has a detrimental impact in clinical settings by contributing to chronic infections and promoting antibiotic resistance. Delving into the molecular mechanisms, the quorum sensing (QS) system involves the release of chemical signals for bacterial cell-to-cell communication, which activates and regulates the expression of various genes and virulence factors, including those related to biofilm formation. Accordingly, the QS system becomes a potential target for combating biofilm-associated concerns. Natural products derived from plants have a long history of treating infectious diseases in humans due to their antimicrobial properties, making them valuable resources for screening anti-biofilm agents. This review aims to discover the mechanisms by which phytochemical agents inhibit QS, potentially offering promising new therapies for treating biofilm-associated infections. By targeting the QS system, these phytochemical agents can prevent bacterial aggregation and biofilm formation while also diminishing other bacterial virulence factors. Additionally, it is important to focus on the advancement of techniques and experiments to investigate their molecular mechanisms. A thorough understanding of these mechanisms may encourage further studies to evaluate the safety and efficacy of phytochemical agents used alone or in combination with other strategies.

Pan-transcriptome analysis provide insights into resistance and fruit quality breeding of pear (Pyrus pyrifolia)

Pyrus pyrifolia, commonly known as sand pear, is a key economic fruit tree in temperate regions and possesses highly diverse germplasm resources for pear quality improvement. However, research on the relationship between resistance and fruit quality traits during the breeding of fruit species like pear is limited. Pan-transcriptomes effectively capture genetic information from coding regions and reflect variations in gene expression between individuals. Here, we constructed a pan-transcriptome composed of 506 samples from different tissues of sand pear, and explored the intrinsic relationships among phenotypes and the selection for disease resistance during improvement based on expression presence/absence variation (ePAVs). The pan-transcriptome contains 156,744 transcripts, among which the novel transcripts show significant enrichment in the defense response. Interestingly, disease resistance genes are highly expressed in landraces of pear but have been selected against during the improvement of this perennial tree species. We found that genetically diverse landraces can be divided into two subgroups and inferred that they have undergone different dispersal processes. Through co-expression network analysis, we confirmed that the formation of stone cells in pears, the deposition of fruit anthocyanins, and the ability to resist stress are interrelated. They are jointly regulated by several modules, and the expression of regulatory genes has a significant correlation. Moreover, we identified candidate genes such as HKL1 that may affect sugar content and are missing from the reference genome. Our study provides insights into the associations between complex fruit traits, while providing a database resource for pear disease resistance and fruit quality breeding.

Structural variation reshapes population gene expression and trait variation in 2,105 Brassica napus accessions

Although individual genomic structural variants (SVs) are known to influence gene expression and trait variation, the extent and scale of SV impact across a species remain unknown. In the present study, we constructed a reference library of 334,461 SVs from genome assemblies of 16 representative morphotypes of neopolyploid Brassica napus accessions and detected 258,865 SVs in 2,105 resequenced genomes. Coupling with 5 tissue population transcriptomes, we uncovered 285,976 SV-expression quantitative trait loci (eQTLs) that associate with altered expression of 73,580 genes. We developed a pipeline for the high-throughput joint analyses of SV-genome-wide association studies (SV-GWASs) and transcriptome-wide association studies of phenomic data, eQTLs and eQTL-GWAS colocalization, and identified 726 SV–gene expression–trait variation associations, some of which were verified by transgenics. The pervasive SV impact on how SV reshapes trait variation was demonstrated with the glucosinolate biosynthesis and transport pathway. The study highlighting the impact of genome-wide and species-scale SVs provides a powerful methodological strategy and valuable resources for studying evolution, gene discovery and breeding.

Microevolution of longan flowering behaviors from primitive to KClO3 induction is driven by artificial selection for agricultural benefit

Dimocarpus longan is an economically important fruit tree primarily cultivated in various Asian countries and the Indochina region encompassing Thailand, China, Vietnam, Laos and Cambodia. Here, we report an extensive dataset of 606.72Gb of the whole-genome resequencing data obtained from 31 Dimocarpus longan and 1 Dimocarpus obtusus samples performed on Illumina HiSeq PE150, providing the coverage over 30-fold of the reference genome (averaging 15Gb per sample) with an addition of 12 longan accessions from China. A phylogenomic tree inferred from this analyzed dataset (19,270,513 SNPs) was constructed and results indicated that there are three separate monophyletic clades of China-USA (1), Thai (2) and Vietnam longan (3). Interestingly, however, Thai and Vietnam clades appear to be closely linked in genetic relationship based on analyses of karyotype, genomic data and certain morphological characteristics. De novo transcriptome assembly further revealed variations in gene sequences, divergent gene expression and candidate genes associated with different flowering behaviors. The greatest change in differential gene expression (51.87%) observed between natural independent and seasonal flowering implied the microevolutionary shift observed in longan is influenced by artificial selection, resulting in gradual changes in flowering behaviors from natural independent flowering to artificial flowering induction using potassium chlorate (KClO3). In summary, the data obtained from this study serves as the essential evidence for elucidating the microevolution of longan and shedding new light on an agronomical application of artificial flowering induction via modulation of KClO3 responsive genes in longan or other fruits in the future.

A DNA base-specific sequence interposed between CRX and NRL contributes to RHODOPSIN expression

Gene expression emerges from DNA sequences through the interaction of transcription factors (TFs) with DNA cis-regulatory sequences. In eukaryotes, TFs bind to transcription factor binding sites (TFBSs) with differential affinities, enabling cell-specific gene expression. In this view, DNA enables TF binding along a continuum ranging from low to high affinity depending on its sequence composition; however, it is not known whether evolution has entailed a further level of entanglement between DNA-protein interaction. Here we found that the composition and length (22 bp) of the DNA sequence interposed between the CRX and NRL retinal TFs in the proximal promoter of RHODOPSIN (RHO) largely controls the expression levels of RHO. Mutagenesis of CRX-NRL DNA linking sequences (here termed “DNA-linker”) results in uncorrelated gene expression variation. In contrast, mutual exchange of naturally occurring divergent human and mouse Rho cis-regulatory elements conferred similar yet species-specific Rho expression levels. Two orthogonal DNA-binding proteins targeted to the DNA-linker either activate or repress the expression of Rho depending on the DNA-linker orientation relative to the CRX and NRL binding sites. These results argue that, in this instance, DNA itself contributes to CRX and NRL activities through a code based on specific base sequences of a defined length, ultimately determining optimal RHO expression levels.

Evolutionary innovations in germline biology of placental mammals identified by transcriptomics of first-wave spermatogenesis in opossum.

Mammalian spermatogenesis is a highly stereotyped and conserved developmental process that is essential for fitness. At the same time, gene expression in spermatogenic cells is rapidly evolving. This combination of features has been suggested to drive rapid fixation of new gene expression patterns. Using a high-resolution dataset comprising bulk and single-cell data from juvenile and adult testes of the opossum Monodelphis domestica, a model marsupial, we define the developmental timing of the spermatogenic first wave in opossum and delineate conserved and divergent gene expression programs across the placental-marsupial split by comparison to equivalent data from mouse, a model placental mammal. Epigenomic data confirmed divergent regulation at the level of transcription, and comparison to data from four additional amniote species identified hundreds of genes with evidence of rapid fixation of expression. This gene set encompasses known and previously undescribed regulators of spermatogenic development.

Pan-genome analysis reveals genomic variations during enoki mushroom domestication, with emphasis on genetic signatures of cap color and stipe length

IntroductionThe domestication of edible mushrooms, including Flammulina filiformis, offers valuable insights into the genetic changes driven by artificial selection. Understanding these changes is crucial for uncovering the mechanisms behind genome evolution in domesticated mushrooms. ObjectivesThis study aims to investigate the population structure, genetic diversity, and domestication-related genomic changes in F. filiformis. By comparing the genome sequences of 199 wild and cultivated strains, we aim to elucidate the impact of domestication on F. filiformis. MethodsWe performed de novo genome assembly and gene-based pan-genome analysis on the 199 strains, which included both wild and cultivated strains. We also conducted genome-wide association studies (GWAS) using presence-absence variation (PAV) and SNP data, combined with RNA sequencing, to identify genes associated with domestication traits, such as cap color and stipe length. Gene functional confirmation was achieved through genetic transformation experiments. ResultsOur analysis grouped the strains into four distinct populations, which correlated with varying intensities of artificial selection. The three cultivated populations exhibited smaller genome sizes, fewer genes, lower genetic variation, reduced gene expression diversity, and lower heterozygosity compared to the wild population. The analysis revealed the loss of genes related to the beta-lactam antibiotic catabolic process and specific MAPK pathway genes during domestication, rendering domesticated strains more susceptible to diseases. Four genes closely associated with cap color and stipe length were identified, but genetic transformation experiments confirmed the functional relevance of only two (FfB and FfD) identified through PAV-based GWAS. ConclusionThis study uncovered significant genomic variations between cultivated and wild F. filiformis populations, including the loss of pathogen resistance genes during domestication. We also identified key genes linked to cap color and stipe length, demonstrating for the first time the important role of PAV variation in mushroom domestication. These insights provide a foundation for future mushroom breeding and evolutionary research.

Transcriptomic and bioinformatic analysis of Histone H2A-genes in Lates calcarifer and Oreochromis niloticus upon Vibrio vulnificus challenge

Histone H2A-derived antimicrobial peptides (AMPs) have been recognized for their antimicrobial properties and role in host defense mechanisms. The present study investigated the transcriptional responses of two histone H2A-derived AMPs, Lc-Histone from Lates calcarifer and GIFT-Histone from Oreochromis niloticus, in response to Vibrio vulnificus infection. The genes encoding these peptides were amplified and sequenced from fingerling stage fishes, followed by in silico analysis and phylogenetic characterization using ClustalW. Computational methods were used to assess their functional attributes, and in vivo bacterial challenge assay was employed to evaluate the efficacy. Sequence alignment identified conserved histone motifs, supporting their evolutionary conservation. The antimicrobial peptide database predicted membrane interaction potential, suggesting antimicrobial activity, alongside biofilm-active properties and efficient expression in heterologous protein systems. Expression profiling of Histone H2A-derived AMP genes under Gram-negative bacterial stimulation revealed distinct temporal dynamics, indicative of their innate immune response involvement. Significant gene expression variations were observed between L. calcarifer and O. niloticus. These baseline findings enhance understanding of host defense mechanisms and bacterial interactions in aquaculture. The antimicrobial potential of histone H2A-derived AMPs highlights their significance as therapeutic agents and warrants further investigation into their mechanisms of action and optimal usage in aquaculture settings.

Transcriptomic analysis identifies enrichment of cAMP/PKA/CREB signaling in invasive lobular breast cancer

ObjectiveInvasive lobular breast cancer (ILC) is the most common special type of breast cancer and has unique clinicopathological and molecular hallmarks that differentiate it from the more common invasive carcinoma—no special type (NST). Despite these differences, ILC and NST are treated as a single entity and there is a lack of ILC-targeted therapies. To fill this gap, we sought to identify novel molecular alterations in ILC that could be exploited for targeted therapies.MethodsDifferential gene expression and Geneset Enrichment and Variation analyses were performed on RNA-seq data from three large public breast cancer databases—the Sweden Cancerome Analysis Network-Breast (SCAN-B; luminal A ILC N = 263, luminal A NST N = 1162), The Cancer Genome Atlas (TCGA; luminal A ILC N = 157, luminal A NST N = 307) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC; luminal A ILC N = 65, luminal A NST N = 533). Pathways enriched in overlapping differentially expressed genes from these datasets were clustered using Jaccard similarity to identify pathways enriched in ILC. The cAMP/PKA/CREB signaling was studied in ILC, ILC-like and NST cell lines and patient-derived organoids (PDOs) using forskolin, an activator of the pathway.ResultsClinicopathological features of patients with ILC and NST in SCAN-B were similar to prior population-based studies. There was a consistent pattern of up-regulation of cAMP/PKA/CREB related signaling in ILC compared to NST in SCAN-B, TCGA and METABRIC. Treatment with forskolin resulted in a greater increase in phospho-CREB in ILC cell lines and organoids than NST. CRISPR deletion of CDH1 in NST cell lines did not alter response of cells to forskolin as measured by phospho-CREB. Forskolin treatment caused growth inhibition in ILC and NST, with ILC cell lines being more sensitive to forskolin-mediated growth inhibition.ConclusionIn three separate datasets, cAMP/PKA/CREB signaling was identified to be higher in ILC than NST. This in silico finding was validated in cell line and organoid models. Loss of CDH1 was not sufficient to mediate this phenotype. Future studies should investigate the mechanisms for differential cAMP/PKA/CREB signaling and the potential for therapeutic targeting in patients with ILC.Graphical abstract

Lineage tracing of nuclei in skeletal myofibers uncovers distinct transcripts and interplay between myonuclear populations

Multinucleated skeletal muscle cells need to acquire additional nuclei through fusion with activated skeletal muscle stem cells when responding to both developmental and adaptive growth stimuli. A fundamental question in skeletal muscle biology has been the reason underlying this need for new nuclei in cells that already harbor hundreds of nuclei. Here we utilize nuclear RNA-sequencing approaches and develop a lineage tracing strategy capable of defining the transcriptional state of recently fused nuclei and distinguishing this state from that of pre-existing nuclei. Our findings reveal the presence of conserved markers of newly fused nuclei both during development and after a hypertrophic stimulus in the adult. However, newly fused nuclei also exhibit divergent gene expression that is determined by the myogenic environment to which they fuse. Moreover, accrual of new nuclei through fusion is required for nuclei already resident in adult myofibers to mount a normal transcriptional response to a load-inducing stimulus. We propose a model of mutual regulation in the control of skeletal muscle development and adaptations, where newly fused and pre-existing myonuclear populations influence each other to maintain optimal functional growth.

Human-specific protein-coding and lncRNA genes cast sex-biased genes in the brain and their relationships with brain diseases

BackgroundGene expression shows sex bias in the brain as it does in other organs. Since female and male humans exhibit noticeable differences in emotions, logical thinking, movement, spatial orientation, and even the incidence of neurological disorders, sex biases in the brain are especially interesting, but how they are determined, whether they are conserved or lineage specific, and what the consequences of the biases are, remain poorly explored and understood.MethodsBased on RNA-seq datasets from 16 and 14 brain regions in humans and macaques across developmental periods and from patients with brain diseases, we used linear mixed models (LMMs) to differentiate variations in gene expression caused by factors of interest and confounding factors and identify four types of sex-biased genes. Effect size and confidence in each effect were measured upon the local false sign rate (LFSR). We utilized the biomaRt R package to acquire orthologous genes in humans and macaques from the BioMart Ensembl website. Transcriptional regulation of sex-biased genes by sex hormones and lncRNAs were analyzed using the CellOracle, GENIE3, and Longtarget programs. Sex-biased genes’ functions were revealed by gene set enrichment analysis using multiple methods.ResultsLineage-specific sex-biased genes greatly determine the distinct sex biases in human and macaque brains. In humans, those encoding proteins contribute directly to immune-related functions, and those encoding lncRNAs intensively regulate the expression of other sex-biased genes, especially genes with immune-related functions. The identified sex-specific differentially expressed genes (ssDEGs) upon gene expression in disease and normal samples also indicate that protein-coding ssDEGs are conserved in humans and macaques but that lncRNA ssDEGs are not conserved. The results answer the above questions, reveal an intrinsic relationship between sex biases in the brain and sex-biased susceptibility to brain diseases, and will help researchers investigate human- and sex-specific ncRNA targets for brain diseases.ConclusionsHuman-specific genes greatly cast sex-biased genes in the brain and their relationships with brain diseases, with protein-coding genes contributing to immune response related functions and lncRNA genes critically regulating sex-biased genes. The high proportions of lineage-specific lncRNAs in mammalian genomes indicate that sex biases may have evolved rapidly in not only the brain but also other organs.

DNA methylation stochasticity is linked to transcriptional variability and identifies convergent epigenetic disruption across genetically-defined subtypes of AML.

Disruption of the epigenetic landscape is of particular interest in acute myeloid leukemia (AML) due to its relatively low mutational burden and frequent occurrence of mutations in epigenetic regulators. Here, we applied an information-theoretic analysis of methylation potential energy landscapes, capturing changes in mean methylation level and methylation entropy, to comprehensively analyze DNA methylation stochasticity in subtypes of AML defined by mutually exclusive genetic mutations. We identified AML subtypes with CEBPA double mutation and those with IDH mutations as distinctly high-entropy subtypes, marked by methylation disruption over a convergent set of genes. We found a core program of epigenetic landscape disruption across all AML subtypes, with discordant methylation stochasticity and transcriptional dysregulation converging on functionally important leukemic signatures, suggesting a genotype-independent role of stochastic disruption of the epigenetic landscape in mediating leukemogenesis. We further established a relationship between methylation entropy and gene expression variability, connecting the disruption of the epigenetic landscape to transcription in AML. This approach identified a convergent program of epigenetic dysregulation in leukemia, clarifying the contribution of specific genetic mutations to stochastic disruption of the epigenetic and transcriptional landscapes of AML.

Transcriptome reanalysis and gene expression of 13 detoxification genes for avermectin and pyridaben resistance in Panonychus citri

Citrus red mites (P. citri) are key pests affecting citrus production worldwide due to pesticide resistance. The resistance mechanisms of ten pesticides are known, but a comprehensive study using transcriptome data is missing. This study employed deeptools, cuffdiff, rmats, bcftools and other software to examine gene expression variation, alternative splicing (AS), and mutations in mite resistance. The research highlighted that pesticides can regulate gene transcription, and red mites with resistance increase cytochrome P450, glutathione S-transferases, carboxylesterase, and acetylcholinesterase expression. Pyridaben also induces new AS events. Fluazinam-induced mites show mRNA splicing peaking earlier than transcription, both peaking at one day and returning to baseline after two days. AS profiles are similar in different mite populations with overlapping pesticide resistances. Lastly, specific mitochondrial SNPs in mites might mediate resistance against select pesticides.

The Mysterious World of Non-Canonical Caps - What We Know and Why We Need New Sequencing Techniques.

It was long believed that viral and eukaryotic mRNA molecules are capped at their 5' end solely by the N7-methylguanosine cap, which regulates various aspects of the RNA life cycle, from its biogenesis to its decay. However, the recent discovery of a variety of non-canonical RNA caps derived from metabolites and cofactors - such as NAD, FAD, CoA, UDP-glucose, UDP-N-acetylglucosamine, and dinucleoside polyphosphates - has expanded the known repertoire of RNA modifications. These non-canonical caps are found across all domains of life and can impact multiple aspects of RNA metabolism, including stability, translation initiation, and cellular stress responses. The study of these modifications has been facilitated by sophisticated methodologies such as liquid chromatography-mass spectrometry, which have unveiled their presence in both prokaryotic and eukaryotic organisms. The identification of these novel RNA caps highlights the need for advanced sequencing techniques to characterize the specific RNA types bearing these modifications and understand their roles in cellular processes. Unravelling the biological role of non-canonical RNA caps will provide insights into their contributions to gene expression, cellular adaptation, and evolutionary diversity. This review emphasizes the importance of these technological advancements in uncovering the complete spectrum of RNA modifications and their implications for living systems.

Cell-Specific Gene Expressions Underlie Selective White Matter Loss Vulnerability in Mild Traumatic Brain Injury.

Traumatic brain injury (TBI), a risk factor for later-life dementia, leads to salient brain atrophy, particularly in the white matter. It is not clear how white matter atrophy progresses or why some brain regions are damaged while others are spared. We hypothesized that spatial variations of cell-specific gene expression contributed to the selective white matter loss vulnerability following mild TBI (mTBI). Gene expression data were sourced from the publicly available Allen Human Brain Atlas, which comprises microarray data spanning nearly the entire brain, derived from six neurologically normal adult donors. A total of 100 patients with acute stage (within 7 days post-injury) mTBI were enrolled. Of these, 60 patients were followed up at 3 months post-injury and 37 were followed up at 6-12 months post-injury. In addition, 59 healthy controls (HCs), matched for age, gender, and education, were included for comparative analysis. White matter volume changes were analyzed at both the acute stage, 3 months, and 6-12 months follow-up in mTBI patients compared with HCs. Patients with mTBI exhibited significant white matter atrophy in the frontal, parietal, and temporal cortices at 3 months post-injury, which even persisted at 6-12 months follow-up. In addition, mTBI patients with cognitive deficits showed more severe brain atrophy compared with those without cognitive deficits. Crucially, the gene expression marking endothelial cells and S1 pyramidal neurons were associated with increased brain atrophy, whereas the gene expression marking microglia and CA1 pyramidal neurons were associated with decreased brain atrophy in mTBI patients at 3 months post-injury. Microglia and endothelial cells can explain 23.6% of regional variations in the white matter atrophy. These findings suggested that modulating cellular activation, especially by promoting microglial activation at 3 months post-injury, might be a promising approach to prevent white matter atrophy, enhance cognitive outcomes, and reduce the risk of later-life dementia.

Genotyping sequence-resolved copy-number variations using pangenomes reveals paralog-specific global diversity and expression divergence of duplicated genes.

Copy-number variable (CNV) genes are important in evolution and disease, yet sequence variation in CNV genes is a blindspot for large-scale studies. We present a method, ctyper, that leverages pangenomes to produce copy-number maps with allele-specific sequences containing locally phased variants of CNV genes from NGS reads. We extensively characterized accuracy and efficiency on a database of 3,351 CNV genes including HLA , SMN , and CYP2D6 as well as 212 non-CNV medically-relevant challenging genes. The genotypes capture 96.5% of underlying variants in new genomes, requiring 0.9 seconds per gene. Expression analysis of ctyper genotypes explains more variance than known eQTL variants. Comparing allele-specific expression quantified divergent expression on 7.94% of paralogs and tissue-specific biases on 4.7% of paralogs. We found reduced expression of SMN-1 converted from SMN-2, which potentially affects diagnosis of spinal muscular atrophy, and increased expression of a duplicative translocation of AMY2B . Overall, ctyper enables biobank-scale genotyping of CNV and challenging genes.

Genomic Localization Bias of Secondary Metabolite Gene Clusters and Association with Histone Modifications in Aspergillus.

Fungi are well-known producers of bioactive secondary metabolites (SMs), which have been exploited for decades by humankind for various medical applications like therapeutics and antibiotics. SMs are synthesized by biosynthetic gene clusters (BGCs)-physically co-localized and co-regulated genes. Because BGCs are often regulated by histone post-translational modifications (PTMs), it was suggested that their chromosomal location is important for their expression. Studies in a few fungal species indicated an enrichment of BGCs in sub-telomeric regions; however, there is no evidence that BGCs with distinct genomic localization are regulated by different histone PTMs. Here, we used 174 Aspergillus species covering 22 sections to determine the correlation between BGC genomic localization, gene expression, and histone PTMs. We found a high abundance and diversity of SM backbone genes across the Aspergillus genus, with notable unique genes within sections. Being unique or conserved in many species, BGCs showed a strong bias for being localized in low-synteny regions, regardless of their position in chromosomes. Using chromosome-level assemblies, we also confirmed a significantly biased localization in sub-telomeric regions. Notably, SM backbone genes in sub-telomeric regions and about half of those in low-synteny regions exhibit higher gene expression variability, likely due to the similar higher variability in H3K4me3 and H3K36me3 histone PTMs; while variations in histone H3 acetylation and H3K9me3 are not correlated to genomic localization and expression variation, as analyzed in two Aspergillus species. Expression variability across four Aspergillus species further supports that BGCs tend to be located in low-synteny regions and that regulation of expression in those regions likely involves different histone PTMs than the most commonly studied modifications.