Genome-wide DNA Methylation Profiling in Lyme Neuroborreliosis Reveals Altered Methylation Patterns of HLA Genes.

Lyme neuroborreliosis (LNB) is a complex neuroinflammatory disorder caused by Borrelia burgdorferi, which is transmitted through tick bites. Epigenetic alterations, specifically DNA methylation (DNAm), could play a role in the host immune response during infection. In this study, we present the first genome-wide analysis of DNAm in peripheral blood mononuclear cells from patients with LNB and those without LNB. Using a network-based approach, we highlighted HLA genes at the core of these DNAm changes, which were found to be enriched in immune-related pathways. These findings shed light on the role of epigenetic modifications in the LNB pathogenesis that should be confirmed and further expanded upon in future studies.

Multimodal classification of molecular subtypes in pediatric acute lymphoblastic leukemia.

Genomic analyses have redefined the molecular subgrouping of pediatric acute lymphoblastic leukemia (ALL). Molecular subgroups guide risk-stratification and targeted therapies, but outcomes of recently identified subtypes are often unclear, owing to limited cases with comprehensive profiling and cross-protocol studies. We developed a machine learning tool (ALLIUM) for the molecular subclassification of ALL in retrospective cohorts as well as for up-front diagnostics. ALLIUM uses DNA methylation and gene expression data from 1131 Nordic ALL patients to predict 17 ALL subtypes with high accuracy. ALLIUM was used to revise and verify the molecular subtype of 281 B-cell precursor ALL (BCP-ALL) cases with previously undefined molecular phenotype, resulting in a single revised subtype for 81.5% of these cases. Our study shows the power of combining DNA methylation and gene expression data for resolving ALL subtypes and provides a comprehensive population-based retrospective cohort study of molecular subtype frequencies in the Nordic countries.

Feasibility to use whole-genome sequencing as a sole diagnostic method to detect genomic aberrations in pediatric B-cell acute lymphoblastic leukemia.

Introduction: The suitability of whole-genome sequencing (WGS) as the sole method to detect clinically relevant genomic aberrations in B-cell acute lymphoblastic leukemia (ALL) was investigated with the aim of replacing current diagnostic methods. Methods: For this purpose, we assessed the analytical performance of 150 bp paired-end WGS (90x leukemia/30x germline). A set of 88 retrospective B-cell ALL samples were selected to represent established ALL subgroups as well as ALL lacking stratifying markers by standard-of-care (SoC), so-called B-other ALL. Results: Both the analysis of paired leukemia/germline (L/N)(n=64) as well as leukemia-only (L-only)(n=88) detected all types of aberrations mandatory in the current ALLTogether trial protocol, i.e., aneuploidies, structural variants, and focal copy-number aberrations. Moreover, comparison to SoC revealed 100% concordance and that all patients had been assigned to the correct genetic subgroup using both approaches. Notably, WGS could allocate 35 out of 39 B-other ALL samples to one of the emerging genetic subgroups considered in the most recent classifications of ALL. We further investigated the impact of high (90x; n=58) vs low (30x; n=30) coverage on the diagnostic yield and observed an equally perfect concordance with SoC; low coverage detected all relevant lesions. Discussion: The filtration of the WGS findings with a short list of genes recurrently rearranged in ALL was instrumental to extract the clinically relevant information efficiently. Nonetheless, the detection of DUX4 rearrangements required an additional customized analysis, due to multiple copies of this gene embedded in the highly repetitive D4Z4 region. We conclude that the diagnostic performance of WGS as the standalone method was remarkable and allowed detection of all clinically relevant genomic events in the diagnostic setting of B-cell ALL.

A frame-shift mutation in COMTD1 is associated with impaired pheomelanin pigmentation in chicken.

The biochemical pathway regulating the synthesis of yellow/red pheomelanin is less well characterized than the synthesis of black/brown eumelanin. Inhibitor of gold (IG phenotype) is a plumage colour variant in chicken that provides an opportunity to further explore this pathway since the recessive allele (IG) at this locus is associated with a defect in the production of pheomelanin. IG/IG homozygotes display a marked dilution of red pheomelanin pigmentation, whilst black pigmentation (eumelanin) is only slightly affected. Here we show that a 2-base pair insertion (frame-shift mutation) in the 5th exon of the Catechol-O-methyltransferase containing domain 1 gene (COMTD1), expected to cause a complete or partial loss-of-function of the COMTD1 enzyme, shows complete concordance with the IG phenotype within and across breeds. We show that the COMTD1 protein is localized to mitochondria in pigment cells. Knockout of Comtd1 in a mouse melanocytic cell line results in a reduction in pheomelanin metabolites and significant alterations in metabolites of glutamate/glutathione, riboflavin, and the tricarboxylic acid cycle. Furthermore, COMTD1 overexpression enhanced cellular proliferation following chemical-induced transfection, a potential inducer of oxidative stress. These observations suggest that COMTD1 plays a protective role for melanocytes against oxidative stress and that this supports their ability to produce pheomelanin.

The spectrum of tuberculosis described as differential DNA methylation patterns in alveolar macrophages and alveolar T cells.

BACKGROUND: Host innate immune cells have been identified as key players in the early eradication of Mycobacterium tuberculosis and in the maintenance of an anti-mycobacterial immune memory, which we and others have shown are induced through epigenetic reprogramming. Studies on human tuberculosis immunity are dominated by those using peripheral blood as surrogate markers for immunity. We aimed to investigate DNA methylation patterns in immune cells of the lung compartment by obtaining induced sputum from M. tuberculosis- exposed subjects including symptom-free subjects testing positively and negatively for latent tuberculosis as well as patients diagnosed with active tuberculosis. Alveolar macrophages and alveolar T cells were isolated from the collected sputum and DNA methylome analyses performed (Illumina Infinium Human Methylation 450 k). RESULTS: Multidimensional scaling analysis revealed that DNA methylomes of cells from the tuberculosis-exposed subjects and controls appeared as separate clusters. The numerous genes that were differentially methylated between the groups were functionally connected and overlapped with previous findings of trained immunity and tuberculosis. In addition, analysis of the interferon-gamma release assay (IGRA) status of the subjects demonstrated that the IGRA status was reflected in the DNA methylome by a unique signature. CONCLUSIONS: This pilot study suggests that M. tuberculosis induces epigenetic reprogramming in immune cells of the lung compartment, reflected as a specific DNA methylation pattern. The DNA methylation signature emerging from the comparison of IGRA-negative and IGRA-positive subjects revealed a spectrum of signature strength with the TB patients grouping together at one end of the spectrum, both in alveolar macrophages and T cells. DNA methylation-based biosignatures could be considered for further development towards a clinically useful tool for determining tuberculosis infection status and the level of tuberculosis exposure.

Defining post-acute COVID-19 syndrome (PACS) by an epigenetic biosignature in peripheral blood mononuclear cells.

Post-acute COVID-19 syndrome (PACS) has been defined as symptoms persisting after clearance of a COVID-19 infection. We have previously demonstrated that alterations in DNA methylation (DNAm) status persist in individuals who recovered from a COVID-19 infection, but it is currently unknown if PACS is associated with epigenetic changes. We compared DNAm patterns in patients with PACS with those in controls and in healthy COVID-19 convalescents and found a unique DNAm signature in PACS patients. This signature unravelled modified pathways that regulate angiotensin II and muscarinic receptor signalling and protein-protein interaction networks that have bearings on vesicle formation and mitochondrial function.

Epigenetic rewiring of pathways related to odour perception in immune cells exposed to SARS-CoV-2 in vivo and in vitro.

A majority of SARS-CoV-2 recoverees develop only mild-to-moderate symptoms, while some remain completely asymptomatic. Although viruses, including SARS-CoV-2, may evade host immune responses by epigenetic mechanisms including DNA methylation, little is known about whether these modifications are important in defence against and healthy recovery from COVID-19 in the host. To this end, epigenome-wide DNA methylation patterns from COVID-19 convalescents were compared to uninfected controls from before and after the pandemic. Peripheral blood mononuclear cell (PBMC) DNA was extracted from uninfected controls, COVID-19 convalescents, and symptom-free individuals with SARS-CoV-2-specific T cell-responses, as well as from PBMCs stimulated in vitro with SARS-CoV-2. Subsequently, the Illumina MethylationEPIC 850K array was performed, and statistical/bioinformatic analyses comprised differential DNA methylation, pathway over-representation, and module identification analyses. Differential DNA methylation patterns distinguished COVID-19 convalescents from uninfected controls, with similar results in an experimental SARS-CoV-2 infection model. A SARS-CoV-2-induced module was identified in vivo, comprising 66 genes of which six (TP53, INS, HSPA4, SP1, ESR1, and FAS) were present in corresponding in vitro analyses. Over-representation analyses revealed involvement in Wnt, muscarinic acetylcholine receptor signalling, and gonadotropin-releasing hormone receptor pathways. Furthermore, numerous differentially methylated and network genes from both settings interacted with the SARS-CoV-2 interactome. Altered DNA methylation patterns of COVID-19 convalescents suggest recovery from mild-to-moderate SARS-CoV-2 infection leaves longstanding epigenetic traces. Both in vitro and in vivo exposure caused epigenetic modulation of pathways thataffect odour perception. Future studies should determine whether this reflects host-induced protective antiviral defense or targeted viral hijacking to evade host defence.

Hop Mice Display Synchronous Hindlimb Locomotion and a Ventrally Fused Lumbar Spinal Cord Caused by a Point Mutation in Ttc26.

Identifying the spinal circuits controlling locomotion is critical for unravelling the mechanisms controlling the production of gaits. Development of the circuits governing left-right coordination relies on axon guidance molecules such as ephrins and netrins. To date, no other class of proteins have been shown to play a role during this process. Here, we have analyzed hop mice, which walk with a characteristic hopping gait using their hindlimbs in synchrony. Fictive locomotion experiments suggest that a local defect in the ventral spinal cord contributes to the aberrant locomotor phenotype. Hop mutant spinal cords had severe morphologic defects, including the absence of the ventral midline and a poorly defined border between white and gray matter. The hop mice represent the first model where, exclusively found in the lumbar domain, the left and right components of the central pattern generators (CPGs) are fused with a synchronous hindlimb gait as a functional consequence. These defects were associated with abnormal developmental processes, including a misplaced notochord and reduced induction of ventral progenitor domains. Whereas the underlying mutation in hop mice has been suggested to lie within the Ttc26 gene, other genes in close vicinity have been associated with gait defects. Mouse embryos carrying a CRISPR replicated point mutation within Ttc26 displayed an identical morphologic phenotype. Thus, our data suggest that the assembly of the lumbar CPG network is dependent on fully functional TTC26 protein.

A differential DNA methylome signature of pulmonary immune cells from individuals converting to latent tuberculosis infection.

Tuberculosis (TB), caused by Mycobacterium tuberculosis, spreads via aerosols and the first encounter with the immune system is with the pulmonary-resident immune cells. The role of epigenetic regulations in the immune cells is emerging and we have previously shown that macrophages capacity to kill M. tuberculosis is reflected in the DNA methylome. The aim of this study was to investigate epigenetic modifications in alveolar macrophages and T cells in a cohort of medical students with an increased risk of TB exposure, longitudinally. DNA methylome analysis revealed that a unique DNA methylation profile was present in healthy subjects who later developed latent TB during the study. The profile was reflected in a different overall DNA methylation distribution as well as a distinct set of differentially methylated genes (DMGs). The DMGs were over-represented in pathways related to metabolic reprogramming of macrophages and T cell migration and IFN-γ production, pathways previously reported important in TB control. In conclusion, we identified a unique DNA methylation signature in individuals, with no peripheral immune response to M. tuberculosis antigen who later developed latent TB. Together the study suggests that the DNA methylation status of pulmonary immune cells can reveal who will develop latent TB infection.

Mutational patterns and clonal evolution from diagnosis to relapse in pediatric acute lymphoblastic leukemia.

The mechanisms driving clonal heterogeneity and evolution in relapsed pediatric acute lymphoblastic leukemia (ALL) are not fully understood. We performed whole genome sequencing of samples collected at diagnosis, relapse(s) and remission from 29 Nordic patients. Somatic point mutations and large-scale structural variants were called using individually matched remission samples as controls, and allelic expression of the mutations was assessed in ALL cells using RNA-sequencing. We observed an increased burden of somatic mutations at relapse, compared to diagnosis, and at second relapse compared to first relapse. In addition to 29 known ALL driver genes, of which nine genes carried recurrent protein-coding mutations in our sample set, we identified putative non-protein coding mutations in regulatory regions of seven additional genes that have not previously been described in ALL. Cluster analysis of hundreds of somatic mutations per sample revealed three distinct evolutionary trajectories during ALL progression from diagnosis to relapse. The evolutionary trajectories provide insight into the mutational mechanisms leading relapse in ALL and could offer biomarkers for improved risk prediction in individual patients.

A genomic map of clinal variation across the European rabbit hybrid zone.

Speciation is a process proceeding from weak to complete reproductive isolation. In this continuum, naturally hybridizing taxa provide a promising avenue for revealing the genetic changes associated with the incipient stages of speciation. To identify such changes between two subspecies of rabbits that display partial reproductive isolation, we studied patterns of allele frequency change across their hybrid zone using whole-genome sequencing. To connect levels and patterns of genetic differentiation with phenotypic manifestations of subfertility in hybrid rabbits, we further investigated patterns of gene expression in testis. Geographic cline analysis revealed 253 regions characterized by steep changes in allele frequency across their natural region of contact. This catalog of regions is likely to be enriched for loci implicated in reproductive barriers and yielded several insights into the evolution of hybrid dysfunction in rabbits: (i) incomplete reproductive isolation is likely governed by the effects of many loci, (ii) protein-protein interaction analysis suggest that genes within these loci interact more than expected by chance, (iii) regulatory variation is likely the primary driver of incompatibilities, and (iv) large chromosomal rearrangements appear not to be a major mechanism underlying incompatibilities or promoting isolation in the face of gene flow. We detected extensive misregulation of gene expression in testis of hybrid males, but not a statistical overrepresentation of differentially expressed genes in candidate regions. Our results also did not support an X chromosome-wide disruption of expression as observed in mice and cats, suggesting variation in the mechanistic basis of hybrid male reduced fertility among mammals.

Whole-Genome Sequencing of a Canine Family Trio Reveals a FAM83G Variant Associated with Hereditary Footpad Hyperkeratosis.

Over 250 Mendelian traits and disorders, caused by rare alleles have been mapped in the canine genome. Although each disease is rare in the dog as a species, they are collectively common and have major impact on canine health. With SNP-based genotyping arrays, genome-wide association studies (GWAS) have proven to be a powerful method to map the genomic region of interest when 10-20 cases and 10-20 controls are available. However, to identify the genetic variant in associated regions, fine-mapping and targeted resequencing is required. Here we present a new approach using whole-genome sequencing (WGS) of a family trio without prior GWAS. As a proof-of-concept, we chose an autosomal recessive disease known as hereditary footpad hyperkeratosis (HFH) in Kromfohrländer dogs. To our knowledge, this is the first time this family trio WGS-approach has been used successfully to identify a genetic variant that perfectly segregates with a canine disorder. The sequencing of three Kromfohrländer dogs from a family trio (an affected offspring and both its healthy parents) resulted in an average genome coverage of 9.2X per individual. After applying stringent filtering criteria for candidate causative coding variants, 527 single nucleotide variants (SNVs) and 15 indels were found to be homozygous in the affected offspring and heterozygous in the parents. Using the computer software packages ANNOVAR and SIFT to functionally annotate coding sequence differences, and to predict their functional effect, resulted in seven candidate variants located in six different genes. Of these, only FAM83G:c155G > C (p.R52P) was found to be concordant in eight additional cases, and 16 healthy Kromfohrländer dogs.

Evaluation of whole-genome sequencing of four Chinese crested dogs for variant detection using the ion proton system.

BACKGROUND: Next generation sequencing (NGS) has traditionally been performed by large genome centers, but in recent years, the costs for whole-genome sequencing (WGS) have decreased substantially. With the introduction of smaller and less expensive "desktop" systems, NGS is now moving into the general laboratory. To evaluate the Ion Proton system for WGS we sequenced four Chinese Crested dogs and analyzed the data quality in terms of genome and exome coverage, the number of detected single nucleotide variants (SNVs) and insertions and deletions (INDELs) and the genotype concordance with the Illumina HD canine SNP array. For each of the four dogs, a 200 bp fragment library was constructed from genomic DNA and sequenced on two Ion PI chips per dog to reach mean coverage of 6-8x of the canine genome (genome size ≈ 2.4 Gb). RESULTS: On average, each Ion PI chip yielded approximately 73.3 million reads with a mean read length of 130 bp (~9.5 Gb sequence data) of which 98.5 % could be aligned to the canine reference genome (CanFam3.1). By sequencing a single dog using one fragment library and two Ion PI chips, on average 80 % of the genome and 77 % exome was covered by at least four reads. After removing duplicate reads (20.7 %) the mean coverage across the whole genome was 6x. Using sequence data from all four individuals (four fragment libraries and eight Ion PI chips) the genome and exome coverage could be further increased to 97.2 and 94.3 %, respectively. We detected 4.83 million unique SNPs and 6.10 million unique INDEL positions across all individuals. A comparison between SNP genotypes detected with the WGS and the 170 K Illumina HD canine SNP array showed 90 % concordance. CONCLUSIONS: We have evaluated whole-genome sequencing on the Ion Proton system for genetic variant detection in four Chinese crested dogs. Even though INDEL calling with Ion Proton data is challenging due to specific platform errors, in case of SNP calling it can serve as an alternative to other next-generation sequencing platforms and SNP genotyping arrays, in studies aiming to identify causative mutations for rare monogenic diseases. In addition, we have identified new genetic variants of the Chinese Crested dog that will contribute to further whole-genome sequencing studies aimed to identify mutations associated with monogenic diseases with autosomal recessive inheritance.

Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication.

The genetic changes underlying the initial steps of animal domestication are still poorly understood. We generated a high-quality reference genome for the rabbit and compared it to resequencing data from populations of wild and domestic rabbits. We identified more than 100 selective sweeps specific to domestic rabbits but only a relatively small number of fixed (or nearly fixed) single-nucleotide polymorphisms (SNPs) for derived alleles. SNPs with marked allele frequency differences between wild and domestic rabbits were enriched for conserved noncoding sites. Enrichment analyses suggest that genes affecting brain and neuronal development have often been targeted during domestication. We propose that because of a truly complex genetic background, tame behavior in rabbits and other domestic animals evolved by shifts in allele frequencies at many loci, rather than by critical changes at only a few domestication loci.

Strong signatures of selection in the domestic pig genome.

Domestication of wild boar (Sus scrofa) and subsequent selection have resulted in dramatic phenotypic changes in domestic pigs for a number of traits, including behavior, body composition, reproduction, and coat color. Here we have used whole-genome resequencing to reveal some of the loci that underlie phenotypic evolution in European domestic pigs. Selective sweep analyses revealed strong signatures of selection at three loci harboring quantitative trait loci that explain a considerable part of one of the most characteristic morphological changes in the domestic pig--the elongation of the back and an increased number of vertebrae. The three loci were associated with the NR6A1, PLAG1, and LCORL genes. The latter two have repeatedly been associated with loci controlling stature in other domestic animals and in humans. Most European domestic pigs are homozygous for the same haplotype at these three loci. We found an excess of derived nonsynonymous substitutions in domestic pigs, most likely reflecting both positive selection and relaxed purifying selection after domestication. Our analysis of structural variation revealed four duplications at the KIT locus that were exclusively present in white or white-spotted pigs, carrying the Dominant white, Patch, or Belt alleles. This discovery illustrates how structural changes have contributed to rapid phenotypic evolution in domestic animals and how alleles in domestic animals may evolve by the accumulation of multiple causative mutations as a response to strong directional selection.