A corpus of GA4GH Phenopackets: case-level phenotyping for genomic diagnostics and discovery.

The Global Alliance for Genomics and Health (GA4GH) Phenopacket Schema was released in 2022 and approved by ISO as a standard for sharing clinical and genomic information about an individual, including phenotypic descriptions, numerical measurements, genetic information, diagnoses, and treatments. A phenopacket can be used as an input file for software that supports phenotype-driven genomic diagnostics and for algorithms that facilitate patient classification and stratification for identifying new diseases and treatments. There has been a great need for a collection of phenopackets to test software pipelines and algorithms. Here, we present phenopacket-store. Version 0.1.12 of phenopacket-store includes 4916 phenopackets representing 277 Mendelian and chromosomal diseases associated with 236 genes, and 2872 unique pathogenic alleles curated from 605 different publications. This represents the first large-scale collection of case-level, standardized phenotypic information derived from case reports in the literature with detailed descriptions of the clinical data and will be useful for many purposes, including the development and testing of software for prioritizing genes and diseases in diagnostic genomics, machine learning analysis of clinical phenotype data, patient stratification, and genotype-phenotype correlations. This corpus also provides best-practice examples for curating literature-derived data using the GA4GH Phenopacket Schema.

Array-based dynamic allele specific hybridization (Array-DASH): Optimization-free microarray processing for multiple simultaneous genomic assays.

We report proof-of-principle experiments regarding a dynamic microarray protocol enabling accurate and semi-quantitative DNA analysis for re-sequencing, fingerprinting and genotyping. Single-stranded target molecules hybridise to surface-bound probes during initial gradual cooling with high-fidelity. Real-time tracking of target denaturation (via fluorescence) during a 'dynamic' gradual heating phase permits 'melt-curve' analysis. The probe most closely matching the target sequence is identified based on the highest melting temperature. We demonstrated a >99% re-sequencing accuracy and a potential detection rate of 1% for SNPs. Experiments employing Hypericum ribosomal ITS regions and HIV genomes illustrated a reliable detection level of 5% plus simultaneous re-sequencing and genotyping. Such performance suggests a range of potential real-world applications involving rapid sequence interrogation, for example, in the Covid-19 pandemic. Guidance is offered towards the development of a commercial platform and dedicated software required to bring this technique into mainstream science.

Verifying nomenclature of DNA variants in submitted manuscripts: Guidance for journals.

Documenting variation in our genomes is important for research and clinical care. Accuracy in the description of DNA variants is therefore essential. To address this issue, the Human Variome Project convened a committee to evaluate the feasibility of requiring authors to verify that all variants submitted for publication complied with a widely accepted standard for description. After a pilot study of two journals, the committee agreed that requiring authors to verify that variants complied with Human Genome Variation Society nomenclature is a reasonable step toward standardizing the worldwide inventory of human variation.

The first reported case of the rare mitochondrial haplotype H4a1 in ancient Egypt.

Takabuti, was a female who lived in ancient Egypt during the 25th Dynasty, c.660 BCE. Her mummified remains were brought to Belfast, Northern Ireland, in 1834 and are currently displayed in the Ulster Museum. To gain insight into Takabuti's ancestry, we used deep sampling of vertebral bone, under X-ray control, to obtain non-contaminated bone tissue from which we extracted ancient DNA (aDNA) using established protocols. We targeted the maternally inherited mitochondrial DNA (mtDNA), known to be highly informative for human ancestry, and identified 38 single nucleotide variants using next generation sequencing. The specific combination of these SNVs suggests that Takabuti belonged to mitochondrial haplogroup H4a1. Neither H4 nor H4a1 have been reported in ancient Egyptian samples, prior to this study. The modern distribution of H4a1 is rare and sporadic and has been identified in areas including the Canary Islands, southern Iberia and the Lebanon. H4a1 has also been reported in ancient samples from Bell Beaker and Unetice contexts in Germany, as well as Bronze Age Bulgaria. We believe that this is an important finding because first, it adds to the depth of knowledge about the distribution of the H4a1 haplogroup in existing mtDNA, thus creating a baseline for future occurrences of this haplogroup in ancient Egyptian remains. Second, it is of great importance for archaeological sciences, since a predominantly European haplogroup has been identified in an Egyptian individual in Southern Egypt, prior to the Roman and Greek influx (332BCE).

hgvs: A Python package for manipulating sequence variants using HGVS nomenclature: 2018 Update.

The Human Genome Variation Society (HGVS) nomenclature guidelines encourage the accurate and standard description of DNA, RNA, and protein sequence variants in public variant databases and the scientific literature. Inconsistent application of the HGVS guidelines can lead to misinterpretation of variants in clinical settings. Reliable software tools are essential to ensure consistent application of the HGVS guidelines when reporting and interpreting variants. We present the hgvs Python package, a comprehensive tool for manipulating sequence variants according to the HGVS nomenclature guidelines. Distinguishing features of the hgvs package include: (1) parsing, formatting, validating, and normalizing variants on genome, transcript, and protein sequences; (2) projecting variants between aligned sequences, including those with gapped alignments; (3) flexible installation using remote or local data (fully local installations eliminate network dependencies); (4) extensive automated tests; and (5) open source development by a community from eight organizations worldwide. This report summarizes recent and significant updates to the hgvs package since its original release in 2014, and presents results of extensive validation using clinical relevant variants from ClinVar and HGMD.

SMYD2 promoter DNA methylation is associated with abdominal aortic aneurysm (AAA) and SMYD2 expression in vascular smooth muscle cells.

Background: Abdominal aortic aneurysm (AAA) is a deadly cardiovascular disease characterised by the gradual, irreversible dilation of the abdominal aorta. AAA is a complex genetic disease but little is known about the role of epigenetics. Our objective was to determine if global DNA methylation and CpG-specific methylation at known AAA risk loci is associated with AAA, and the functional effects of methylation changes. Results: We assessed global methylation in peripheral blood mononuclear cell DNA from 92 individuals with AAA and 93 controls using enzyme-linked immunosorbent assays, identifying hyper-methylation in those with large AAA and a positive linear association with AAA diameter (P < 0.0001, R2 = 0.3175).We then determined CpG methylation status of regulatory regions in genes located at AAA risk loci identified in genome-wide association studies, using bisulphite next-generation sequencing (NGS) in vascular smooth muscle cells (VSMCs) taken from aortic tissues of 44 individuals (24 AAAs and 20 controls). In IL6R, 2 CpGs were hyper-methylated (P = 0.0145); in ERG, 13 CpGs were hyper-methylated (P = 0.0005); in SERPINB9, 6 CpGs were hypo-methylated (P = 0.0037) and 1 CpG was hyper-methylated (P = 0.0098); and in SMYD2, 4 CpGs were hypo-methylated (P = 0.0012).RT-qPCR was performed for each differentially methylated gene on mRNA from the same VSMCs and compared with methylation. This analysis revealed downregulation of SMYD2 and SERPINB9 in AAA, and a direct linear relationship between SMYD2 promoter methylation and SMYD2 expression (P = 0.038). Furthermore, downregulation of SMYD2 at the site of aneurysm in the aortic wall was further corroborated in 6 of the same samples used for methylation and gene expression analysis with immunohistochemistry. Conclusions: This study is the first to assess DNA methylation in VSMCs from individuals with AAA using NGS, and provides further evidence there is an epigenetic basis to AAA. Our study shows that methylation status of the SMYD2 promoter may be linked with decreased SMYD2 expression in disease pathobiology. In support of our work, downregulated SMYD2 has previously been associated with adverse cardiovascular physiology and inflammation, which are both hallmarks of AAA. The identification of such adverse epigenetic modifications could potentially contribute towards the development of epigenetic treatment strategies in the future.

VariantValidator: Accurate validation, mapping, and formatting of sequence variation descriptions.

The Human Genome Variation Society (HGVS) variant nomenclature is widely used to describe sequence variants in scientific publications, clinical reports, and databases. However, the HGVS recommendations are complex and this often results in inaccurate variant descriptions being reported. The open-source hgvs Python package (https://github.com/biocommons/hgvs) provides a programmatic interface for parsing, manipulating, formatting, and validating of variants according to the HGVS recommendations, but does not provide a user-friendly Web interface. We have developed a Web-based variant validation tool, VariantValidator (https://variantvalidator.org/), which utilizes the hgvs Python package and provides additional functionality to assist users who wish to accurately describe and report sequence-level variations that are compliant with the HGVS recommendations. VariantValidator was designed to ensure that users are guided through the intricacies of the HGVS nomenclature, for example, if the user makes a mistake, VariantValidator automatically corrects the mistake if it can, or provides helpful guidance if it cannot. In addition, VariantValidator has the facility to interconvert genomic variant descriptions in HGVS and Variant Call Format with a degree of accuracy that surpasses most competing solutions.

A mechanistic basis for amplification differences between samples and between genome regions.

BACKGROUND: For many analytical methods the efficiency of DNA amplification varies across the genome and between samples. The most affected genome regions tend to correlate with high C + G content, however this relationship is complex and does not explain why the direction and magnitude of effects varies considerably between samples. RESULTS: Here, we provide evidence that sequence elements that are particularly high in C + G content can remain annealed even when aggressive melting conditions are applied. In turn, this behavior creates broader 'Thermodynamically Ultra-Fastened' (TUF) regions characterized by incomplete denaturation of the two DNA strands, so reducing amplification efficiency throughout these domains. CONCLUSIONS: This model provides a mechanistic explanation for why some genome regions are particularly difficult to amplify and assay in many procedures, and importantly it also explains inter-sample variability of this behavior. That is, DNA samples of varying quality will carry more or fewer nicks and breaks, and hence their intact TUF regions will have different lengths and so be differentially affected by this amplification suppression mechanism - with 'higher' quality DNAs being the most vulnerable. A major practical consequence of this is that inter-region and inter-sample variability can be largely overcome by employing routine fragmentation methods (e.g. sonication or restriction enzyme digestion) prior to sample amplification.