Unusual sequence characteristics of human chromosome 19 are conserved across 11 nonhuman primates

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.

Range-limited Heaps' law for functional DNA words in the human genome

Heaps' or Herdan-Heaps' law is a linguistic law describing the relationship between the vocabulary/dictionary size (type) and word counts (token) to be a power-law function. Its existence in genomes with certain definition of DNA words is unclear partly because the dictionary size in genome could be much smaller than that in a human language. We define a DNA word as a coding region in a genome that codes for a protein domain. Using human chromosomes and chromosome arms as individual samples, we establish the existence of Heaps' law in the human genome within limited range. Our definition of words in a genomic or proteomic context is different from other definitions such as over-represented k-mers which are much shorter in length. Although an approximate power-law distribution of protein domain sizes due to gene duplication and the related Zipf's law is well known, their translation to the Heaps' law in DNA words is not automatic. Several other animal genomes are shown herein also to exhibit range-limited Heaps' law with our definition of DNA words, though with various exponents. When tokens were randomly sampled and sample sizes reach to the maximum level, a deviation from the Heaps' law was observed, but a quadratic regression in log-log type-token plot fits the data perfectly. Investigation of type-token plot and its regression coefficients could provide an alternative narrative of reusage and redundancy of protein domains as well as creation of new protein domains from a linguistic perspective.

ENT3C: an entropy-based similarity measure for Hi-C and micro-C derived contact matrices

Hi-C and micro-C sequencing have shed light on the profound importance of 3D genome organization in cellular function by probing 3D contact frequencies across the linear genome. The resulting contact matrices are extremely sparse and susceptible to technical- and sequence-based biases, making their comparison challenging. The development of reliable, robust and efficient methods for quantifying similarity between contact matrices is crucial for investigating variations in the 3D genome organization in different cell types or under different conditions, as well as evaluating experimental reproducibility. We present a novel method, ENT3C, which measures the change in pattern complexity in the vicinity of contact matrix diagonals to quantify their similarity. ENT3C provides a robust, user-friendly Hi-C or micro-C contact matrix similarity metric and a characteristic entropy signal that can be used to gain detailed biological insights into 3D genome organization.

Cell-free DNA from germline TP53 mutation carriers reflect cancer-like fragmentation patterns

Germline pathogenic TP53 variants predispose individuals to a high lifetime risk of developing multiple cancers and are the hallmark feature of Li-Fraumeni syndrome (LFS). Our group has previously shown that LFS patients harbor shorter plasma cell-free DNA fragmentation; independent of cancer status. To understand the functional underpinning of cfDNA fragmentation in LFS, we conducted a fragmentomic analysis of 199 cfDNA samples from 82 TP53 mutation carriers and 30 healthy TP53-wildtype controls. We find that LFS individuals exhibit an increased prevalence of A/T nucleotides at fragment ends, dysregulated nucleosome positioning at p53 binding sites, and loci-specific changes in chromatin accessibility at development-associated transcription factor binding sites and at cancer-associated open chromatin regions. Machine learning classification resulted in robust differentiation between TP53 mutant versus wildtype cfDNA samples (AUC-ROC = 0.710-1.000) and intra-patient longitudinal analysis of ctDNA fragmentation signal enabled early cancer detection. These results suggest that cfDNA fragmentation may be a useful diagnostic tool in LFS patients and provides an important baseline for cancer early detection.

Ancestry-related differences in chromatin accessibility and gene expression of APOE ε4 are associated with Alzheimer's disease risk

INTRODUCTION

European local ancestry (ELA) surrounding apolipoprotein E (APOE) ε4 confers higher risk for Alzheimer's disease (AD) compared to African local ancestry (ALA). We demonstrated significantly higher APOE ε4  expression in ELA versus ALA in AD brains from APOE ε4/ε4 carriers. Chromatin accessibility differences could contribute to these expression changes.

METHODS

We performed single nuclei assays for transposase accessible chromatin sequencing from the frontal cortex of six ALA and six ELA AD brains, homozygous for local ancestry and APOE ε4.

RESULTS

Our results showed an increased chromatin accessibility at the APOE ε4  promoter area in ELA versus ALA astrocytes. This increased accessibility in ELA astrocytes extended genome wide. Genes with increased accessibility in ELA in astrocytes were enriched for synapsis, cholesterol processing, and astrocyte reactivity.

DISCUSSION

Our results suggest that increased chromatin accessibility of APOE ε4  in ELA astrocytes contributes to the observed elevated APOE ε4  expression, corresponding to the increased AD risk in ELA versus ALA APOE ε4/ε4 carriers.

The human placenta exhibits a unique transcriptomic void

The human placenta exhibits a unique genomic architecture with an unexpectedly high mutation burden and many uniquely expressed genes. The aim of this study is to identify transcripts that are uniquely absent or depleted in the placenta. Here, we show that 40 of 46 of the other organs have no selectively depleted transcripts and that, of the remaining six, the liver has the largest number, with 26. In contrast, the term placenta has 762 depleted transcripts. Gene Ontology analysis of this depleted set highlighted multiple pathways reflecting known unique elements of placental physiology. For example, transcripts associated with neuronal function are in the depleted set-as expected given the lack of placental innervation. However, this demonstrated overrepresentation of genes involved in mitochondrial function (p = 5.8 × 10-10), including PGC-1α, the master regulator of mitochondrial biogenesis, and genes involved in polyamine metabolism (p = 2.1 × 10-4).

Age-related methylation changes in the human sperm epigenome

Advanced paternal age is associated with increased risks for reproductive and offspring medical problems. Accumulating evidence suggests age-related changes in the sperm epigenome as one underlying mechanism. Using reduced representation bisulfite sequencing on 73 sperm samples of males attending a fertility center, we identified 1,162 (74%) regions which were significantly (FDR-adjusted) hypomethylated and 403 regions (26%) being hypermethylated with age. There were no significant correlations with paternal BMI, semen quality, or ART outcome. The majority (1,152 of 1,565; 74%) of age-related differentially methylated regions (ageDMRs) were located within genic regions, including 1,002 genes with symbols. Hypomethylated ageDMRs were closer to transcription start sites than hypermethylated DMRs, half of which reside in gene-distal regions. In this and conceptually related genome-wide studies, so far 2,355 genes have been reported with significant sperm ageDMRs, however most (90%) of them in only one study. The 241 genes which have been replicated at least once showed significant functional enrichments in 41 biological processes associated with development and the nervous system and in 10 cellular components associated with synapses and neurons. This supports the hypothesis that paternal age effects on the sperm methylome affect offspring behaviour and neurodevelopment. It is interesting to note that sperm ageDMRs were not randomly distributed throughout the human genome; chromosome 19 showed a highly significant twofold enrichment with sperm ageDMRs. Although the high gene density and CpG content have been conserved, the orthologous marmoset chromosome 22 did not appear to exhibit an increased regulatory potential by age-related DNA methylation changes.

SVhound: detection of regions that harbor yet undetected structural variation

BACKGROUND

Recent population studies are ever growing in number of samples to investigate the diversity of a population or species. These studies reveal new polymorphism that lead to important insights into the mechanisms of evolution, but are also important for the interpretation of these variations. Nevertheless, while the full catalog of variations across entire species remains unknown, we can predict which regions harbor additional not yet detected variations and investigate their properties, thereby enhancing the analysis for potentially missed variants.

RESULTS

To achieve this we developed SVhound ( https://github.com/lfpaulin/SVhound ), which based on a population level SVs dataset can predict regions that harbor unseen SV alleles. We tested SVhound using subsets of the 1000 genomes project data and showed that its correlation (average correlation of 2800 tests r = 0.7136) is high to the full data set. Next, we utilized SVhound to investigate potentially missed or understudied regions across 1KGP and CCDG. Lastly we also apply SVhound on a small and novel SV call set for rhesus macaque (Macaca mulatta) and discuss the impact and choice of parameters for SVhound.

CONCLUSIONS

SVhound is a unique method to identify potential regions that harbor hidden diversity in model and non model organisms and can also be potentially used to ensure high quality of SV call sets.

epialleleR: an R/Bioconductor package for sensitive allele-specific methylation analysis in NGS data

Low-level mosaic epimutations within the BRCA1 gene promoter occur in 5-8% of healthy individuals and are associated with a significantly elevated risk of breast and ovarian cancer. Similar events may also affect other tumor suppressor genes, potentially being a significant contributor to cancer burden. While this opens a new area for translational research, detection of low-level mosaic epigenetic events requires highly sensitive and robust methodology for methylation analysis. We here present epialleleR, a computational framework for sensitive detection, quantification, and visualization of mosaic epimutations in methylation sequencing data. Analyzing simulated and real data sets, we provide in-depth assessments of epialleleR performance and show that linkage to epihaplotype data is necessary to detect low-level methylation events. The epialleleR is freely available at https://github.com/BBCG/epialleleR and https://bioconductor.org/packages/epialleleR/ as an open-source R/Bioconductor package.

Maternal Copy Number Imbalances in Non-Invasive Prenatal Testing: Do They Matter?

Non-invasive prenatal testing (NIPT) has become a routine practice in screening for common aneuploidies of chromosomes 21, 18, and 13 and gonosomes X and Y in fetuses worldwide since 2015 and has even expanded to include smaller subchromosomal events. In fact, the fetal fraction represents only a small proportion of cell-free DNA on a predominant background of maternal DNA. Unlike fetal findings that have to be confirmed using invasive testing, it has been well documented that NIPT provides information on maternal mosaicism, occult malignancies, and hidden health conditions due to copy number variations (CNVs) with diagnostic resolution. Although large duplications or deletions associated with certain medical conditions or syndromes are usually well recognized and easy to interpret, very little is known about small, relatively common copy number variations on the order of a few hundred kilobases and their potential impact on human health. We analyzed data from 6422 NIPT patient samples with a CNV detection resolution of 200 kb for the maternal genome and identified 942 distinct CNVs; 328 occurred repeatedly. We defined them as multiple occurring variants (MOVs). We scrutinized the most common ones, compared them with frequencies in the gnomAD SVs v2.1, dbVar, and DGV population databases, and analyzed them with an emphasis on genomic content and potential association with specific phenotypes.

Active RB causes visible changes in nuclear organization

RB restricts G1/S progression by inhibiting E2F. Here, we show that sustained expression of active RB, and prolonged G1 arrest, causes visible changes in chromosome architecture that are not directly associated with E2F inhibition. Using FISH probes against two euchromatin RB-associated regions, two heterochromatin domains that lack RB-bound loci, and two whole-chromosome probes, we found that constitutively active RB (ΔCDK-RB) promoted a more diffuse, dispersed, and scattered chromatin organization. These changes were RB dependent, were driven by specific isoforms of monophosphorylated RB, and required known RB-associated activities. ΔCDK-RB altered physical interactions between RB-bound genomic loci, but the RB-induced changes in chromosome architecture were unaffected by dominant-negative DP1. The RB-induced changes appeared to be widespread and influenced chromosome localization within nuclei. Gene expression profiles revealed that the dispersion phenotype was associated with an increased autophagy response. We infer that, after cell cycle arrest, RB acts through noncanonical mechanisms to significantly change nuclear organization, and this reorganization correlates with transitions in cellular state.

Crosstalk between Hepatitis B Virus and the 3D Genome Structure

Viruses that transcribe their DNA within the nucleus have to adapt to the existing cellular mechanisms that govern transcriptional regulation. Recent technological breakthroughs have highlighted the highly hierarchical organization of the cellular genome and its role in the regulation of gene expression. This review provides an updated overview on the current knowledge on how the hepatitis B virus interacts with the cellular 3D genome and its consequences on viral and cellular gene expression. We also briefly discuss the strategies developed by other DNA viruses to co-opt and sometimes subvert cellular genome spatial organization.