The UCSC genome browser and associated tools

Harnessing machine learning technique to authenticate differentially expressed genes in oral squamous cell carcinoma

OBJECTIVE

Advancements in early detection of the disease, prognosis and the development of therapeutic strategies necessitate tumor-specific biomarkers. Despite continuous efforts, no molecular marker has been proven to be an effective therapeutic tool for the early detection of cancer. The study aims to determine diagnostic and prognostic signature genes that may be involved in cancer pathology and hence, may serve as molecular markers.

STUDY DESIGN

Eight candidate genes were selected based on our prior study of transcriptomic sequencing and validated in 100 matched pair samples of oral squamous cell carcinoma (OSCC). We further utilized machine learning approaches and examined the diagnostic presentation and predictive ability of the OSCC genes retrieved from publicly available The Cancer Genome Atlas (TCGA) database and compared with our results.

RESULTS

We conducted qPCR analysis to validate the expression of each gene and observed that each gene was present in the majority of OSCC samples. The predictive ability of selected genes was stable (with an average accuracy of 84%) across different classifiers. However, on validation with our dataset, it showed 75% accuracy, which might be because of the demographic variation of the samples.

CONCLUSIONS

The present research outlines cancer-associated molecular biomarkers that might eventually contribute to an enhanced prognosis of cancer patient by identifying novel therapeutic targets.

Concurrent genotyping of mitochondrial DNA and nuclear DNA in rootless hair shafts and blood samples for enhanced analysis

Hair is an important type of biological evidence at crime scenes. However, the highly degraded nature of DNA fragments in hair shafts poses challenges for the detection of nuclear DNA (nuDNA) through capillary electrophoresis-based short tandem repeat (STR) genotyping. In this study, an all-in-one multiplex system named MGIEasy Signature Identification Library Prep Kit (MGI Tech, China) was employed to the simultaneous genotyping of both mitochondrial DNA (mtDNA) and nuDNA in hair shafts. This system is based on massively parallel sequencing (MPS) technology and encompasses Amelogenin, STRs, single nucleotide polymorphisms (SNPs) and mtDNA hypervariable regions (HVRs) in a single reaction. A total of 370 hair shafts, together with 180 blood samples as the references, were examined. The mtDNA analysis of 110 unrelated blood samples unveiled a total of 150 homoplasmic variants and 105 distinct haplotypes, revealing population polymorphisms in the Guangdong Han Chinese. The study also delved into the detection of mtDNA heteroplasmy, revealing 8.18 % and 16.36 % of individuals with point heteroplasmies (PHPs) in blood and hair shaft samples, respectively. Additionally, hair shafts with DNA extracted using the Investigator method yielded higher average depth of coverage (DoC), lower drop-out rate for SNP genotyping, higher nuDNA genotyped rates and success rates, than those using the MinElute method. In the longitudinal research, a gradual decrease in the total DoC of mtDNA fragments was observed along the length of the hair shaft, from the proximal root to the distal end. In contrast, the DoC of nuDNA exhibited a relatively stable pattern along the length of the hair shafts. The study contributes valuable insights into the simultaneous detection of nuDNA and mtDNA in hair shafts, emphasizing the need for optimized DNA extraction and detection methods for these highly degraded samples.

Rare Variant Analyses in Ancestrally Diverse Cohorts Reveal Novel ADHD Risk Genes

Attention-deficit/hyperactivity disorder (ADHD) is a highly heritable neurodevelopmental disorder, but its genetic architecture remains incompletely characterized. Rare coding variants, which can profoundly impact gene function, represent an underexplored dimension of ADHD risk. In this study, we analyzed large-scale DNA sequencing datasets from ancestrally diverse cohorts and observed significant enrichment of rare protein-truncating and deleterious missense variants in highly evolutionarily constrained genes. This analysis identified 15 high-confidence ADHD risk genes, including the previously implicated KDM5B . Integrating these findings with genome-wide association study (GWAS) data revealed nine enriched pathways, with strong involvement in synapse organization, neuronal development, and chromatin regulation. Protein-protein interaction analyses identified chromatin regulators as central network hubs, and single-cell transcriptomic profiling confirmed their expression in neurons and glial cells, with distinct patterns in oligodendrocyte subtypes. These findings advance our understanding of the genetic architecture of ADHD, uncover core molecular mechanisms, and provide promising directions for future therapeutic development.

Transcriptome-wide association study identifies genes associated with bladder cancer risk

Genome-wide association studies (GWAS) have detected several susceptibility variants for urinary bladder cancer, but how gene regulation affects disease development remains unclear. To extend GWAS findings, we conducted a transcriptome-wide association study (TWAS) using PrediXcan to predict gene expression levels in whole blood using genome-wide genotype data for 6180 bladder cancer cases and 5699 controls included in the database of Genotypes and Phenotypes (dbGaP). Logistic regression was used to estimate adjusted gene-level odds ratios (OR) per 1-standard deviation higher expression with 95% confidence intervals (CI) for bladder cancer risk. We further assessed associations for individual single-nucleotide polymorphisms (SNPs) used to predict expression levels and proximal loci for genes identified in gene-level analyses with false-discovery rate (FDR) correction. TWAS identified four genes for which expression levels were associated with bladder cancer risk: SLC39A3 (OR = 0.91, CI = 0.87-0.95, FDR = 0.015), ZNF737 (OR = 0.91, CI = 0.88-0.95, FDR = 0.016), FAM53A (OR = 1.09, CI = 1.05-1.14, FDR = 0.022), and PPP1R2 (OR = 1.09, CI = 1.05-1.13, FDR = 0.049). Findings from this TWAS enhance our understanding of how genetically-regulated gene expression affects bladder cancer development and point to potential prevention and treatment targets.

Gene regulation by convergent promoters

Convergent transcription, that is, the collision of sense and antisense transcription, is ubiquitous in mammalian genomes and believed to diminish RNA expression. Recently, antisense transcription downstream of promoters was found to be surprisingly prevalent. However, functional characteristics of affected promoters are poorly investigated. Here we show that convergent transcription marks an unexpected positively co-regulated promoter constellation. By assessing transcriptional dynamic systems, we identified co-regulated constituent promoters connected through a distinct chromatin structure. Within these cis-regulatory domains, transcription factors can regulate both constituting promoters by binding to only one of them. Convergent promoters comprise about a quarter of all active transcript start sites and initiate 5'-overlapping antisense RNAs-an RNA class believed previously to be rare. Visualization of nascent RNA molecules reveals convergent cotranscription at these loci. Together, our results demonstrate that co-regulated convergent promoters substantially expand the cis-regulatory repertoire, reveal limitations of the transcription interference model and call for adjusting the promoter concept.

Profiling Genome-Wide Methylation Patterns in Cattle Infected with Ostertagia ostertagi

DNA methylation (DNAm) regulates gene expression and genomic imprinting. This study aimed to investigate the effect of gastrointestinal (GI) nematode infection on host DNAm. Helminth-free Holstein steers were either infected with Ostertagia ostertagi (the brown stomach worm) or given tap water only as a control. Animals were euthanized 30 days post-infection, and tissues were collected at necropsy. We conducted epigenome-wide profiling using a mammalian methylation array to explore the impact of infection on methylation patterns in the mucosa from abomasal fundus (FUN), pylorus (PYL), draining lymph nodes (dLNs), and the duodenum (DUO). The analysis covered 31,107 cattle CpGs of 5082 genes and revealed infection-driven, tissue-specific, differential methylation patterns. A total of 389 shared and 2770 tissue-specific, differentially methylated positions (DMPs) were identified in dLN and FUN, particularly in genes associated with immune responses. The shared DMPs were found in 263 genes, many of which are involved in immune responses. Furthermore, 282, 244, 52, and 24 differentially methylated regions (DMRs) were observed in dLN, FUN, PYL, and DUO, respectively. More hypomethylated DMRs were detected in dLN and FUN, while more hypermethylated DMRs were found in PYL and DUO. Genes carrying DMPs and DMRs and enriched pathways relating to immune functions/responses were detected in infected animals, indicating a link between DNA methylation and the infection. The data may implicate a crucial role of DNAm in regulating the nature/strength of host immunity to infection and contribute to a deeper understanding of the epigenetic regulatory landscape in cattle infected by GI nematodes.

DDX39A resolves replication fork-associated RNA-DNA hybrids to balance fork protection and cleavage for genomic stability maintenance

Safeguarding replication fork stability in transcriptionally active regions is crucial for precise DNA replication and mutation prevention. Here, we discover the pervasive existence of replication fork-associated RNA-DNA hybrids (RF-RDs) in transcriptionally active regions of human cells. These hybrids function as protective barriers, preventing DNA2-mediated nascent DNA degradation and replication fork collapse under replication stress. We also identify DDX39A as a RAD51-associated protein that binds to stalled forks and resolves RF-RDs, facilitating proper DNA2-mediated DNA resection and replication fork restart. Excessive dissolution of RF-RDs causes replication fork collapse and genomic instability, while insufficient dissolution of RF-RDs under replication stress increases fork stability, resulting in chemoresistance that can be reversed by eliminating RF-RDs. In summary, we elucidated the prevalence of RF-RDs at replication forks within transcriptionally active regions, revealed their pivotal role in safeguarding replication fork stability, and proposed that targeting RF-RDs holds promise for augmenting chemotherapeutic efficacy.

Phylogeny, antiquity, and niche occupancy of Trinomia (Hymenoptera: Halictidae), an Afrotropical endemic genus of Nomiinae

The Afrotropical region is home to many endemic bee groups, yet almost none have been studied from an integrated, holistic perspective. Among them, the halictid subfamily Nomiinae contains exceptional African diversity with variable distributions and life histories. Here, we combine phylogenomics, molecular dating, and distributional modelling to explore the evolutionary ecology of the genus Trinomia. We analyzed a matrix of 59 species of Nomiinae using ultraconserved element (UCE) and whole genome data, including all six species of Trinomia, and estimated divergence times for the subfamily. We then generated distribution models for all six species of Trinomia using Maximum Entropy models (MaxEnt) and 671 spatial data points. From these methods, we discovered a monophyletic Trinomia with an unexpected sister group relationship to the Asian-endemic genus Gnathonomia, as well as a recent origin of Trinomia in the late Miocene (∼5.8 million years ago). From our results, we found hints of phylogenetic conservatism in distribution among sister-groups of Trinomia, however, our results also highlight the need for additional efforts inventorying, identifying, and sharing data on African bees. This study represents an exemplary first step into studying bee spatial phylogenomics of African endemic bees.

Uncovering the complexity of structural variants in four individuals with autism spectrum disorder

Autism spectrum disorder (ASD) is an increasingly recognized childhood developmental disorder. Despite extensive study, causal variants and molecular diagnosis remain elusive. There is both heterogeneity of the phenotype, as well as the genetic landscape associated with phenotype, which includes both inherited and de novo mutations. Currently, diagnosis is complex and behaviourally based, oftentimes occurring years after the ideal 1-2 years of age. Structural variants (SVs) are large and sometimes complex genomic variants that are likely underrepresented contributors to ASD due to the limitations of short-read DNA sequencing, such as alignment in repetitive regions and regions with GC bias. Here, we performed long-read sequencing (LRS) on four individuals with autism spectrum disorder to delineate SV complexity and determine precise breakpoints for SVs, which was not possible with short-read whole-genome sequencing (SRS). We use LRS to interrogate the methylation pattern associated with the SVs and phase the SV haplotypes to further clarify their contribution to disorder. LRS allows insight into the genome and methylome that allow us to uncover variant complexity and contribution that was previously unseen with SRS. Ultimately, this furthers precision diagnosis and contributes to individualized treatment for affected individuals and their families within the clinic.

Pharmacogenetic insights into ABCB1, ABCC2, CYP1A2, and CYP2B6 variants with epilepsy susceptibility among Egyptian Children: A retrospective case-control study

BACKGROUND

Pediatric epilepsy is a complicated neuropsychiatric disorder that is characterized by recurrent seizures and unusual synchronized electrical activities within brain tissues. It has a substantial effect on the quality of life of children, thus understanding of the hereditary considerations influencing epilepsy susceptibility and the response to antiepileptic medications is crucial. This study focuses on assessing the correlation of the ABCB1, ABCC2, CYP1A2, and CYP2B6 genetic polymorphisms with the susceptibility to epileptic seizures and their contributions to antiepileptic medication throughout the course of the disease.

METHODS

This study included 134 Egyptian epileptic children, comprising 67 drug-responsive and 67 drug-resistant patients, along with 124 healthy controls matching for age, gender, and geographical district. Genotyping of the rs2032582, rs717620, rs2273697, rs762551, and rs3745274 variants was performed using the PCR technique. Statistical analyses, including haplotype, multivariate, logistic regression, and bioinformatics approaches, were conducted to evaluate the associations within the disease.

RESULTS

The ABCC2*rs717620 (T allele) revealed an increased risk of epilepsy compared to healthy controls (OR = 2.12, p-value < 0.001), with the rs717620 (C/T + T/T genotypes) showing significant differences between drug-responsive and drug-resistant patients (p-value < 0.05). Moreover, the ABCC2*rs2273697 (A allele) indicated a decreased risk of epileptic seizures compared to healthy controls (OR = 0.51, p-value = 0.033), with the rs2273697 (G/A + A/A genotypes) indicating a significant association with drug-resistant patients (OR = 0.21, p-value = 0.002). The rs717620*T/rs2273697*G haplotype was significantly correlated with an elevated risk of epileptic seizures within drug-responsive patients (OR = 2.26, p-value = 0.019). Additionally, the CYP1A2*rs762551 (A allele) represented a protective effect against epilepsy susceptibility (OR = 0.50, p-value < 0.001), with the rs762551 (G/A + A/A genotypes) disclosing a substantial association with a decreased risk of epileptic seizures among drug-resistant patients compared to drug-responsive patients (OR = 0.07, p-value < 0.001). Conversely, the ABCB1*rs2032582 (G allele) and the CYP2B6*rs3745274 (T allele) did not attain a significant difference with the epilepsy risk compared to healthy controls (p-value > 0.05).

CONCLUSIONS

The findings of our study emphasize the importance of pharmacogenetic screening in epilepsy research, particularly regarding to drug-resistant patients. The ABCC2*rs717620 variant conferred a significant correlation with elevated risk of epileptic seizures, while the ABCC2*rs2273697 and CYP1A2*rs762551 variants confirmed their contributions as protective markers against epilepsy development. Conversely, the ABCB1*rs2032582 and CYP2B6*rs3745274 alleles were not considered as independent risk factors with the course of epilepsy disease.

Non-coding genetic variants underlying higher prostate cancer risk in men of African ancestry

Incidence and severity of prostate cancer (PrCa) substantially varies across ancestries. American men of African ancestry (AA) are more likely to be diagnosed with and die from PrCa than the those of European ancestry (EA). Published polygenic risk scores for developing prostate cancer, even those based on multi-ancestry genome-wide association studies, do not address population-specific genetic mechanisms underlying PrCa risk in men of African ancestry. Specifically, the role of non-coding regulatory polymorphisms in driving inter-ancestry variation in PrCa has not been suficiently explored. Here, by employing a sequence-based deep learning model of prostate regulatory enhancers, we identified ~ 2,000 SNPs with higher alternate allele frequency in AA men that potentially affect enhancer function associated with PrCa susceptibility, as supported by our experimental validation. The identified enhancer SNPs (eSNPs) may influence PrCa development through two complementary mechanisms: 1) the alternate allele that increase enhancer activity result in immune suppression and telomere elongation, and 2) the alternate alleles that decrease enhancer activity, lead to de-differentiation and inhibition of apoptosis. Notably, the eSNPs tend to disrupt the binding of known prostate transcription factors including FOX, AR and HOX families. Lastly, the identified eSNPs can be combined into a polygenic risk score that adds value to current GWAS-based risk variants in assessing PrCa risk in independent cohorts.

Comprehensive multi-tissue epigenome atlas in sheep: A resource for complex traits, domestication, and breeding

Comprehensive functional genome annotation is crucial to elucidate the molecular mechanisms of agronomic traits in livestock, yet systematic functional annotation of the sheep genome is lacking. Here, we generated 92 transcriptomic and epigenomic data sets from nine major tissues, along with whole-genome data from 2357 individuals across 29 breeds worldwide, and 4006 phenotypic data related to tail fat weight. We constructed the first multi-tissue epigenome atlas in terms of functional elements, chromatin states, and their functions and explored the utility of the functional elements in interpreting phenotypic variation during sheep domestication and improvement. Particularly, we identified a total of 753,723 nonredundant functional elements, with over 60% being novel. We found tissue-specific promoters and enhancers related to sensory abilities and immune response that were highly enriched in genomic regions influenced by domestication, while longissimus dorsi tissue-specific active enhancers and tail fat tissue-specific active promoters were highly enriched in genomic regions influenced by breeding and improvement. Notably, a variant, Chr13:51760995A>C, located in an enhancer region, was identified as a causal variant for tail fat deposition based on multi-layered data sets. Overall, this research provides foundational resources and a successful case for future investigations of complex traits in sheep through the integration of multi-omics data sets.

Patterns of chromosome evolution in ruminants

Studying when and where gross genomic rearrangements occurred during evolution is key to understanding changes in genome structure with functional consequences that might eventually lead to speciation. Here we identified chromosome rearrangements in ruminants, a clade characterized by large chromosome differences. Using 26 genome assemblies, we reconstructed five ancestral karyotypes and classified the rearrangement events occurring in each lineage. With these reconstructions, we then identified evolutionary breakpoints regions (EBRs) and synteny fragments. Ruminant karyotype evolution is characterized by inversions, while interchromosomal rearrangements occurred preferentially in the oldest ancestor of ruminants. We found that EBRs are depleted of protein coding genes, including housekeeping genes. Similarly, EBRs are not enriched in high GC regions, suggesting that meiotic double strand breaks might not be their origin. Overall, our results characterize at fine detail the location of chromosome rearrangements in ruminant evolution and provide new insights into the formation of EBRs.

Genetic linkage disequilibrium of deleterious mutations in threatened mammals

The impact of negative selection against deleterious mutations in endangered species remains underexplored. Recent studies have measured mutation load by comparing the accumulation of deleterious mutations, however, this method is most effective when comparing within and between populations of phylogenetically closely related species. Here, we introduced new statistics, LDcor, and its standardized form nLDcor, which allows us to detect and compare global linkage disequilibrium of deleterious mutations across species using unphased genotypes. These statistics measure averaged pairwise standardized covariance and standardize mutation differences based on the standard deviation of alleles to reflect selection intensity. We then examined selection strength in the genomes of seven mammals. Tigers exhibited an over-dispersion of deleterious mutations, while gorillas, giant pandas, and golden snub-nosed monkeys displayed negative linkage disequilibrium. Furthermore, the distribution of deleterious mutations in threatened mammals did not reveal consistent trends. Our results indicate that these newly developed statistics could help us understand the genetic burden of threatened species.

Genetic Architecture of Immune Cell DNA Methylation in the Rhesus Macaque

Genetic variation that impacts gene regulation, rather than protein function, can have strong effects on trait variation both within and between species. Epigenetic mechanisms, such as DNA methylation, are often an important intermediate link between genotype and phenotype, yet genetic effects on DNA methylation remain understudied in natural populations. To address this gap, we used reduced representation bisulfite sequencing to measure DNA methylation levels at 555,856 CpGs in peripheral whole blood of 573 samples collected from free-ranging rhesus macaques (Macaca mulatta) living on the island of Cayo Santiago, Puerto Rico. We used allele-specific methods to map cis-methylation quantitative trait loci (meQTL) and tested for effects of 243,389 single nucleotide polymorphisms (SNPs) on local DNA methylation levels. Of 776,092 tested SNP-CpG pairs, we identified 516,213 meQTL, with 69.12% of CpGs having at least one meQTL (FDR < 5%). On average, meQTL explained 21.2% of nearby methylation variance, significantly more than age or sex. meQTL were enriched in genomic compartments where methylation is likely to impact gene expression, for example, promoters, enhancers and binding sites for methylation-sensitive transcription factors. In support, using mRNA-seq data from 172 samples, we confirmed 332 meQTL as whole blood cis-expression QTL (eQTL) in the population, and found meQTL-eQTL genes were enriched for immune response functions, like antigen presentation and inflammation. Overall, our study takes an important step towards understanding the genetic architecture of DNA methylation in natural populations, and more generally points to the biological mechanisms driving phenotypic variation in our close relatives.

Chromosome-level genome assemblies and genetic maps reveal heterochiasmy and macrosynteny in endangered Atlantic Acropora

BACKGROUND

Over their evolutionary history, corals have adapted to sea level rise and increasing ocean temperatures, however, it is unclear how quickly they may respond to rapid change. Genome structure and genetic diversity contained within may highlight their adaptive potential.

RESULTS

We present chromosome-scale genome assemblies and linkage maps of the critically endangered Atlantic acroporids, Acropora palmata and A. cervicornis. Both assemblies and linkage maps were resolved into 14 chromosomes with their gene content and colinearity. Repeats and chromosome arrangements were largely preserved between the species. The family Acroporidae and the genus Acropora exhibited many phylogenetically significant gene family expansions. Macrosynteny decreased with phylogenetic distance. Nevertheless, scleractinians shared six of the 21 cnidarian ancestral linkage groups as well as numerous fission and fusion events compared to other distantly related cnidarians. Genetic linkage maps were constructed from one A. palmata family and 16 A. cervicornis families using a genotyping array. The consensus maps span 1,013.42 cM and 927.36 cM for A. palmata and A. cervicornis, respectively. Both species exhibited high genome-wide recombination rates (3.04 to 3.53 cM/Mb) and pronounced sex-based differences, known as heterochiasmy, with 2 to 2.5X higher recombination rates estimated in the female maps.

CONCLUSIONS

Together, the chromosome-scale assemblies and genetic maps we present here are the first detailed look at the genomic landscapes of the critically endangered Atlantic acroporids. These data sets revealed that adaptive capacity of Atlantic acroporids is not limited by their recombination rates. The sister species maintain macrosynteny with few genes with high sequence divergence that may act as reproductive barriers between them. In the Atlantic Acropora, hybridization between the two sister species yields an F1 hybrid with limited fertility despite the high levels of macrosynteny and gene colinearity of their genomes. Together, these resources now enable genome-wide association studies and discovery of quantitative trait loci, two tools that can aid in the conservation of these species.

Detection and analysis of complex structural variation in human genomes across populations and in brains of donors with psychiatric disorders

Complex structural variations (cxSVs) are often overlooked in genome analyses due to detection challenges. We developed ARC-SV, a probabilistic and machine-learning-based method that enables accurate detection and reconstruction of cxSVs from standard datasets. By applying ARC-SV across 4,262 genomes representing all continental populations, we identified cxSVs as a significant source of natural human genetic variation. Rare cxSVs have a propensity to occur in neural genes and loci that underwent rapid human-specific evolution, including those regulating corticogenesis. By performing single-nucleus multiomics in postmortem brains, we discovered cxSVs associated with differential gene expression and chromatin accessibility across various brain regions and cell types. Additionally, cxSVs detected in brains of psychiatric cases are enriched for linkage with psychiatric GWAS risk alleles detected in the same brains. Furthermore, our analysis revealed significantly decreased brain-region- and cell-type-specific expression of cxSV genes, specifically for psychiatric cases, implicating cxSVs in the molecular etiology of major neuropsychiatric disorders.

Comparative single-cell analyses identify shared and divergent features of human and mouse kidney development

The mammalian kidney maintains fluid homeostasis through diverse epithelial cell types generated from nephron and ureteric progenitor cells. To extend a developmental understanding of the kidney's epithelial networks, we compared chromatin organization (single-nuclear assay for transposase-accessible chromatin sequencing [ATAC-seq]; 112,864 nuclei) and gene expression (single-cell/nuclear RNA sequencing [RNA-seq]; 109,477 cells/nuclei) in the developing human (10.6-17.6 weeks; n = 10) and mouse (post-natal day [P]0; n = 10) kidney, supplementing analysis with published mouse datasets from earlier stages. Single-cell/nuclear datasets were analyzed at a species level, and then nephron and ureteric cellular lineages were extracted and integrated into a common, cross-species, multimodal dataset. Comparative computational analyses identified conserved and divergent features of chromatin organization and linked gene activity, identifying species-specific and cell-type-specific regulatory programs. In situ validation of human-enriched gene activity points to human-specific signaling interactions in kidney development. Further, human-specific enhancer regions were linked to kidney diseases through genome-wide association studies (GWASs), highlighting the potential for clinical insight from developmental modeling.

StableLift: Optimized Germline and Somatic Variant Detection Across Genome Builds

Reference genomes are foundational to modern genomics. Our growing understanding of genome structure leads to continual improvements in reference genomes and new genome "builds" with incompatible coordinate systems. We quantified the impact of genome build on germline and somatic variant calling by analyzing tumour-normal whole-genome pairs against the two most widely used human genome builds. The average individual had a build-discordance of 3.8% for germline SNPs, 8.6% for germline SVs, 25.9% for somatic SNVs and 49.6% for somatic SVs. Build-discordant variants are not simply false-positives: 47% were verified by targeted resequencing. Build-discordant variants were associated with specific genomic and technical features in variant- and algorithm-specific patterns. We leveraged these patterns to create StableLift, an algorithm that predicts cross-build stability with AUROCs of 0.934 ± 0.029. These results call for significant caution in cross-build analyses and for use of StableLift as a computationally efficient solution to mitigate inter-build artifacts.

Mapping the epigenomic landscape of post-traumatic stress disorder in human cortical neurons

The study conducted a comprehensive genome-wide analysis of differential 5mC and 5hmC modifications at both CpG and non-CpG sites in postmortem orbitofrontal neurons from 25 PTSD cases and 13 healthy controls. It was observed that PTSD patients exhibit a greater number of differential 5hmC sites compared to 5mC sites. Specifically, individuals with PTSD tend to show hyper-5mC/5hmC at CpG sites, particularly within CpG islands and promoter regions, and hypo-5mC/5hmC at non-CpG sites, especially within intragenic regions. Functional enrichment analysis indicated distinct yet interconnected roles for 5mC and 5hmC in PTSD. The 5mC marks primarily regulate cell-cell adhesion processes, whereas 5hmC marks are involved in embryonic morphogenesis and cell fate commitment. By integrating published PTSD findings from central and peripheral tissues through multi-omics approaches, several biological mechanisms were prioritized, including developmental processes, HPA axis regulation, and immune responses. Based on the consistent enrichment in developmental processes, we hypothesize that if epigenetic changes occur during early developmental stages, they may increase the risk of developing PTSD following trauma exposure. Conversely, if these epigenetic changes occur in adulthood, they may influence neuronal apoptosis and survival mechanisms.

Genomic adaptation to small population size and saltwater consumption in the critically endangered Cat Ba langur

Many mammal species have declining populations, but the consequences of small population size on the genomic makeup of species remain largely unknown. We investigated the evolutionary history, genetic load and adaptive potential of the Cat Ba langur (Trachypithecus poliocephalus), a primate species endemic to Vietnam's famous Ha Long Bay and with less than 100 living individuals one of the most threatened primates in the world. Using high-coverage whole genome data of four wild individuals, we revealed the Cat Ba langur as sister species to its conspecifics of the northern limestone langur clade and found no evidence for extensive secondary gene flow after their initial separation. Compared to other primates and mammals, the Cat Ba langur showed low levels of genetic diversity, long runs of homozygosity, high levels of inbreeding and an excess of deleterious mutations in homozygous state. On the other hand, genetic diversity has been maintained in protein-coding genes and on the gene-rich human chromosome 19 ortholog, suggesting that the Cat Ba langur retained most of its adaptive potential. The Cat Ba langur also exhibits several unique non-synonymous variants that are related to calcium and sodium metabolism, which may have improved adaptation to high calcium intake and saltwater consumption.

Genetic Regulatory Perturbation of Gene Expression Impacted by Genomic Introgression in Fiber Development of Allotetraploid Cotton

Interspecific genomic introgression is an important evolutionary process with respect to the generation of novel phenotypic diversity and adaptation. A key question is how gene flow perturbs gene expression networks and regulatory interactions. Here, an introgression population of two species of allopolyploid cotton (Gossypium) to delineate the regulatory perturbations of gene expression regarding fiber development accompanying fiber quality change is utilized. De novo assembly of the recipient parent (G. hirsutum Emian22) genome allowed the identification of genomic variation and introgression segments (ISs) in 323 introgression lines (ILs) from the donor parent (G. barbadense 3-79). It documented gene expression dynamics by sequencing 1,284 transcriptomes of developing fibers and characterized genetic regulatory perturbations mediated by genomic introgression using a multi-locus model. Introgression of individual homoeologous genes exhibiting extreme low or high expression bias can lead to a parallel expression bias in their non-introgressed duplicates, implying a shared yet divergent regulatory fate of duplicated genes following allopolyploidy. Additionally, the IL N182 with improved fiber quality is characterized, and the candidate gene GhFLAP1 related to fiber length is validated. This study outlines a framework for understanding introgression-mediated regulatory perturbations in polyploids, and provides insights for targeted breeding of superior upland cotton fiber.

Identification and analysis of short indels inducing exon extension/shrinkage events

The search for genetic variants that act as causative factors in human diseases by disrupting the normal splicing process has primarily focused on single nucleotide variants (SNVs). It is worth noting that insertions or deletions (indels) have also been sporadically reported as causative disease variants through their potential impact on the splicing process. In this study, to perform identification of indels inducing exon extension/shrinkage events, we used individual-specific genomes and RNA sequencing (RNA-seq) data pertaining to the corresponding individuals and identified 12 exon extension/shrinkage events that were potentially induced by indels that disrupted authentic splice sites or created novel splice sites in 235 normal individuals. By evaluating the impact of these abnormal splicing events on the resulting transcripts, we found that five events led to the generation of premature termination codons (PTCs), including those occurring within genes associated with genetic disorders. Our analysis revealed that the potential functions of indels have been underexamined, and it is worth considering the possibility that indels may affect splice site usage, using RNA-seq data to discover novel potentially disease-associated mutations.

Identification of the VLDLR locus associated with giant cell arteritis and the possible causal role of low-density lipoprotein cholesterol in its pathogenesis

OBJECTIVES

To elucidate the association between genetic variants and the risk of GCA via large-scale genome-wide association studies (GWAS). In addition, to assess the causal effect of a specific molecule by employing the obtained GWAS results as genetic epidemiological tools.

METHODS

We applied additional variant quality control to the publicly available GWAS results from the biobanks of the UK (UKBB) and Finland (FinnGen), which comprised 532 cases vs 408 565 controls and 884 cases vs 332 115 controls, respectively. We further meta-analysed these two sets of results. We performed two-sample Mendelian randomization (MR) to test the causal effect of low-density lipoprotein (LDL) cholesterol on the risk of GCA.

RESULTS

The MHC class II region showed significant associations in UKBB, FinnGen and the meta-analysis. The VLDLR region was associated with GCA risk in the meta-analysis. The T allele of rs7044155 increased the expression of VLDLR, decreased the LDL cholesterol level and decreased the disease risk. The subsequent MR results indicated that a 1 s.d. increase in LDL cholesterol was associated with an increased risk of GCA (odds ratio 1.21, 95% CI 1.01-1.45; P = 0.04).

CONCLUSIONS

Our study identified associations between GCA risk and the MHC class II and VLDLR regions. Moreover, LDL cholesterol was suggested to have a causal effect on the risk of developing GCA.

p53motifDB: integration of genomic information and tumor suppressor p53 binding motifs

The tumor suppressor gene TP53 encodes the DNA binding transcription factor p53 and is one of the most commonly mutated genes in human cancer. Tumor suppressor activity requires binding of p53 to its DNA response elements and subsequent transcriptional activation of a diverse set of target genes. Despite decades of close study, the logic underlying p53 interactions with its numerous potential genomic binding sites and target genes is not yet fully understood. Here, we present a database of DNA and chromatin-based information focused on putative p53 binding sites in the human genome to allow users to generate and test new hypotheses related to p53 activity in the genome. Users can query genomic locations based on experimentally observed p53 binding, regulatory element activity, genetic variation, evolutionary conservation, chromatin modification state, and chromatin structure. We present multiple use cases demonstrating the utility of this database for generating novel biological hypotheses, such as chromatin-based determinants of p53 binding and potential cell type-specific p53 activity. All database information is also available as a precompiled sqlite database for use in local analysis or as a Shiny web application.

Epigenetic programming of host lipid metabolism associated with resistance to TST/IGRA conversion after exposure to Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) exposure leads to a range of outcomes including clearance, latent TB infection (LTBI), and pulmonary tuberculosis (TB). Some heavily exposed individuals resist tuberculin skin test (TST) and interferon-gamma (IFNγ) release assay (IGRA) conversion (RSTR), which suggests that they employ IFNγ-independent mechanisms of Mtb control. Here, we compare monocyte epigenetic profiles of RSTR and LTBI from a Ugandan household contact cohort. Chromatin accessibility did not differ between uninfected RSTR and LTBI monocytes. By contrast, methylation significantly differed at 174 CpG sites and across 63 genomic regions. Consistent with previous transcriptional findings in this cohort, differential methylation was enriched in lipid- and cholesterol-associated pathways including the genes APOC3, KCNQ1, and PLA2G3. In addition, methylation was enriched in Hippo signaling, which is associated with cholesterol homeostasis and includes CIT and SHANK2. Lipid export and Hippo signaling pathways were also associated with gene expression in response to Mtb in RSTR as well as IFN stimulation in monocyte-derived macrophages (MDMs) from an independent healthy donor cohort. Moreover, serum-derived high-density lipoprotein from RSTR had elevated ABCA1-mediated cholesterol efflux capacity (CEC) compared to LTBI. Our findings suggest that resistance to TST/IGRA conversion is linked to regulation of lipid accumulation in monocytes, which could facilitate early Mtb clearance among RSTR subjects through IFNγ-independent mechanisms.IMPORTANCETuberculosis (TB) remains an enduring global health challenge with millions of deaths and new cases each year. Despite recent advances in TB treatment, we lack an effective vaccine or a durable cure. While heavy exposure to Mycobacterium tuberculosis often results in latent TB latent infection (LTBI), subpopulations exist that are either resistant to infection or contain Mtb with interferon-gamma (IFNγ)-independent mechanisms not indicative of LTBI. These resisters provide an opportunity to investigate the mechanisms of TB disease and discover novel therapeutic targets. Here, we compare monocyte epigenetic profiles of RSTR and LTBI from a Ugandan household contact cohort. We identify methylation signatures in host lipid and cholesterol pathways with potential relevance to early TB clearance before the sustained IFN responses indicative of LTBI. This adds to a growing body of literature linking TB disease outcomes to host lipids.

primerForge: a Python program for identifying primer pairs capable of distinguishing groups of genomes from each other

In both molecular epidemiology and microbial ecology, it is useful to be able to categorize specific strains of microorganisms in either an ingroup or an outgroup in a given population, e.g. to distinguish a pathogenic strain of interest from its non-virulent relatives. An "ingroup" refers to a group of microbes that are the primary focus of study or interest. Conversely, an "outgroup" consists of microbes that are closely-related to, but have evolved separately from, the ingroup. While whole genome sequencing and downstream phylogenetic analyses can be employed to do this, these techniques are often slow and can be resource intensive. Additionally, the laboratory would have to sequence the whole genome to use these tools to determine whether or not a new sample is part of the ingroup or outgroup. Alternatively, polymerase chain reaction (PCR) can be used to amplify regions of genetic material that are specific to the strain(s) of interest. PCR is faster, less expensive, and more accessible than whole genome sequencing, so having a PCR-based approach can accelerate the detection of specific strain(s) of microbes and facilitate diagnoses and/or population studies.

Transcriptional function of E2A, Ebf1, Pax5, Ikaros and Aiolos analyzed by in vivo acute protein degradation in early B cell development

Early B cell lymphopoiesis depends on E2A, Ebf1, Pax5 and Ikaros family members. In the present study, we used acute protein degradation in mice to identify direct target genes of these transcription factors in pro-B, small pre-B and immature B cells. E2A, Ebf1 and Pax5 predominantly function as transcriptional activators by inducing open chromatin at their target genes, have largely unique functions and are essential for early B cell maintenance. Ikaros and Aiolos act as dedicated repressors to cooperatively control early B cell development. The surrogate light-chain genes Igll1 and Vpreb1 are directly activated by Ebf1 and Pax5 in pro-B cells and directly repressed by Ikaros and Aiolos in small pre-B cells. Pax5 and E2A contribute to V(D)J recombination by activating Rag1, Rag2, Dntt, Irf4 and Irf8. Similar to Pax5, Ebf1 also represses the cohesin-release factor gene Wapl to mediate prolonged loop extrusion across the Igh locus. In summary, in vivo protein degradation has provided unprecedented insight into the control of early B cell lymphopoiesis by five transcription factors.

Clusters of mammalian conserved RNA structures in UTRs associate with RBP binding sites

RNA secondary structures play essential roles in the formation of the tertiary structure and function of a transcript. Recent genome-wide studies highlight significant potential for RNA structures in the mammalian genome. However, a major challenge is assigning functional roles to these structured RNAs. In this study, we conduct a guilt-by-association analysis of clusters of computationally predicted conserved RNA structure (CRSs) in human untranslated regions (UTRs) to associate them with gene functions. We filtered a broad pool of ∼500 000 human CRSs for UTR overlap, resulting in 4734 and 24 754 CRSs from the 5' and 3' UTR of protein-coding genes, respectively. We separately clustered these CRSs for both sets using RNAscClust, obtaining 793 and 2403 clusters, each containing an average of five CRSs per cluster. We identified overrepresented binding sites for 60 and 43 RNA-binding proteins co-localizing with the clustered CRSs. Furthermore, 104 and 441 clusters from the 5' and 3' UTRs, respectively, showed enrichment for various Gene Ontologies, including biological processes such as 'signal transduction', 'nervous system development', molecular functions like 'transferase activity' and the cellular components such as 'synapse' among others. Our study shows that significant functional insights can be gained by clustering RNA structures based on their structural characteristics.

Regulatory transposable elements in the encyclopedia of DNA elements

Transposable elements (TEs) comprise ~50% of our genome, but knowledge of how TEs affect genome evolution remains incomplete. Leveraging ENCODE4 data, we provide the most comprehensive study to date of TE contributions to the regulatory genome. We find 236,181 (~25%) human candidate cis-regulatory elements (cCREs) are TE-derived, with over 90% lineage-specific since the human-mouse split, accounting for 8-36% of lineage-specific cCREs. Except for SINEs, cCRE-associated transcription factor (TF) motifs in TEs are derived from ancestral TE sequence more than expected by chance. We show that TEs may adopt similar regulatory activities of elements near their integration site. Since human-mouse divergence, TEs have contributed 3-56% of TF binding site turnover events across 30 examined TFs. Finally, TE-derived cCREs are similar to non-TE cCREs in terms of MPRA activity and GWAS variant enrichment. Overall, our results substantiate the notion that TEs have played an important role in shaping the human regulatory genome.

Systematic discovery of gene-environment interactions underlying the human plasma proteome in UK Biobank

Understanding how gene-environment interactions (GEIs) influence the circulating proteome could aid in biomarker discovery and validation. The presence of GEIs can be inferred from single nucleotide polymorphisms that associate with phenotypic variability - termed variance quantitative trait loci (vQTLs). Here, vQTL association studies are performed on plasma levels of 1463 proteins in 52,363 UK Biobank participants. A set of 677 independent vQTLs are identified across 568 proteins. They include 67 variants that lack conventional additive main effects on protein levels. Over 1100 GEIs are identified between 101 proteins and 153 environmental exposures. GEI analyses uncover possible mechanisms that explain why 13/67 vQTL-only sites lack corresponding main effects. Additional analyses also highlight how age, sex, epistatic interactions and statistical artefacts may underscore associations between genetic variation and variance heterogeneity. This study establishes the most comprehensive database yet of vQTLs and GEIs for the human proteome.

PHF6 cooperates with SWI/SNF complexes to facilitate transcriptional progression

Genes encoding subunits of SWI/SNF (BAF) chromatin remodeling complexes are mutated in nearly 25% of cancers. To gain insight into the mechanisms by which SWI/SNF mutations drive cancer, we contributed ten rhabdoid tumor (RT) cell lines mutant for SWI/SNF subunit SMARCB1 to a genome-scale CRISPR-Cas9 depletion screen performed across 896 cell lines. We identify PHF6 as specifically essential for RT cell survival and demonstrate that dependency on Phf6 extends to Smarcb1-deficient cancers in vivo. As mutations in either SWI/SNF or PHF6 can cause the neurodevelopmental disorder Coffin-Siris syndrome, our findings of a dependency suggest a previously unrecognized functional link. We demonstrate that PHF6 co-localizes with SWI/SNF complexes at promoters, where it is essential for maintenance of an active chromatin state. We show that in the absence of SMARCB1, PHF6 loss disrupts the recruitment and stability of residual SWI/SNF complex members, collectively resulting in the loss of active chromatin at promoters and stalling of RNA Polymerase II progression. Our work establishes a mechanistic basis for the shared syndromic features of SWI/SNF and PHF6 mutations in CSS and the basis for selective dependency on PHF6 in SMARCB1-mutant cancers.

Quantitative analysis of cis-regulatory elements in transcription with KAS-ATAC-seq

Cis-regulatory elements (CREs) are pivotal in orchestrating gene expression throughout diverse biological systems. Accurate identification and in-depth characterization of functional CREs are crucial for decoding gene regulation networks during cellular processes. In this study, we develop Kethoxal-Assisted Single-stranded DNA Assay for Transposase-Accessible Chromatin with Sequencing (KAS-ATAC-seq) to quantitatively analyze the transcriptional activity of CREs. A main advantage of KAS-ATAC-seq lies in its precise measurement of ssDNA levels within both proximal and distal ATAC-seq peaks, enabling the identification of transcriptional regulatory sequences. This feature is particularly adept at defining Single-Stranded Transcribing Enhancers (SSTEs). SSTEs are highly enriched with nascent RNAs and specific transcription factors (TFs) binding sites that define cellular identity. Moreover, KAS-ATAC-seq provides a detailed characterization and functional implications of various SSTE subtypes. Our analysis of CREs during mouse neural differentiation demonstrates that KAS-ATAC-seq can effectively identify immediate-early activated CREs in response to retinoic acid (RA) treatment. Our findings indicate that KAS-ATAC-seq provides more precise annotation of functional CREs in transcription. Future applications of KAS-ATAC-seq would help elucidate the intricate dynamics of gene regulation in diverse biological processes.

Population structure in Mixornis tit-babblers across Sunda Shelf matches interfluvia of paleo-rivers

Rivers constitute an important biogeographic divide in vast areas of tropical rainforest, such as the Amazon and Congo Basins. Southeast Asia's rainforests are currently fragmented across islands divided by sea, which has long obscured their extensive history of terrestrial connectivity as part of a vast (but now submerged) subcontinent - Sundaland - during most of the Quaternary. The role of paleo-rivers in determining population structure in Sundaic rainforests at a time when these forests were connected remains little understood. We examined the coloration of museum skins and used the genomic DNA of museum samples and freshly-collected blood tissue of a pair of Sundaic songbird species, the pin-striped and bold-striped tit-babblers (Mixornis gularis and M. bornensis, respectively), to assess the genetic affinity of populations on small Sundaic islands that have largely been ignored by modern research. Our genomic and morphological results place the populations from the Anambas and Natuna Islands firmly within M. gularis from the Malay Peninsula in western Sundaland, even though some of these islands are geographically much closer to Borneo, where M. bornensis resides. Our results reveal genetic structure consistent with the course of Sundaic paleo-rivers and the location of the interfluvia they formed, and add to a small but growing body of evidence that rivers would have been of equal biogeographic importance in Sundaland's former connected forest landscape as they are in Amazonia and the Congo Basin today.

Whole genome sequencing study of identical twins discordant for psychosis

Monozygotic (MZ) twins are often thought to have identical genomes, but recent work has shown that early post-zygotic events can result in a spectrum of DNA variants that are different between MZ twins. Such variants may explain phenotypic discordance and contribute to disease etiology. Here we performed whole genome sequencing in 17 pairs of MZ twins discordant for a psychotic disorder (schizophrenia, schizoaffective disorder or bipolar disorder). We examined various classes of rare variants that are discordant within a twin pair. We identified four genes harboring rare, predicted deleterious missense variants that were private to an affected individual in the cohort. Variants in FOXN1 and FLOT2 would have been categorized as damaging from recent schizophrenia and bipolar exome sequencing studies. Additionally, we identified four rare genic copy number variants (CNVs) private to an affected sample, two of which overlapped genes that have shown evidence for association with schizophrenia or bipolar disorder. One such CNV was a 3q29 duplication previously implicated in autism and developmental delay. We have performed the largest MZ twin study for discordant psychotic phenotypes to date. These findings warrant further investigation using other analytical approaches.

A multiomic atlas of the aging hippocampus reveals molecular changes in response to environmental enrichment

Aging involves the deterioration of organismal function, leading to the emergence of multiple pathologies. Environmental stimuli, including lifestyle, can influence the trajectory of this process and may be used as tools in the pursuit of healthy aging. To evaluate the role of epigenetic mechanisms in this context, we have generated bulk tissue and single cell multi-omic maps of the male mouse dorsal hippocampus in young and old animals exposed to environmental stimulation in the form of enriched environments. We present a molecular atlas of the aging process, highlighting two distinct axes, related to inflammation and to the dysregulation of mRNA metabolism, at the functional RNA and protein level. Additionally, we report the alteration of heterochromatin domains, including the loss of bivalent chromatin and the uncovering of a heterochromatin-switch phenomenon whereby constitutive heterochromatin loss is partially mitigated through gains in facultative heterochromatin. Notably, we observed the multi-omic reversal of a great number of aging-associated alterations in the context of environmental enrichment, which was particularly linked to glial and oligodendrocyte pathways. In conclusion, our work describes the epigenomic landscape of environmental stimulation in the context of aging and reveals how lifestyle intervention can lead to the multi-layered reversal of aging-associated decline.

Relationships between nine neuropsychiatric disorders and cervical cancer: insights from genetics, causality and shared gene expression patterns

BACKGROUND

Neuropsychiatric disorders and cervical cancer exert substantial influences on women's health. Furthermore, neuropsychiatric disorders frequently manifest as common symptoms in cancer patients, potentially increasing the risk of malignant neoplasms. This study aimed to identify neuropsychiatric disorders that are genetically and causally related to cervical cancer and to investigate the molecular mechanisms underlying these associations.

METHODS

GWAS data related to nine neuropsychiatric disorders, namely, schizophrenia, bipolar disorder, autism spectrum disorder, Parkinson's disease, anxiety, Alzheimer's disease, mood disorders, depression, and alcohol dependence, were obtained to calculate heritability (h2) and genetic correlation (rg) with cervical cancer using linkage disequilibrium score regression (LDSC). Mendelian randomization (MR) analysis of the two cohorts was employed to assess the causal effects. Shared gene expression pattern analysis was subsequently conducted to investigate the molecular mechanism underlying these significant associations.

RESULTS

Anxiety, mood disorders, depression, and alcohol dependence were genetically correlated with cervical cancer (all adjusted P < 0.05). Only depression was causally related to cervical cancer in both the discovery (ORIVW: 1.41, PIVW = 0.02) and replication cohorts (ORIVW: 1.80, PIVW = 0.03) in the MR analysis. Gene expression pattern analysis revealed that 270 genes related to depression and cervical cancer, including tumour necrosis factor (TNF), were significantly upregulated in cervical cancer patients, while vascular endothelial growth factor A (VEGFA), transcription factor AP-1 (JUN), and insulin-like growth factor I (IGF-I) were associated with prognosis in cervical cancer patients (all P < 0.05). These overlapping genes implicated the involvement of multiple biological mechanisms, such as neuron death, the PI3K-Akt signalling pathway, and human papillomavirus infection.

CONCLUSIONS

Genetic, causal and molecular evidence indicates that depression increases the risk of cervical cancer. The TNF, VEGFA, JUN, and IGF-1 genes and the neuron death, PI3K-Akt, and human papillomavirus infection signalling pathways may possibly explain this association.

Intratumoral immune triads are required for immunotherapy-mediated elimination of solid tumors

Tumor-specific CD8+ T cells are frequently dysfunctional and unable to halt tumor growth. We investigated whether tumor-specific CD4+ T cells can be enlisted to overcome CD8+ T cell dysfunction within tumors. We find that the spatial positioning and interactions of CD8+ and CD4+ T cells, but not their numbers, dictate anti-tumor responses in the context of adoptive T cell therapy as well as immune checkpoint blockade (ICB): CD4+ T cells must engage with CD8+ T cells on the same dendritic cell during the effector phase, forming a three-cell-type cluster (triad) to license CD8+ T cell cytotoxicity and cancer cell elimination. When intratumoral triad formation is disrupted, tumors progress despite equal numbers of tumor-specific CD8+ and CD4+ T cells. In patients with pleural mesothelioma treated with ICB, triads are associated with clinical responses. Thus, CD4+ T cells and triads are required for CD8+ T cell cytotoxicity during the effector phase and tumor elimination.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

Kataegis associated mutational processes linked to adverse prostate cancer presentation in African men

Kataegis, the focal hypermutation of single base substitutions (SBS) in tumour genomes, has received little attention with respect to prostate cancer (PCa) associated molecular and clinical features. Most notably, data is lacking with regards to this tumour evolutionary phenomenon and PCa racial disparities, with African men disproportionately impacted. Here through comparison between African (n = 109) and non-African (n = 79) whole genome sequenced treatment naïve primary tumours, using a single analytical workflow we assessed for shared and unique features of kataegis. Linking kataegis to aggressive presentation, structural variant burden and copy number loss, we attributed APOBEC3 activity through higher rates of SBS2 to high-risk African tumours. While kataegis positive African patients presented with elevated prostate specific antigen levels, their tumours showed evolutionary unique trajectories marked by increased subclonal and structural variant-independent kataegis. The potential to exacerbate tumour heterogeneity emphases the significance of continued exploration of biological behaviours and environmental exposures for African patients.

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae

BACKGROUND

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems.

RESULTS

We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary.

CONCLUSIONS

Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

Control of endothelial cell function and arteriogenesis by MEG3:EZH2 epigenetic regulation of integrin expression

Epigenetic processes involving long non-coding RNAs regulate endothelial gene expression. However, the underlying regulatory mechanisms causing endothelial dysfunction remain to be elucidated. Enhancer of zeste homolog 2 (EZH2) is an important rheostat of histone H3K27 trimethylation (H3K27me3) that represses endothelial targets, but EZH2 RNA binding capacity and EZH2:RNA functional interactions have not been explored in post-ischemic angiogenesis. We used formaldehyde/UV-assisted crosslinking ligation and sequencing of hybrids and identified a new role for maternally expressed gene 3 (MEG3). MEG3 formed the predominant RNA:RNA hybrid structures in endothelial cells. Moreover, MEG3:EZH2 assists recruitment onto chromatin. By EZH2-chromatin immunoprecipitation, following MEG3 depletion, we demonstrated that MEG3 controls recruitment of EZH2/H3K27me3 onto integrin subunit alpha4 (ITGA4) promoter. Both MEG3 knockdown or EZH2 inhibition (A-395) promoted ITGA4 expression and improved endothelial cell migration and adhesion to fibronectin in vitro. The A-395 inhibitor re-directed MEG3-assisted chromatin remodeling, offering a direct therapeutic benefit by increasing endothelial function and resilience. This approach subsequently increased the expression of ITGA4 in arterioles following ischemic injury in mice, thus promoting arteriogenesis. Our findings show a context-specific role for MEG3 in guiding EZH2 to repress ITGA4. Novel therapeutic strategies could antagonize MEG3:EZH2 interaction for pre-clinical studies.

Atypical Modes of CTCF Binding Facilitate Tissue-Specific and Neuronal Activity-Dependent Gene Expression States

Multivalent binding of CTCF to variable DNA sequences is thought to underlie its ability to mediate diverse cellular functions. CTCF typically binds a 20 base-pair consensus DNA sequence, but the full diversity of CTCF binding sites (CBS) within the genome has not been interrogated. We assessed CTCF occupancy in cultured cortical neurons and observed surprisingly that ~ 22% of CBS lack the consensus CTCF motif. We report here that sequence diversity at most of these atypical CBS involves degeneracy at specific nucleotide positions within the consensus CTCF motif, which likely affect the binding of CTCF zinc fingers 6 and 7. This mode of atypical CTCF binding defines most CBS at gene promoters, as well as CBS that are dynamically altered during neural differentiation and following neuronal stimulation, revealing how atypical CTCF binding could influence gene activity. Dynamic CBS are distributed both within and outside loop anchors and TAD boundaries, suggesting both looping-dependent and independent roles for CTCF. Finally, we describe a second mode of atypical CTCF binding to DNA sequences that are completely unrelated to the consensus CTCF motif, which are enriched within the bodies of tissue-specific genes. These tissue-specific atypical CBS are also enriched in H3K27ac, which marks cis-regulatory elements within chromatin, including enhancers. Overall, these results indicate how atypical CBS could dynamically regulate gene activity patterns during differentiation, development, and in response to environmental cues.

DNA mismatch and damage patterns revealed by single-molecule sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other diseases1,2. Most mutations begin as nucleotide mismatches or damage in one of the two strands of the DNA before becoming double-strand mutations if unrepaired or misrepaired3,4. However, current DNA-sequencing technologies cannot accurately resolve these initial single-strand events. Here we develop a single-molecule, long-read sequencing method (Hairpin Duplex Enhanced Fidelity sequencing (HiDEF-seq)) that achieves single-molecule fidelity for base substitutions when present in either one or both DNA strands. HiDEF-seq also detects cytosine deamination-a common type of DNA damage-with single-molecule fidelity. We profiled 134 samples from diverse tissues, including from individuals with cancer predisposition syndromes, and derive from them single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumours deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples that are deficient in only polymerase proofreading. We also define a single-strand damage signature for APOBEC3A. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. As double-strand DNA mutations are only the end point of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable studies of how mutations arise in a variety of contexts, especially in cancer and ageing.

ACMGA: a reference-free multiple-genome alignment pipeline for plant species

BACKGROUND

The short-read whole-genome sequencing (WGS) approach has been widely applied to investigate the genomic variation in the natural populations of many plant species. With the rapid advancements in long-read sequencing and genome assembly technologies, high-quality genome sequences are available for a group of varieties for many plant species. These genome sequences are expected to help researchers comprehensively investigate any type of genomic variants that are missed by the WGS technology. However, multiple genome alignment (MGA) tools designed by the human genome research community might be unsuitable for plant genomes.

RESULTS

To fill this gap, we developed the AnchorWave-Cactus Multiple Genome Alignment (ACMGA) pipeline, which improved the alignment of repeat elements and could identify long (> 50 bp) deletions or insertions (INDELs). We conducted MGA using ACMGA and Cactus for 8 Arabidopsis (Arabidopsis thaliana) and 26 Maize (Zea mays) de novo assembled genome sequences and compared them with the previously published short-read variant calling results. MGA identified more single nucleotide variants (SNVs) and long INDELs than did previously published WGS variant callings. Additionally, ACMGA detected significantly more SNVs and long INDELs in repetitive regions and the whole genome than did Cactus. Compared with the results of Cactus, the results of ACMGA were more similar to the previously published variants called using short-read. These two MGA pipelines identified numerous multi-allelic variants that were missed by the WGS variant calling pipeline.

CONCLUSIONS

Aligning de novo assembled genome sequences could identify more SNVs and INDELs than mapping short-read. ACMGA combines the advantages of AnchorWave and Cactus and offers a practical solution for plant MGA by integrating global alignment, a 2-piece-affine-gap cost strategy, and the progressive MGA algorithm.

A Polygenic Risk Analysis for Identifying Ulcerative Colitis Patients with European Ancestry

The incidence of ulcerative colitis (UC) has increased globally. As a complex disease, the genetic predisposition for UC could be estimated by the polygenic risk score (PRS), which aggregates the effects of a large number of genetic variants in a single quantity and shows promise in identifying individuals at higher lifetime risk of UC. Here, based on a cohort of 2869 UC cases and 2900 controls with genotype array datasets, we used PRSice-2 to calculate PRS, and systematically analyzed factors that could affect the power of PRS, including GWAS summary statistics, population stratification, and impact of variants. After leveraging a stepwise condition analysis, we eventually established the best PRS model, achieving an AUC of 0.713. Meanwhile, samples in the top 20% of the PRS distribution had a risk of UC more than ten times higher than samples in the lowest 20% (OR = 10.435, 95% CI 8.571-12.703). Our analyses demonstrated that including population-enriched, more disease-associated SNPs and using GWAS summary statistics from similar ethnic background can improve the power of PRS. Strictly following the principle of focusing on one population in all aspects of generating PRS can be a cost-effective way to apply genotype-array-derived PRS to practical risk estimation.

Breaking down causes, consequences, and mediating effects of telomere length variation on human health

BACKGROUND

Telomeres form repeated DNA sequences at the ends of chromosomes, which shorten with each cell division. Yet, factors modulating telomere attrition and the health consequences thereof are not fully understood. To address this, we leveraged data from 326,363 unrelated UK Biobank participants of European ancestry.

RESULTS

Using linear regression and bidirectional univariable and multivariable Mendelian randomization (MR), we elucidate the relationships between leukocyte telomere length (LTL) and 142 complex traits, including diseases, biomarkers, and lifestyle factors. We confirm that telomeres shorten with age and show a stronger decline in males than in females, with these factors contributing to the majority of the 5.4% of LTL variance explained by the phenome. MR reveals 23 traits modulating LTL. Smoking cessation and high educational attainment associate with longer LTL, while weekly alcohol intake, body mass index, urate levels, and female reproductive events, such as childbirth, associate with shorter LTL. We also identify 24 traits affected by LTL, with risk for cardiovascular, pulmonary, and some autoimmune diseases being increased by short LTL, while longer LTL increased risk for other autoimmune conditions and cancers. Through multivariable MR, we show that LTL may partially mediate the impact of educational attainment, body mass index, and female age at childbirth on proxied lifespan.

CONCLUSIONS

Our study sheds light on the modulators, consequences, and the mediatory role of telomeres, portraying an intricate relationship between LTL, diseases, lifestyle, and socio-economic factors.

Genetic Variants Underlying Plasticity in Natural Populations of Spadefoot Toads: Environmental Assessment versus Phenotypic Response

Many organisms facultatively produce different phenotypes depending on their environment, yet relatively little is known about the genetic bases of such plasticity in natural populations. In this study, we describe the genetic variation underlying an extreme form of plasticity--resource polyphenism--in Mexican spadefoot toad tadpoles, Spea multiplicata. Depending on their environment, these tadpoles develop into one of two drastically different forms: a carnivore morph or an omnivore morph. We collected both morphs from two ponds that differed in which morph had an adaptive advantage and performed genome-wide association studies of phenotype (carnivore vs. omnivore) and adaptive plasticity (adaptive vs. maladaptive environmental assessment). We identified four quantitative trait loci associated with phenotype and nine with adaptive plasticity, two of which exhibited signatures of minor allele dominance and two of which (one phenotype locus and one adaptive plasticity locus) did not occur as minor allele homozygotes. Investigations into the genetics of plastic traits in natural populations promise to provide novel insights into how such complex, adaptive traits arise and evolve.

Sex-specific circulating unconventional neutrophils determine immunological outcome of auto-inflammatory Behçet's uveitis

Neutrophils are the most abundant immune cells that first respond to insults in circulation. Although associative evidence suggests that differences in neutrophils may be linked to the sex-specific vulnerability of inflammatory diseases, mechanistic links remain elusive. Here, we identified extensive sex-specific heterogeneity in neutrophil composition under normal and auto-inflammatory conditions at single-cell resolution. Using a combination of single-cell RNA sequencing analysis, neutrophil-specific genetic knockouts and transfer experiments, we discovered dysregulation of two unconventional (interferon-α responsive and T cell regulatory) neutrophil subsets leading to male-biased incidence, severity and poor prognosis of auto-inflammatory Behçet's uveitis. Genome-wide association study (GWAS) and exosome study revealed that male-specific negative effects of both genetic factors and circulating exosomes on unconventional neutrophil subsets contributed to male-specific vulnerability to disease. Collectively, our findings identify sex-specifically distinct neutrophil subsets and highlight unconventional neutrophil subsets as sex-specific therapeutic targets to limit inflammatory diseases.