Combined approaches, including long-read sequencing, address the diagnostic challenge of HYDIN in primary ciliary dyskinesia.

Primary ciliary dyskinesia (PCD), a disorder of the motile cilia, is now recognised as an underdiagnosed cause of bronchiectasis. Accurate PCD diagnosis comprises clinical assessment, analysis of cilia and the identification of biallelic variants in one of 50 known PCD-related genes, including HYDIN. HYDIN-related PCD is underdiagnosed due to the presence of a pseudogene, HYDIN2, with 98% sequence homology to HYDIN. This presents a significant challenge for Short-Read Next Generation Sequencing (SR-NGS) and analysis, and many diagnostic PCD gene panels do not include HYDIN. We have used a combined approach of SR-NGS with bioinformatic masking of HYDIN2, and state-of-the-art long-read Nanopore sequencing (LR_NGS), together with analysis of respiratory cilia including transmission electron microscopy and immunofluorescence to address the underdiagnosis of HYDIN as a cause of PCD. Bioinformatic masking of HYDIN2 after SR-NGS facilitated the detection of biallelic HYDIN variants in 15 of 437 families, but compromised the detection of copy number variants. Supplementing testing with LR-NGS detected HYDIN deletions in 2 families, where SR-NGS had detected a single heterozygous HYDIN variant. LR-NGS was also able to confirm true homozygosity in 2 families when parental testing was not possible. Utilising a combined genomic diagnostic approach, biallelic HYDIN variants were detected in 17 families from 242 genetically confirmed PCD cases, comprising 7% of our PCD cohort. This represents the largest reported HYDIN cohort to date and highlights previous underdiagnosis of HYDIN-associated PCD. Moreover this provides further evidence for the utility of LR-NGS in diagnostic testing, particularly for regions of high genomic complexity.

A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development.

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.

CLN2 disease resulting from a novel homozygous deep intronic splice variant in TPP1 discovered using long-read sequencing.

Neuronal ceroid lipofuscinosis type 2 (CLN2) is an autosomal recessive neurodegenerative disorder with enzyme replacement therapy available. We present two siblings with a clinical diagnosis of CLN2 disease, but no identifiable TPP1 variants after standard clinical testing. Long-read sequencing identified a homozygous deep intronic variant predicted to affect splicing, confirmed by clinical DNA and RNA sequencing. This case demonstrates how traditional laboratory assays can complement emerging molecular technologies to provide a precise molecular diagnosis.

Diagnostic Utility of Genome-wide DNA Methylation Analysis in Genetically Unsolved Developmental and Epileptic Encephalopathies and Refinement of a CHD2 Episignature.

Sequence-based genetic testing currently identifies causative genetic variants in ∼50% of individuals with developmental and epileptic encephalopathies (DEEs). Aberrant changes in DNA methylation are implicated in various neurodevelopmental disorders but remain unstudied in DEEs. Rare epigenetic variations ("epivariants") can drive disease by modulating gene expression at single loci, whereas genome-wide DNA methylation changes can result in distinct "episignature" biomarkers for monogenic disorders in a growing number of rare diseases. Here, we interrogate the diagnostic utility of genome-wide DNA methylation array analysis on peripheral blood samples from 516 individuals with genetically unsolved DEEs who had previously undergone extensive genetic testing. We identified rare differentially methylated regions (DMRs) and explanatory episignatures to discover causative and candidate genetic etiologies in 10 individuals. We then used long-read sequencing to identify DNA variants underlying rare DMRs, including one balanced translocation, three CG-rich repeat expansions, and two copy number variants. We also identify pathogenic sequence variants associated with episignatures; some had been missed by previous exome sequencing. Although most DEE genes lack known episignatures, the increase in diagnostic yield for DNA methylation analysis in DEEs is comparable to the added yield of genome sequencing. Finally, we refine an episignature for CHD2 using an 850K methylation array which was further refined at higher CpG resolution using bisulfite sequencing to investigate potential insights into CHD2 pathophysiology. Our study demonstrates the diagnostic yield of genome-wide DNA methylation analysis to identify causal and candidate genetic causes as ∼2% (10/516) for unsolved DEE cases.

Genome Report: chromosome-scale genome assembly of the African spiny mouse (Acomys cahirinus).

There is increasing interest in the African spiny mouse (Acomys cahirinus) as a model organism because of its ability for regeneration of tissue after injury in skin, muscle, and internal organs such as the kidneys. A high-quality reference genome is needed to better understand these regenerative properties at the molecular level. Here, we present an improved reference genome for A. cahirinus generated from long Nanopore sequencing reads. We confirm the quality of our annotations using RNA sequencing data from 4 different tissues. Our genome is of higher contiguity and quality than previously reported genomes from this species and will facilitate ongoing efforts to better understand the regenerative properties of this organism.

GENOME REPORT: Chromosome-scale genome assembly of the African spiny mouse ( Acomys cahirinus ).

There is increasing interest in the African spiny mouse ( Acomys cahirinus ) as a model organism because of its ability for regeneration of tissue after injury in skin, muscle, and internal organs such as the kidneys. A high-quality reference genome is needed to better understand these regenerative properties at the molecular level. Here, we present an improved reference genome for A. cahirinus generated from long Nanopore sequencing reads. We confirm the quality of our annotations using RNA sequencing data from four different tissues. Our genome is of higher contiguity and quality than previously reported genomes from this species and will facilitate ongoing efforts to better understand the regenerative properties of this organism.

Heterozygous inversion breakpoints suppress meiotic crossovers by altering recombination repair outcomes.

Heterozygous chromosome inversions suppress meiotic crossover (CO) formation within an inversion, potentially because they lead to gross chromosome rearrangements that produce inviable gametes. Interestingly, COs are also severely reduced in regions nearby but outside of inversion breakpoints even though COs in these regions do not result in rearrangements. Our mechanistic understanding of why COs are suppressed outside of inversion breakpoints is limited by a lack of data on the frequency of noncrossover gene conversions (NCOGCs) in these regions. To address this critical gap, we mapped the location and frequency of rare CO and NCOGC events that occurred outside of the dl-49 chrX inversion in D. melanogaster. We created full-sibling wildtype and inversion stocks and recovered COs and NCOGCs in the syntenic regions of both stocks, allowing us to directly compare rates and distributions of recombination events. We show that COs outside of the proximal inversion breakpoint are distributed in a distance-dependent manner, with strongest suppression near the inversion breakpoint. We find that NCOGCs occur evenly throughout the chromosome and, importantly, are not suppressed near inversion breakpoints. We propose a model in which COs are suppressed by inversion breakpoints in a distance-dependent manner through mechanisms that influence DNA double-strand break repair outcome but not double-strand break formation. We suggest that subtle changes in the synaptonemal complex and chromosome pairing might lead to unstable interhomolog interactions during recombination that permits NCOGC formation but not CO formation.

Targeted Long-Read Sequencing Identifies a Retrotransposon Insertion as a Cause of Altered GNAS Exon A/B Methylation in a Family With Autosomal Dominant Pseudohypoparathyroidism Type 1b (PHP1B).

Pseudohypoparathyroidism type Ib (PHP1B) is characterized predominantly by resistance to parathyroid hormone (PTH) leading to hypocalcemia and hyperphosphatemia. These laboratory abnormalities are caused by maternal loss-of-methylation (LOM) at GNAS exon A/B, which reduces in cis expression of the stimulatory G protein α-subunit (Gsα). Paternal Gsα expression in proximal renal tubules is silenced through unknown mechanisms, hence LOM at exon A/B reduces further Gsα protein in this kidney portion, leading to PTH resistance. In a previously reported PHP1B family, affected members showed variable LOM at exon A/B, yet no genetic defect was found by whole-genome sequencing despite linkage to GNAS. Using targeted long-read sequencing (T-LRS), we discovered an approximately 2800-bp maternally inherited retrotransposon insertion nearly 1200 bp downstream of exon XL not found in public databases or in 13,675 DNA samples analyzed by short-read whole-genome sequencing. T-LRS data furthermore confirmed normal methylation at exons XL, AS, and NESP and showed that LOM comprising exon A/B is broader than previously thought. The retrotransposon most likely causes the observed epigenetic defect by impairing function of a maternally derived NESP transcript, consistent with findings in mice lacking full-length NESP mRNA and in PHP1B patients with deletion of exon NESP and adjacent intronic sequences. In addition to demonstrating that T-LRS is an effective strategy for identifying a small disease-causing variant that abolishes or severely reduces exon A/B methylation, our data demonstrate that this sequencing technology has major advantages for simultaneously identifying structural defects and altered methylation. © 2022 American Society for Bone and Mineral Research (ASBMR).

Targeted long-read sequencing identifies missing pathogenic variants in unsolved Werner syndrome cases.

BACKGROUND: Werner syndrome (WS) is an autosomal recessive progeroid syndrome caused by variants in WRN. The International Registry of Werner Syndrome has identified biallelic pathogenic variants in 179/188 cases of classical WS. In the remaining nine cases, only one heterozygous pathogenic variant has been identified. METHODS: Targeted long-read sequencing (T-LRS) on an Oxford Nanopore platform was used to search for a second pathogenic variant in WRN. Previously, T-LRS was successfully used to identify missing variants and analyse complex rearrangements. RESULTS: We identified a second pathogenic variant in eight of nine unsolved WS cases. In five cases, T-LRS identified intronic splice variants that were confirmed by either RT-PCR or exon trapping to affect splicing; in one case, T-LRS identified a 339 kbp deletion, and in two cases, pathogenic missense variants. Phasing of long reads predicted all newly identified variants were on a different haplotype than the previously known variant. Finally, in one case, RT-PCR previously identified skipping of exon 20; however, T-LRS did not detect a pathogenic DNA sequence variant. CONCLUSION: T-LRS is an effective method for identifying missing pathogenic variants. Although limitations with computational prediction algorithms can hinder the interpretation of variants, T-LRS is particularly effective in identifying intronic variants.

Enrichment of psychrophilic and acidophilic sulfate-reducing bacterial consortia - a solution toward acid mine drainage treatment in cold regions.

Failure of sulfate-reducing bacteria (SRB)-mediated treatment of acid mine drainage (AMD) in cold regions due to inhibition of bacteria by acidic pH and low temperature can be overcome by enriching psychrophilic and acidophilic microbial consortia from local metal-rich sediments. In this study, we enriched microbial consortia from Arctic mine sediments at varying pH (3-7) and temperatures (15-37 °C) under anaerobic conditions with repeated sub-culturing in three successive stages, and analyzed the microbial community using 16S rRNA gene sequencing. The enriched SRB genera resulted in high sulfate reduction (85-88%), and significant metal removal (49-99.9%) during the initial stages (stage 1 and 2). Subsequently, sub-culturing the inoculum at pH 3-4.5 resulted in lower sulfate reduction (9-34%) due to the inhibition of SRB by accumulated acetic acid (0.3-9 mM). The microbial metabolic interactions for successful sulfate and metal removal involved initial glycerol co-fermentation to acetic acid at acidic pH (by Desulfosporosinus, Desulfotomaculum, Desulfurospora, and fermentative bacteria including Cellulomonas and Anaerovorax), followed by acetic acid oxidation to CO2 and H2 (by Desulfitobacterium) at neutral pH, and subsequent H2 utilization (by Desulfosporosinus). The results, including the structural and functional properties of enriched microbial consortia, can inform the development of effective biological treatment strategies for AMD in cold regions.

Probabilistic models of biological enzymatic polymerization.

In this study, hierarchies of probabilistic models are evaluated for their ability to characterize the untemplated addition of adenine and uracil to the 3' ends of mitochondrial mRNAs of the human pathogen Trypanosoma brucei, and for their generative abilities to reproduce populations of these untemplated adenine/uridine "tails". We determined the most ideal Hidden Markov Models (HMMs) for this biological system. While our HMMs were not able to generatively reproduce the length distribution of the tails, they fared better in reproducing nucleotide composition aspects of the tail populations. The HMMs robustly identified distinct states of nucleotide addition that correlate to experimentally verified tail nucleotide composition differences. However they also identified a surprising subclass of tails among the ND1 gene transcript populations that is unexpected given the current idea of sequential enzymatic action of untemplated tail addition in this system. Therefore, these models can not only be utilized to reflect biological states that we already know about, they can also identify hypotheses to be experimentally tested. Finally, our HMMs supplied a way to correct a portion of the sequencing errors present in our data. Importantly, these models constitute rare simple pedagogical examples of applied bioinformatic HMMs, due to their binary emissions.

Ecosystem-bedrock interaction changes nutrient compartmentalization during early oxidative weathering.

Ecosystem-bedrock interactions power the biogeochemical cycles of Earth's shallow crust, supporting life, stimulating substrate transformation, and spurring evolutionary innovation. While oxidative processes have dominated half of terrestrial history, the relative contribution of the biosphere and its chemical fingerprints on Earth's developing regolith are still poorly constrained. Here, we report results from a two-year incipient weathering experiment. We found that the mass release and compartmentalization of major elements during weathering of granite, rhyolite, schist and basalt was rock-specific and regulated by ecosystem components. A tight interplay between physiological needs of different biota, mineral dissolution rates, and substrate nutrient availability resulted in intricate elemental distribution patterns. Biota accelerated CO2 mineralization over abiotic controls as ecosystem complexity increased, and significantly modified the stoichiometry of mobilized elements. Microbial and fungal components inhibited element leaching (23.4% and 7%), while plants increased leaching and biomass retention by 63.4%. All biota left comparable biosignatures in the dissolved weathering products. Nevertheless, the magnitude and allocation of weathered fractions under abiotic and biotic treatments provide quantitative evidence for the role of major biosphere components in the evolution of upper continental crust, presenting critical information for large-scale biogeochemical models and for the search for stable in situ biosignatures beyond Earth.