SweGen: a whole-genome data resource of genetic variability in a cross-section of the Swedish population

An integrated approach for genetic risk profiling of typhoid, tuberculosis, and cholera in local population of tehsil Haripur

Despite notable progress in public health throughout the 21st century, infectious diseases like tuberculosis, typhoid, and cholera remain serious threats to global health, particularly in high-risk regions. Understanding the genetic factors that influence susceptibility and resistance to these diseases is essential for developing more effective strategies for their prevention and treatment. This study investigates the genetic variations associated with these infectious diseases with a focus on regions where these diseases are most prevalent. The aim of this study is to identify genetic variants that may influence susceptibility to tuberculosis, typhoid, and cholera. A thorough analysis of genetic variants associated with susceptibility and resistance to tuberculosis, typhoid, and cholera was conducted. Using publicly available genetic data from various ethnic groups. The allele frequency of the prioritized variants was calculated to assess their distribution within the different populations, including Middle Eastern, Ashkenazi Jewish, European (Non-Finnish), Latino/Admixed American, South Asian, East Asian, European (Finnish), North Asian, Southeast Asian, African American, and Swedish populations. The variants of the IL1β gene exhibiting the highest allele frequencies in the South Asian population were identified and subsequently examined within the local population. Specifically, the variants rs1143627 and rs1143629, which demonstrate the highest allele frequencies and are associated with typhoid, tuberculosis, and cholera, were subjected to detailed analysis. To determine their distribution and potential impact on disease susceptibility. In the local population, statistical analysis of the available sample revealed allele frequencies of 0.1128 % for IL1β (rs1143627) and 0.18 % for IL1β (rs1143629). Furthermore, these findings revealed that certain genetic profiles may play a role in the population's overall response to infectious diseases such as tuberculosis, typhoid, and cholera. This research has the potential to guide future public health strategies for more effective management and prevention of these diseases.

UYSD: a novel data repository accessible via public website for worldwide population frequencies of Y-SNP haplogroups

For decades, there has been scientific interest in the variation and geographic distribution of paternal lineages associated with the human Y chromosome. However, the relevant data have been dispersed across numerous publications, making it difficult to consolidate. Additionally, understanding the relationships between different variants, and the tools used to analyze them, have evolved over time, further complicating efforts to harmonize this information. The Universal Y-SNP Database (UYSD) marks a substantial advancement by providing a comprehensive and accessible platform for Y-SNP and haplogroup data from populations around the world. UYSD harmonizes diverse datasets into a unified repository, facilitating the exploration of global Y-chromosomal variation. The platform handles data generated with both high- and low-throughput technology and is compatible with the automated analysis software tool, Yleaf v3. Key functionalities include the ability to: i) visualize haplogroup distributions on an interactive world map, ii) estimate haplogroup frequencies in geographic regions with sparse data through interpolation, and iii) display detailed phylogenetic trees of Y-chromosomal haplogroups. Currently, UYSD encompasses data from over 6,600 males across 27 populations. This dataset largely aligns with known global Y-haplogroup patterns, but also reveals unexplored finer-scale geographic variations. While the present dataset is largely European-centered, UYSD is designed for ongoing expansion by the scientific community, aiming to include more global data and higher-resolution population sequencing data. The platform thus offers valuable insights into human genetic diversity and migration patterns, serving several fields of research such as: human population genetics, genetic anthropology, ancient DNA analysis and forensic genetics.

Egypt Genome: Towards an African new genomic era

BACKGROUND

Studying the human genome is crucial to embrace precision medicine through tailoring treatment and prevention strategies to the unique genetic makeup and molecular information of individuals. After human genome project (1990-2003) generated the first full sequence of a human genome, there have been concerns towards Northern bias due to underrepresentation of other populations. Multiple countries have now established national genome projects aiming at the genomic knowledge that can be harnessed from their populations, which in turn can serve as a basis for their health care policies in the near future.

AIM OF REVIEW

The intention is to introduce the recently established Egypt Genome (EG) to delineate the genomics and genetics of both the modern and Ancient Egyptian populations. Leveraging genomic medicine to improve precision medicine strategies while building a solid foundation for large-scale genomic research capacity is the fundamental focus of EG.

KEY SCIENTIFIC CONCEPTS

EG generated genomic knowledge is predicted to enrich the existing human genome and to expand its diversity by studying the underrepresented African/Middle Eastern populations. The insightful impact of EG goes beyond Egypt and Africa as it fills the knowledge gaps in health and disease genomics towards improved and sustainable genomic-driven healthcare systems for better outcomes. Promoting the integration of genomics into clinical practice and spearheading the implementation of genomic-driven healthcare and precision medicine is therefore a key focus of EG. Mining into the wealth of Ancient Egyptian Genomics to delineate the genetic bridge between the contemporary and Ancient Egyptian populations is another excitingly unique area of EG to realize the global vision of human genome.

Lessons from national biobank projects utilizing whole-genome sequencing for population-scale genomics

Large-scale national biobank projects utilizing whole-genome sequencing have emerged as transformative resources for understanding human genetic variation and its relationship to health and disease. These initiatives, which include the UK Biobank, All of Us Research Program, Singapore's PRECISE, Biobank Japan, and the National Project of Bio-Big Data of Korea, are generating unprecedented volumes of high-resolution genomic data integrated with comprehensive phenotypic, environmental, and clinical information. This review examines the methodologies, contributions, and challenges of major WGS-based national genome projects worldwide. We first discuss the landscape of national biobank initiatives, highlighting their distinct approaches to data collection, participant recruitment, and phenotype characterization. We then introduce recent technological advances that enable efficient processing and analysis of large-scale WGS data, including improvements in variant calling algorithms, innovative methods for creating multi-sample VCFs, optimized data storage formats, and cloud-based computing solutions. The review synthesizes key discoveries from these projects, particularly in identifying expression quantitative trait loci and rare variants associated with complex diseases. Our review introduces the latest findings from the National Project of Bio-Big Data of Korea, which has advanced our understanding of population-specific genetic variation and rare diseases in Korean and East Asian populations. Finally, we discuss future directions and challenges in maximizing the impact of these resources on precision medicine and global health equity. This comprehensive examination demonstrates how large-scale national genome projects are revolutionizing genetic research and healthcare delivery while highlighting the importance of continued investment in diverse, population-specific genomic resources.

Genetic and environmental contributions to individual differences in visual attention and oculomotor control in early infancy

Infants differ in their level of eye movement control, which at the extreme could be linked to autism. We assessed eye movements in 450 twins (225 pairs, 57% monozygotic, 46% female, aged 5-6 months) using the gap-overlap eye-tracking task. Shorter latency in the gap condition was associated with having more parent-rated autistic traits at 2 years. Latency across the task's three conditions was primarily explained by one highly heritable latent factor likely representing individual differences in basic oculomotor efficiency and/or in visual information processing. Additionally, disengagement of attention was linked to unique genetic factors, suggesting that genetic factors involved in visual attention are different from those involved in basic visual information processing and oculomotor efficiency.

Advanced algorithm for step detection in single-entity electrochemistry: a comparative study of wavelet transforms and convolutional neural networks

Single-entity electrochemistry (SEE) is an emerging field within electrochemistry focused on investigating individual entities such as nanoparticles, bacteria, cells, or single molecules. Accurate identification and analysis of SEE signals require effective data processing methods for unbiased and automated feature extraction. In this study, we apply and compare two approaches for step detection in SEE data: discrete wavelet transforms (DWT) and convolutional neural networks (CNN).

An integrated approach to predict genetic risk for Mosquito-Borne diseases in the local Population of Tehsil Haripur, Khyber Pakhtunkhwa, Pakistan

Highly variable response shown by individuals against mosquito-borne infections suggests that host genetic factors play an important role in determining mosquito-borne disease onset. Therefore, it is necessary to determine the genetic risk of these diseases in specific populations. The current study aimed to determine the percentage of individuals in the general population carrying mosquito-borne disease susceptibility and protection-related variants. This study initially aggregated mosquito-borne disease susceptibility and protection-related variants from all publically available data and literature. Afterward, the allele frequency was calculated 1009 genetic variants of 366 genes associated with susceptibility and protection to estimate the global prevalence in multiple ethnicities (Middle Eastern, Ashkenazi Jewish, European (Non-Finnish), Latino/Admixed American, South Asian, East Asian, European (Finnish), North Asian, Southeast Asian, African American, and Swedish population). Furthermore, the cumulative allele frequency of all susceptibility and protection-related variants was calculated in diverse ethnic groups and the relationship with mosquito-borne disease-associated morbidity and mortality was examined to determine whether results are consistent with founder effect in these populations. Two prioritized genetic variants of IL-10 (rs1800871) and FcγRIIA (rs1801274) were examined in the Tehsil Haripur population to assess the genetic risks linked to susceptibility and protection against mosquito-borne diseases. The findings of this study revealed overlapping genes most implicated in mosquito-borne disease linked with susceptibility and protection across different ethnic ancestries. In the available sample size, the percentage of TC and TT genotypes in IL-10 genetic variant (rs1800871) was 12% and 88%, respectively and GA and GG genotypes in FcγRIIA(rs1801274) genetic variant were 6% and 94% respectively. Based on statistical analysis, the percentage allele frequency of IL-10 (rs1800871) variant was 0.2112% and the FcγRIIA (rs1801274) variant is 0.1128% in the current study. Additionally, this study reflects that screening of genetic variants associated with susceptibility and protection in a population gives better insights into organizing public health awareness campaigns to control diseases.

Y chromosome sequencing data suggest dual paths of haplogroup N1a1 into Finland

The paternally inherited Y chromosome is highly informative of genetic ancestry, therefore making it useful in studies of population history. In Finland, two Y-chromosomal haplogroups reveal the major substructure of the population: N1a1 enriched in the northeast and I1a in the southwest, suggested to reflect eastern and western ancestry contributions to the population. Yet, beyond these major Y-chromosomal lineages, the distribution of finer-scale Y-chromosomal variation has not been assessed in Finland. Here, we provide the most comprehensive Y-chromosomal study among the Finns to date, exploiting sequences for 1802 geographically mapped Finnish Y chromosomes from the FINRISK project. We assessed the distribution of common Y-chromosomal haplogroups (frequency ≥1%) throughout 19 Finnish regions and compared the autosomal genetic backgrounds of the Y-chromosomal haplogroups. With such high-resolution data, we were able to find previously unreported sublineages and resolve phylogenetic relationships within haplogroups N1a1 (64%), I1a (25%), R1a (4.3%), and R1b (4.8%). We further find novel geographical enrichment patterns among these Y-chromosomal haplogroups, most notably observed for haplogroup N1a1 dividing into two lineages with differing distributions. While sublineage N-Z1934 (42%) followed a northeastern enrichment pattern observed for all N1a1 carriers in general, sublineage N-VL29 (22%) displayed an enrichment in the southwest. Further, the carriers of N-VL29 showed a higher proportion of southwestern autosomal ancestry compared to carriers of N-Z1934. Collectively, these results point to distinct demographics within haplogroup N1a1, possibly induced by two distinct arrival routes into Finland. Overall, our study suggests a more complex genetic population history for Finns than previously proposed.

Regional Variations in the Incidence of Lichen Sclerosus in Sweden: Insights from a Nationwide Register Study (2001-2021)

Background: Lichen Sclerosus (LSc) is a chronic inflammatory skin condition predominantly affecting the anogenital regions, with a well-recognised potential for malignancy. This study examines the incidence, demographic characteristics, and regional distribution of LSc in Sweden over a 20-year period. The analysis is based on data from the Swedish National Patient Register (NPR), with a focus on cases diagnosed in specialist care settings. Methods: A nationwide register-based study was conducted using data from the NPR, identifying cases of LSc diagnosed between 1 January 2001 and 1 January 2021. Data analysis explored incidence by region, sex, age, and diagnostic care setting. A total of 154,424 patients with LSc were included, and the control group consisted of the general Swedish population without known LSc. Results: The mean annual incidence of LSc was 0.81 per 1000 individuals across Sweden, with higher rates in females (1.14 per 1000) compared to males (0.47 per 1000). Incidence varied significantly across regions, with Blekinge, Kalmar, and Gotland exhibiting the highest rates. This study analysed the distribution of LSc diagnoses across medical specialties, finding that 29.8% of cases were managed by dermatology and venereology, while 17.2% were handled by gynaecology and obstetrics. The analysis of marital status revealed that the proportion of married and divorced LSc patients was significantly lower than the national averages for men and women. Conclusions: This study highlights significant regional variations in LSc incidence. Future research should investigate whether environmental factors, genetic predisposition, socioeconomic disparities, or variations in healthcare access contribute to the variations in incidence. Such insights could lead to more targeted public health strategies for managing LSc across different regions.

Diagnosing Monogenic Stroke at Younger Age

BACKGROUND

An increasing number of monogenic conditions underlying stroke are being identified. We explored the possibilities of increasing the diagnostic yield of monogenic stroke in a population under 56 years of age.

METHODS

Fifty probands ≤55 years at their first stroke episode were characterized clinically and investigated by whole genome sequencing. Probands had one or more of: (1) one or more first to second degree relatives with stroke under 60 years or same stroke-causing condition/disease; (2) no hypertension, hypercholesterolemia, diabetes, heart disease, or smoking; or (3) either multiple stroke episodes or multiple arterial dissections. Variants with minor allele frequency under 0.01, identified by using our stroke gene panels, were assessed. The stroke subtypes, including large artery atherosclerotic, large artery nonatherosclerotic (tortuosity, dolichoectasia, aneurysm, nonatherosclerotic dissection, or occlusion), cerebral small vessel disease, cardioembolic (arrhythmia, heart defect, or cardiomyopathy), coagulation dysfunctions (venous thrombosis, arterial thrombosis, or bleeding tendency), intracerebral hemorrhage, vascular malformations (cavernoma or arteriovenous malformations), metabolic disorders, or cryptogenic embolic, were used for genotype-phenotype correlation. In a final step, we combined genetic and clinical information to determine if the genetic variant likely was the cause of stroke in the patients.

RESULTS

Whole genome sequencing of younger patients with stroke identified 17 clinically matching genetic variants in 15 of 50 (30%) patients, while a stronger clinical correlation with stroke was established in only 6 (12%) of them. Stroke-related genetic variants were identified in 4 of 5 (80%) patients with cardioembolic stroke subtype, 3 of 4 (75%) with intracerebral hemorrhage, 7 of 18 (39%) with cryptogenic embolic stroke, 1 of 6 (17%) with small vessel disease, and 3 of 15 (20%) of patients with nonatherosclerotic large artery stroke, including 1 of 11 (9%) with cervical dissection stroke.

CONCLUSIONS

Careful clinical interpretation of whole genome data using stroke gene panels can detect monogenic causes of early stroke, allowing individualized follow-up and opening new possibilities for potential treatment.

A national long-read sequencing study on chromosomal rearrangements uncovers hidden complexities

Clinical genetic laboratories often require a comprehensive analysis of chromosomal rearrangements/structural variants (SVs), from large events like translocations and inversions to supernumerary ring/marker chromosomes and small deletions or duplications. Understanding the complexity of these events and their clinical consequences requires pinpointing breakpoint junctions and resolving the derivative chromosome structure. This task often surpasses the capabilities of short-read sequencing technologies. In contrast, long-read sequencing techniques present a compelling alternative for clinical diagnostics. Here, Genomic Medicine Sweden-Rare Diseases has explored the utility of HiFi Revio long-read genome sequencing (lrGS) for digital karyotyping of SVs nationwide. The 16 samples from 13 families were collected from all Swedish healthcare regions. Prior investigations had identified 16 SVs, ranging from simple to complex rearrangements, including inversions, translocations, and copy number variants. We have established a national pipeline and a shared variant database for variant calling and filtering. Using lrGS, 14 of the 16 known SVs are detected. Of these, 13 are mapped at nucleotide resolution, and one complex rearrangement is only visible by read depth. Two Chromosome 21 rearrangements, one mosaic, remain undetected. Average read lengths are 8.3-18.8 kb with coverage exceeding 20× for all samples. De novo assembly results in a limited number of phased contigs per individual (N50 6-86 Mb), enabling direct characterization of the chromosomal rearrangements. In a national pilot study, we demonstrate the utility of HiFi Revio lrGS for analyzing chromosomal rearrangements. Based on our results, we propose a 5-year plan to expand lrGS use for rare disease diagnostics in Sweden.

Leveraging the T2T assembly to resolve rare and pathogenic inversions in reference genome gaps

Chromosomal inversions (INVs) are particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage-sensitive genes in cis Short-read genome sequencing (srGS) can only resolve ∼70% of cytogenetically visible inversions referred to clinical diagnostic laboratories, likely due to breakpoints in repetitive regions. Here, we study 12 inversions by long-read genome sequencing (lrGS) (n = 9) or srGS (n = 3) and resolve nine of them. In four cases, the inversion breakpoint region was missing from at least one of the human reference genomes (GRCh37, GRCh38, T2T-CHM13) and a reference agnostic analysis was needed. One of these cases, an INV9 mappable only in de novo assembled lrGS data using T2T-CHM13 disrupts EHMT1 consistent with a Mendelian diagnosis (Kleefstra syndrome 1; MIM#610253). Next, by pairwise comparison between T2T-CHM13, GRCh37, and GRCh38, as well as the chimpanzee and bonobo, we show that hundreds of megabases of sequence are missing from at least one human reference, highlighting that primate genomes contribute to genomic diversity. Aligning population genomic data to these regions indicated that these regions are variable between individuals. Our analysis emphasizes that T2T-CHM13 is necessary to maximize the value of lrGS for optimal inversion detection in clinical diagnostics. These results highlight the importance of leveraging diverse and comprehensive reference genomes to resolve unsolved molecular cases in rare diseases.

GBA1 T369M and Parkinson's disease - Further evidence of a lack of association in the Swedish population

Variants in GBA1 are important genetic risk factors in Parkinson's disease (PD). GBA1 T369M has been linked to an ∼80 % increased PD risk but the reports are conflicting and the relevance of GBA1 variants in different populations varies. A lack of association between T369M and PD in the Swedish population was recently reported but needs further validation. We therefore investigated T369M in 1,808 PD patients and 2,183 controls and our results support that T369M is not a risk factor for PD in the Swedish population.

Genetic determinants of renal scarring in children with febrile UTI

BACKGROUND

Febrile urinary tract infections (UTIs) are among the most severe bacterial infections in infants, in which a subset of patients develops complications. Identifying infants at risk of recurrent infections or kidney damage based on clinical signs is challenging. Previous observations suggest that genetic factors influence UTI outcomes and could serve as predictors of disease severity. In this study, we conducted a nationwide survey of infant genotypes to develop a strategy for infection management based on individual genetic risk. Our aims were to identify genetic susceptibility variants for renal scarring (RS) and genetic host factors predisposing to dilating vesicoureteral reflux (VUR) and recurrent UTIs.

METHODS

To assess genetic susceptibility, we collected and analyzed DNA from blood using exome genotyping. Disease-associated genetic variants were identified through bioinformatics analysis, including allelic frequency tests and odds ratio calculations. Kidney involvement was defined using dimercaptosuccinic acid (DMSA) scintigraphy.

RESULTS

In this investigation, a cohort comprising 1087 infants presenting with their first episode of febrile UTI was included. Among this cohort, a subset of 137 infants who underwent DMSA scanning was subjected to gene association analysis. Remarkable genetic distinctions were observed between patients with RS and those exhibiting resolved kidney involvement. Notably, the genetic signature indicative of renal scarring prominently featured mitochondrial genes.

CONCLUSIONS

In this nationwide study of genetic susceptibility to RS after febrile UTIs in infancy, we identified a profile dominated by mitochondrial polymorphisms. This profile can serve as a predictor of future complications, including RS and recurrent UTIs.

Spinocerebellar ataxia type 4 is caused by a GGC expansion in the ZFHX3 gene and is associated with prominent dysautonomia and motor neuron signs

BACKGROUND

Spinocerebellar ataxia 4 (SCA4), characterized in 1996, features adult-onset ataxia, polyneuropathy, and linkage to chromosome 16q22.1; its underlying mutation has remained elusive.

OBJECTIVE

To explore the radiological and neuropathological abnormalities in the entire neuroaxis in SCA4 and search for its mutation.

METHODS

Three Swedish families with undiagnosed ataxia went through clinical, neurophysiological, and neuroimaging tests, including PET studies and genetic investigations. In four cases, neuropathological assessments of the neuroaxis were performed. Genetic testing included short read whole genome sequencing, short tandem repeat analysis with ExpansionHunter de novo, and long read sequencing.

RESULTS

Novel features for SCA4 include dysautonomia, motor neuron affection, and abnormal eye movements. We found evidence of anticipation; neuroimaging demonstrated atrophy in the cerebellum, brainstem, and spinal cord. [18F]FDG-PET demonstrated brain hypometabolism and [11C]Flumazenil-PET reduced binding in several brain lobes, insula, thalamus, hypothalamus, and cerebellum. Moderate to severe loss of Purkinje cells in the cerebellum and of motor neurons in the anterior horns of the spinal cord along with pronounced degeneration of posterior tracts was also found. Intranuclear, mainly neuronal, inclusions positive for p62 and ubiquitin were sparse but widespread in the CNS. This finding prompted assessment for nucleotide expansions. A polyglycine stretch encoding GGC expansions in the last exon of the zink finger homeobox 3 gene was identified segregating with disease and not found in 1000 controls.

CONCLUSIONS

SCA4 is a neurodegenerative disease caused by a novel GGC expansion in the coding region of ZFHX3, and its spectrum is expanded to include dysautonomia and neuromuscular manifestations.

Uncommon Variants in FLG2 and NOD2 Are Associated with Atopic Dermatitis in the Ethiopian Population

Loss-of-function variants in the FLG gene have been identified as the strongest cause of susceptibility to atopic dermatitis (AD) in Europeans and Asians. However, very little is known about the genetic etiology behind AD in African populations, where the prevalence of AD is notably high. We sought to investigate the genetic origins of AD by performing whole-genome sequencing in an Ethiopian family with 12 individuals and several affected in different generations. We identified 2 variants within FLG2 (p.D13Y) and NOD2 (p.A918S) genes cosegregating with AD in the affected individuals. Further genotyping analyses in both Ethiopian and Swedish AD cases and controls revealed a significant association with the FLG2 variant (p.D13Y, P < .0013) only in the Ethiopian cohort. However, the NOD2 variant (p.A918S) did not show any association in our Ethiopian cohort. Instead, 2 previously recognized NOD2 variants (p.A849V, P < .0085 and p.G908R, P < .0036) were significantly associated with AD in our Ethiopian cohort. Our study indicates that the FLG2 and NOD2 genes might be important in the etiology of AD in Ethiopians. Additional genetic and functional studies are needed to confirm the role of these genes and the associated variants into the development of AD.

Predicting human and viral protein variants affecting COVID-19 susceptibility and repurposing therapeutics

The COVID-19 disease is an ongoing global health concern. Although vaccination provides some protection, people are still susceptible to re-infection. Ostensibly, certain populations or clinical groups may be more vulnerable. Factors causing these differences are unclear and whilst socioeconomic and cultural differences are likely to be important, human genetic factors could influence susceptibility. Experimental studies indicate SARS-CoV-2 uses innate immune suppression as a strategy to speed-up entry and replication into the host cell. Therefore, it is necessary to understand the impact of variants in immunity-associated human proteins on susceptibility to COVID-19. In this work, we analysed missense coding variants in several SARS-CoV-2 proteins and their human protein interactors that could enhance binding affinity to SARS-CoV-2. We curated a dataset of 19 SARS-CoV-2: human protein 3D-complexes, from the experimentally determined structures in the Protein Data Bank and models built using AlphaFold2-multimer, and analysed the impact of missense variants occurring in the protein-protein interface region. We analysed 468 missense variants from human proteins and 212 variants from SARS-CoV-2 proteins and computationally predicted their impacts on binding affinities for the human viral protein complexes. We predicted a total of 26 affinity-enhancing variants from 13 human proteins implicated in increased binding affinity to SARS-CoV-2. These include key-immunity associated genes (TOMM70, ISG15, IFIH1, IFIT2, RPS3, PALS1, NUP98, AXL, ARF6, TRIMM, TRIM25) as well as important spike receptors (KREMEN1, AXL and ACE2). We report both common (e.g., Y13N in IFIH1) and rare variants in these proteins and discuss their likely structural and functional impact, using information on known and predicted functional sites. Potential mechanisms associated with immune suppression implicated by these variants are discussed. Occurrence of certain predicted affinity-enhancing variants should be monitored as they could lead to increased susceptibility and reduced immune response to SARS-CoV-2 infection in individuals/populations carrying them. Our analyses aid in understanding the potential impact of genetic variation in immunity-associated proteins on COVID-19 susceptibility and help guide drug-repurposing strategies.

Mind the gap: the relevance of the genome reference to resolve rare and pathogenic inversions

Both long-read genome sequencing (lrGS) and the recently published Telomere to Telomere (T2T) reference genome provide increased coverage and resolution across repetitive regions promising heightened structural variant detection and improved mapping. Inversions (INV), intrachromosomal segments which are rotated 180° and inserted back into the same chromosome, are a class of structural variants particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage sensitive genes in cis . Here we remapped the genome data from six individuals carrying unsolved cytogenetically detected inversions. An INV6 and INV10 were resolved using GRCh38 and T2T-CHM13. Finally, an INV9 required optical genome mapping, de novo assembly of lrGS data and T2T-CHM13. This inversion disrupted intron 25 of EHMT1, confirming a diagnosis of Kleefstra syndrome 1 (MIM#610253). These three inversions, only mappable in specific references, prompted us to investigate the presence and population frequencies of differential reference regions (DRRs) between T2T-CHM13, GRCh37, GRCh38, the chimpanzee and bonobo, and hundreds of megabases of DRRs were identified. Our results emphasize the significance of the chosen reference genome and the added benefits of lrGS and optical genome mapping in solving rearrangements in challenging regions of the genome. This is particularly important for inversions and may impact clinical diagnostics.

Is GBA1 T369M not a risk factor for Parkinson's disease in the Swedish population?

Variants in GBA1 are important genetic risk factors in Parkinson's disease (PD). GBA1 T369M has been linked to an ~80% increased PD risk but the reports are conflicting and the relevance of GBA1 variants in different populations varies. A lack of association between T369M and PD in the Swedish population was recently reported but needs further validation. We therefore investigated T369M in 1,808 PD patients and 2,183 controls and our results support that T369M is not a risk factor for PD in the Swedish population.

Rare coding variants in NOX4 link high ROS levels to psoriatic arthritis mutilans

Psoriatic arthritis mutilans (PAM) is the rarest and most severe form of psoriatic arthritis, characterized by erosions of the small joints and osteolysis leading to joint disruption. Despite its severity, the underlying mechanisms are unknown, and no susceptibility genes have hitherto been identified. We aimed to investigate the genetic basis of PAM by performing massive parallel sequencing in sixty-one patients from the PAM Nordic cohort. We found rare variants in the NADPH oxidase 4 (NOX4) in four patients. In silico predictions show that the identified variants are potentially damaging. NOXs are the only enzymes producing reactive oxygen species (ROS). NOX4 is specifically involved in the differentiation of osteoclasts, the cells implicated in bone resorption. Functional follow-up studies using cell culture, zebrafish models, and measurement of ROS in patients uncovered that these NOX4 variants increase ROS levels both in vitro and in vivo. We propose NOX4 as the first candidate susceptibility gene for PAM. Our study links high levels of ROS caused by NOX4 variants to the development of PAM, offering a potential therapeutic target.

Gustavson syndrome is caused by an in-frame deletion in RBMX associated with potentially disturbed SH3 domain interactions

RNA binding motif protein X-linked (RBMX) encodes the heterogeneous nuclear ribonucleoprotein G (hnRNP G) that regulates splicing, sister chromatid cohesion and genome stability. RBMX knock down experiments in various model organisms highlight the gene's importance for brain development. Deletion of the RGG/RG motif in hnRNP G has previously been associated with Shashi syndrome, however involvement of other hnRNP G domains in intellectual disability remain unknown. In the current study, we present the underlying genetic and molecular cause of Gustavson syndrome. Gustavson syndrome was first reported in 1993 in a large Swedish five-generation family presented with profound X-linked intellectual disability and an early death. Extensive genomic analyses of the family revealed hemizygosity for a novel in-frame deletion in RBMX in affected individuals (NM_002139.4; c.484_486del, p.(Pro162del)). Carrier females were asymptomatic and presented with skewed X-chromosome inactivation, indicating silencing of the pathogenic allele. Affected individuals presented minor phenotypic overlap with Shashi syndrome, indicating a different disease-causing mechanism. Investigation of the variant effect in a neuronal cell line (SH-SY5Y) revealed differentially expressed genes enriched for transcription factors involved in RNA polymerase II transcription. Prediction tools and a fluorescence polarization assay imply a novel SH3-binding motif of hnRNP G, and potentially a reduced affinity to SH3 domains caused by the deletion. In conclusion, we present a novel in-frame deletion in RBMX segregating with Gustavson syndrome, leading to disturbed RNA polymerase II transcription, and potentially reduced SH3 binding. The results indicate that disruption of different protein domains affects the severity of RBMX-associated intellectual disabilities.

Loss of Cardiac Splicing Regulator RBM20 Is Associated With Early-Onset Atrial Fibrillation

We showed an association between atrial fibrillation and rare loss-of-function (LOF) variants in the cardiac splicing regulator RBM20 in 2 independent cohorts. In a rat model with loss of RBM20, we demonstrated altered splicing of sarcomere genes (NEXN, TTN, TPM1, MYOM1, and LDB3), and differential expression in key cardiac genes. We identified altered sarcomere and mitochondrial structure on electron microscopy imaging and found compromised mitochondrial function. Finally, we demonstrated that 3 novel LOF variants in RBM20, identified in patients with atrial fibrillation, lead to significantly reduced splicing activity. Our results implicate alternative splicing as a novel proarrhythmic mechanism in the atria.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than fivefold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here, we sought to define the roles of PAX1 and newly identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); p=7.07E-11, OR = 1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wild-type. By genetic targeting we found that wild-type Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, the latter suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2 or tamoxifen treatment significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a PAX1-COL11a1-MMP3 signaling axis in spinal chondrocytes.

Integrating a Polygenic Risk Score into a clinical setting would impact risk predictions in familial breast cancer

BACKGROUND

Low-impact genetic variants identified in population-based genetic studies are not routinely measured as part of clinical genetic testing in familial breast cancer (BC). We studied the consequences of integrating an established Polygenic Risk Score (PRS) (BCAC 313, PRS313) into clinical sequencing of women with familial BC in Sweden.

METHODS

We developed an add-on sequencing panel to capture 313 risk variants in addition to the clinical screening of hereditary BC genes. Index patients with no pathogenic variant from 87 families, and 1000 population controls, were included in comparative PRS calculations. Including detailed family history, sequencing results and tumour pathology information, we used BOADICEA (Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm) V.6 to estimate contralateral and lifetime risks without and with PRS313.

RESULTS

Women with BC but no pathogenic variants in hereditary BC genes have a higher PRS313 compared with population controls (mean+0.78 SD, p<3e-9). Implementing PRS313 in the clinical risk estimation before their BC diagnosis would have changed the recommended follow-up in 24%-45% of women.

CONCLUSIONS

Our results show the potential impact of incorporating PRS313 directly in the clinical genomic investigation of women with familial BC.

Oxidative stress-related genetic variation and antioxidant vitamin intake in intact and ruptured abdominal aortic aneurysm: a Swedish population-based retrospective cohort study

AIMS

The aim of this study is to investigate how genetic variations in genes related to oxidative stress, intake of antioxidant vitamins, and any potential interactions between these factors affect the incidence of intact abdominal aortic aneurysm (AAA) and its rupture (rAAA), accounting for sex differences where possible.

METHODS AND RESULTS

The present retrospective cohort study (n = 25 252) uses baseline single-nucleotide polymorphisms (SNPs) and total antioxidant vitamin intake data from the large population-based, Malmö Diet and Cancer Study. Cumulative incidence of intact AAA was 1.6% and of rAAA 0.3% after a median follow-up of 24.3 years. A variant in NOX3 (rs3749930) was associated with higher rAAA risk in males [adjusted hazard ratio (aHR): 2.49; 95% confidence interval (CI): 1.36-4.35] and the overall population (aHR: 1.88; 95% CI: 1.05-3.37). Higher intakes of antioxidant vitamins, riboflavin, and folate were associated with 20% and 19% reduced intact AAA incidence, respectively. Interestingly, the inverse associations between riboflavin and vitamin D intake with intact AAA incidence were stronger in the individuals carrying the NOX3 variant as compared with the wild-type recessive genotype, i.e. by 60% and 66%, respectively (P for interaction < 0.05). Higher riboflavin intake was associated with a 33% male-specific intact AAA risk reduction, while higher intake of vitamin B12 intake was associated with 55% female-specific intact AAA risk increase; both these associations were significantly modified by sex (P for interaction < 0.05).

CONCLUSIONS

Our findings highlight the role of oxidative stress genetic variations and antioxidant vitamin intake in AAA. Although a low AAA/rAAA sample size limited some analyses, especially in females, our findings highlight the need for future randomized controlled trials and mechanistic studies, to explore the potential benefits of antioxidant vitamins while accounting for genetic and sex differences.

Exonic trinucleotide repeat expansions in ZFHX3 cause spinocerebellar ataxia type 4: A poly-glycine disease

Autosomal-dominant ataxia with sensory and autonomic neuropathy is a highly specific combined phenotype that we described in two Swedish kindreds in 2014; its genetic cause had remained unknown. Here, we report the discovery of exonic GGC trinucleotide repeat expansions, encoding poly-glycine, in zinc finger homeobox 3 (ZFHX3) in these families. The expansions were identified in whole-genome datasets within genomic segments that all affected family members shared. Non-expanded alleles carried one or more interruptions within the repeat. We also found ZFHX3 repeat expansions in three additional families, all from the region of Skåne in southern Sweden. Individuals with expanded repeats developed balance and gait disturbances at 15 to 60 years of age and had sensory neuropathy and slow saccades. Anticipation was observed in all families and correlated with different repeat lengths determined through long-read sequencing in two family members. The most severely affected individuals had marked autonomic dysfunction, with severe orthostatism as the most disabling clinical feature. Neuropathology revealed p62-positive intracytoplasmic and intranuclear inclusions in neurons of the central and enteric nervous system, as well as alpha-synuclein positivity. ZFHX3 is located within the 16q22 locus, to which spinocerebellar ataxia type 4 (SCA4) repeatedly had been mapped; the clinical phenotype in our families corresponded well with the unique phenotype described in SCA4, and the original SCA4 kindred originated from Sweden. ZFHX3 has known functions in neuronal development and differentiation n both the central and peripheral nervous system. Our findings demonstrate that SCA4 is caused by repeat expansions in ZFHX3.

Clinical and genetic analyses of a Swedish patient series diagnosed with ataxia

Hereditary ataxia is a heterogeneous group of complex neurological disorders. Next-generation sequencing methods have become a great help in clinical diagnostics, but it may remain challenging to determine if a genetic variant is the cause of the patient's disease. We compiled a consecutive single-center series of 87 patients from 76 families with progressive ataxia of known or unknown etiology. We investigated them clinically and genetically using whole exome or whole genome sequencing. Test methods were selected depending on family history, clinical phenotype, and availability. Genetic results were interpreted based on the American College of Medical Genetics criteria. For high-suspicion variants of uncertain significance, renewed bioinformatical and clinical evaluation was performed to assess the level of pathogenicity. Thirty (39.5%) of the 76 families had received a genetic diagnosis at the end of our study. We present the predominant etiologies of hereditary ataxia in a Swedish patient series. In two families, we established a clinical diagnosis, although the genetic variant was classified as "of uncertain significance" only, and in an additional three families, results are pending. We found a pathogenic variant in one family, but we suspect that it does not explain the complete clinical picture. We conclude that correctly interpreting genetic variants in complex neurogenetic diseases requires genetics and clinical expertise. The neurologist's careful phenotyping remains essential to confirm or reject a diagnosis, also by reassessing clinical findings after a candidate genetic variant is suggested. Collaboration between neurology and clinical genetics and combining clinical and research approaches optimizes diagnostic yield.

Copy number variations and their effect on the plasma proteome

Structural variations, including copy number variations (CNVs), affect around 20 million bases in the human genome and are common causes of rare conditions. CNVs are rarely investigated in complex disease research because most CNVs are not targeted on the genotyping arrays or the reference panels for genetic imputation. In this study, we characterize CNVs in a Swedish cohort (N = 1,021) using short-read whole-genome sequencing (WGS) and use long-read WGS for validation in a subcohort (N = 15), and explore their effect on 438 plasma proteins. We detected 184,182 polymorphic CNVs and identified 15 CNVs to be associated with 16 proteins (P < 8.22×10-10). Of these, 5 CNVs could be perfectly validated using long-read sequencing, including a CNV which was associated with measurements of the osteoclast-associated immunoglobulin-like receptor (OSCAR) and located upstream of OSCAR, a gene important for bone health. Two other CNVs were identified to be clusters of many short repetitive elements and another represented a complex rearrangement including an inversion. Our findings provide insights into the structure of common CNVs and their effects on the plasma proteome, and highlights the importance of investigating common CNVs, also in relation to complex diseases.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than five-fold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here we sought to define the roles of PAX1 and newly-identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); P=7.07e-11, OR=1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc (IVD)-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wildtype. By genetic targeting we found that wildtype Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, this suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2, or tamoxifen treatment, significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a Pax1-Col11a1-Mmp3 signaling axis in spinal chondrocytes.

Position statement of the International Society for Gastrointestinal Hereditary Tumours (InSiGHT) on APC I1307K and cancer risk

While constitutional pathogenic variants in the APC gene cause familial adenomatous polyposis, APC c.3920T>A; p.Ile1307Lys (I1307K) has been associated with a moderate increased risk of colorectal cancer (CRC), particularly in individuals of Ashkenazi Jewish descent. However, published data include relatively small sample sizes, generating inconclusive results regarding cancer risk, particularly in non-Ashkenazi populations. This has led to different country/continental-specific guidelines regarding genetic testing, clinical management and surveillance recommendations for I1307K. A multidisciplinary international expert group endorsed by the International Society for Gastrointestinal Hereditary Tumours (InSiGHT), has generated a position statement on the APC I1307K allele and its association with cancer predisposition. Based on a systematic review and meta-analysis of the evidence published, the aim of this document is to summarise the prevalence of the APC I1307K allele and analysed the evidence of the associated cancer risk in different populations. Here we provide recommendations on the laboratory classification of the variant, define the role of predictive testing for I1307K, suggest recommendations for cancer screening in I1307K heterozygous and homozygous individuals and identify knowledge gaps to be addressed in future research studies. Briefly, I1307K, classified as pathogenic, low penetrance, is a risk factor for CRC in individuals of Ashkenazi Jewish origin and should be tested in this population, offering carriers specific clinical surveillance. There is not enough evidence to support an increased risk of cancer in other populations/subpopulations. Therefore, until/unless future evidence indicates otherwise, individuals of non-Ashkenazi Jewish descent harbouring I1307K should be enrolled in national CRC screening programmes for average-risk individuals.

Complete Mitochondrial DNA Genome Variation in the Swedish Population

The development of complete mitochondrial genome (mitogenome) reference data for inclusion in publicly available population databases is currently underway, and the generation of more high-quality mitogenomes will only enhance the statistical power of this forensically useful locus. To characterize mitogenome variation in Sweden, the mitochondrial DNA (mtDNA) reads from the SweGen whole genome sequencing (WGS) dataset were analyzed. To overcome the interference from low-frequency nuclear mtDNA segments (NUMTs), a 10% variant frequency threshold was applied for the analysis. In total, 934 forensic-quality mitogenome haplotypes were characterized. Almost 45% of the SweGen haplotypes belonged to haplogroup H. Nearly all mitogenome haplotypes (99.1%) were assigned to European haplogroups, which was expected based on previous mtDNA studies of the Swedish population. There were signature northern Swedish and Finnish haplogroups observed in the dataset (e.g., U5b1, W1a), consistent with the nuclear DNA analyses of the SweGen data. The complete mitogenome analysis resulted in high haplotype diversity (0.9996) with a random match probability of 0.15%. Overall, the SweGen mitogenomes provide a large mtDNA reference dataset for the Swedish population and also contribute to the effort to estimate global mitogenome haplotype frequencies.

Origin matters: Using a local reference genome improves measures in population genomics

Genome sequencing enables answering fundamental questions about the genetic basis of adaptation, population structure and epigenetic mechanisms. Yet, we usually need a suitable reference genome for mapping population-level resequencing data. In some model systems, multiple reference genomes are available, giving the challenging task of determining which reference genome best suits the data. Here, we compared the use of two different reference genomes for the three-spined stickleback (Gasterosteus aculeatus), one novel genome derived from a European gynogenetic individual and the published reference genome of a North American individual. Specifically, we investigated the impact of using a local reference versus one generated from a distinct lineage on several common population genomics analyses. Through mapping genome resequencing data of 60 sticklebacks from across Europe and North America, we demonstrate that genetic distance among samples and the reference genomes impacts downstream analyses. Using a local reference genome increased mapping efficiency and genotyping accuracy, effectively retaining more and better data. Despite comparable distributions of the metrics generated across the genome using SNP data (i.e. π, Tajima's D and FST ), window-based statistics using different references resulted in different outlier genes and enriched gene functions. A marker-based analysis of DNA methylation distributions had a comparably high overlap in outlier genes and functions, yet with distinct differences depending on the reference genome. Overall, our results highlight how using a local reference genome decreases reference bias to increase confidence in downstream analyses of the data. Such results have significant implications in all reference-genome-based population genomic analyses.

Precision medicine in complex diseases-Molecular subgrouping for improved prediction and treatment stratification

Complex diseases are caused by a combination of genetic, lifestyle, and environmental factors and comprise common noncommunicable diseases, including allergies, cardiovascular disease, and psychiatric and metabolic disorders. More than 25% of Europeans suffer from a complex disease, and together these diseases account for 70% of all deaths. The use of genomic, molecular, or imaging data to develop accurate diagnostic tools for treatment recommendations and preventive strategies, and for disease prognosis and prediction, is an important step toward precision medicine. However, for complex diseases, precision medicine is associated with several challenges. There is a significant heterogeneity between patients of a specific disease-both with regards to symptoms and underlying causal mechanisms-and the number of underlying genetic and nongenetic risk factors is often high. Here, we summarize precision medicine approaches for complex diseases and highlight the current breakthroughs as well as the challenges. We conclude that genomic-based precision medicine has been used mainly for patients with highly penetrant monogenic disease forms, such as cardiomyopathies. However, for most complex diseases-including psychiatric disorders and allergies-available polygenic risk scores are more probabilistic than deterministic and have not yet been validated for clinical utility. However, subclassifying patients of a specific disease into discrete homogenous subtypes based on molecular or phenotypic data is a promising strategy for improving diagnosis, prediction, treatment, prevention, and prognosis. The availability of high-throughput molecular technologies, together with large collections of health data and novel data-driven approaches, offers promise toward improved individual health through precision medicine.

Aagenaes syndrome/lymphedema cholestasis syndrome 1 is caused by a founder variant in the 5'-untranslated region of UNC45A

BACKGROUND & AIMS

Lymphedema cholestasis syndrome 1 or Aagenaes syndrome is a condition characterized by neonatal cholestasis, lymphedema, and giant cell hepatitis. The genetic background of this autosomal recessive disease was unknown up to now.

METHODS

A total of 26 patients with Aagenaes syndrome and 17 parents were investigated with whole-genome sequencing and/or Sanger sequencing. PCR and western blot analyses were used to assess levels of mRNA and protein, respectively. CRISPR/Cas9 was used to generate the variant in HEK293T cells. Light microscopy, transmission electron microscopy and immunohistochemistry for biliary transport proteins were performed in liver biopsies.

RESULTS

One specific variant (c.-98G>T) in the 5'-untranslated region of Unc-45 myosin chaperone A (UNC45A) was identified in all tested patients with Aagenaes syndrome. Nineteen were homozygous for the c.-98G>T variant and seven were compound heterozygous for the variant in the 5'-untranslated region and an exonic loss-of-function variant in UNC45A. Patients with Aagenaes syndrome exhibited lower expression of UNC45A mRNA and protein than controls, and this was reproduced in a CRISPR/Cas9-created cell model. Liver biopsies from the neonatal period demonstrated cholestasis, paucity of bile ducts and pronounced formation of multinucleated giant cells. Immunohistochemistry revealed mislocalization of the hepatobiliary transport proteins BSEP (bile salt export pump) and MRP2 (multidrug resistance-associated protein 2).

CONCLUSIONS

c.-98G>T in the 5'-untranslated region of UNC45A is the causative genetic variant in Aagenaes syndrome.

IMPACT AND IMPLICATIONS

The genetic background of Aagenaes syndrome, a disease presenting with cholestasis and lymphedema in childhood, was unknown until now. A variant in the 5'-untranslated region of the Unc-45 myosin chaperone A (UNC45A) was identified in all tested patients with Aagenaes syndrome, providing evidence of the genetic background of the disease. Identification of the genetic background provides a tool for diagnosis of patients with Aagenaes syndrome before lymphedema is evident.

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.

Transposable element insertions in 1000 Swedish individuals

The majority of rare diseases are genetic, and regardless of advanced high-throughput genomics-based investigations, 60% of patients remain undiagnosed. A major factor limiting our ability to identify disease-causing alterations is a poor understanding of the morbid and normal human genome. A major genomic contributor of which function and distribution remain largely unstudied are the transposable elements (TE), which constitute 50% of our genome. Here we aim to resolve this knowledge gap and increase the diagnostic yield of rare disease patients investigated with clinical genome sequencing. To this end we characterized TE insertions in 1000 Swedish individuals from the SweGen dataset and 2504 individuals from the 1000 Genomes Project (1KGP), creating seven population-specific TE insertion databases. Of note, 66% of TE insertions in SweGen were present at >1% in the 1KGP databases, proving that most insertions are common across populations. Focusing on the rare TE insertions, we show that even though ~0.7% of those insertions affect protein coding genes, they rarely affect known disease casing genes (<0.1%). Finally, we applied a TE insertion identification workflow on two clinical cases where disease causing TE insertions were suspected and could verify the presence of pathogenic TE insertions in both. Altogether we demonstrate the importance of TE insertion detection and highlight possible clinical implications in rare disease diagnostics.

Discovery of non-reference processed pseudogenes in the Swedish population

The vast majority of the human genome is non-coding. There is a diversity of non-coding features, some of which have functional importance. Although the non-coding regions constitute the majority of the genome, they remain understudied, and for a long time, these regions have been referred to as junk DNA. Pseudogenes are one of these features. A pseudogene is a non-functional copy of a protein-coding gene. Pseudogenes may arise through a variety of genetic mechanisms. Processed pseudogenes are formed through reverse transcription of mRNA by LINE elements, after which the cDNA is integrated into the genome. Processed pseudogenes are known to be variable across populations; however, the variability and distribution remains unknown. Herein, we apply a custom-designed processed pseudogene pipeline on the whole genome sequencing data of 3,500 individuals; 2,500 individuals from the thousand genomes dataset, as well as 1,000 Swedish individuals. Through these analyses, we discover over 3,000 pseudogenes missing from the GRCh38 reference. Utilising our pipeline, we position 74% of the detected processed pseudogenes-allowing for analyses of formation. Notably, we find that common structural variant callers, such as Delly, classify the processed pseudogenes as deletion events, which are later predicted to be truncating variants. By compiling lists of non-reference processed pseudogenes and their frequencies, we find a great variability of pseudogenes; indicating that non-reference processed pseudogenes may be useful for DNA testing and as population-specific markers. In summary, our findings highlight a great diversity of processed pseudogenes, that processed pseudogenes are actively formed in the human genome; and that our pipeline may be used to reduce false positive structural variation caused by the misalignment and subsequent misclassification of non-reference processed pseudogenes.

Applying Unique Molecular Indices with an Extensive All-in-One Forensic SNP Panel for Improved Genotype Accuracy and Sensitivity

One of the major challenges in forensic genetics is being able to detect very small amounts of DNA. Massively parallel sequencing (MPS) enables sensitive detection; however, genotype errors may exist and could interfere with the interpretation. Common errors in MPS-based analysis are often induced during PCR or sequencing. Unique molecular indices (UMIs) are short random nucleotide sequences ligated to each template molecule prior to amplification. Applying UMIs can improve the limit of detection by enabling accurate counting of initial template molecules and removal of erroneous data. In this study, we applied the FORCE panel, which includes ~5500 SNPs, with a QIAseq Targeted DNA Custom Panel (Qiagen), including UMIs. Our main objective was to investigate whether UMIs can enhance the sensitivity and accuracy of forensic genotyping and to evaluate the overall assay performance. We analyzed the data both with and without the UMI information, and the results showed that both genotype accuracy and sensitivity were improved when applying UMIs. The results showed very high genotype accuracies (>99%) for both reference DNA and challenging samples, down to 125 pg. To conclude, we show successful assay performance for several forensic applications and improvements in forensic genotyping when applying UMIs.

Exome sequencing in a Swedish family with PMS2 mutation with varying penetrance of colorectal cancer: investigating the presence of genetic risk modifiers in colorectal cancer risk

OBJECTIVE

Lynch syndrome is caused by germline mutations in the mismatch repair (MMR) genes, such as the PMS2 gene, and is characterised by a familial accumulation of colorectal cancer. The penetrance of cancer in PMS2 carriers is still not fully elucidated as a colorectal cancer risk has been shown to vary between PMS2 carriers, suggesting the presence of risk modifiers.

METHODS

Whole exome sequencing was performed in a Swedish family carrying a PMS2 missense mutation [c.2113G>A, p.(Glu705Lys)]. Thirteen genetic sequence variants were further selected and analysed in a case-control study (724 cases and 711 controls).

RESULTS

The most interesting variant was an 18 bp deletion in gene BAG1. BAG1 has been linked to colorectal tumour progression with poor prognosis and is thought to promote colorectal tumour cell survival through increased NF-κB activity.

CONCLUSIONS

We conclude the genetic architecture behind the incomplete penetrance of PMS2 is complicated and must be assessed in a genome wide manner using large families and multifactorial analysis.

Heritable defects in telomere and mitotic function selectively predispose to sarcomas

Cancer genetics has to date focused on epithelial malignancies, identifying multiple histotype-specific pathways underlying cancer susceptibility. Sarcomas are rare malignancies predominantly derived from embryonic mesoderm. To identify pathways specific to mesenchymal cancers, we performed whole-genome germline sequencing on 1644 sporadic cases and 3205 matched healthy elderly controls. Using an extreme phenotype design, a combined rare-variant burden and ontologic analysis identified two sarcoma-specific pathways involved in mitotic and telomere functions. Variants in centrosome genes are linked to malignant peripheral nerve sheath and gastrointestinal stromal tumors, whereas heritable defects in the shelterin complex link susceptibility to sarcoma, melanoma, and thyroid cancers. These studies indicate a specific role for heritable defects in mitotic and telomere biology in risk of sarcomas.

Using coding and non-coding rare variants to target candidate genes in patients with severe tinnitus

Tinnitus is the phantom percept of an internal non-verbal set of noises and tones. It is reported by 15% of the population and it is usually associated with hearing and/or brain disorders. The role of structural variants (SVs) in coding and non-coding regions has not been investigated in patients with severe tinnitus. In this study, we performed whole-genome sequencing in 97 unrelated Swedish individuals with chronic tinnitus (TIGER cohort). Rare single nucleotide variants (SNV), large structural variants (LSV), and copy number variations (CNV) were retrieved to perform a gene enrichment analysis in TIGER and in a subgroup of patients with severe tinnitus (SEVTIN, n = 34), according to the tinnitus handicap inventory (THI) scores. An independent exome sequencing dataset of 147 Swedish tinnitus patients was used as a replication cohort (JAGUAR cohort) and population-specific datasets from Sweden (SweGen) and Non-Finish Europeans (NFE) from gnomAD were used as control groups. SEVTIN patients showed a higher prevalence of hyperacusis, hearing loss, and anxiety when they were compared to individuals in the TIGER cohort. We found an enrichment of rare missense variants in 6 and 8 high-constraint genes in SEVTIN and TIGER cohorts, respectively. Of note, an enrichment of missense variants was found in the CACNA1E gene in both SEVTIN and TIGER. We replicated the burden of missense variants in 9 high-constrained genes in the JAGUAR cohort, including the gene NAV2, when data were compared with NFE. Moreover, LSVs in constrained regions overlapping CACNA1E, NAV2, and TMEM132D genes were observed in TIGER and SEVTIN.

Strong Association of Combined Genetic Deficiencies in the Classical Complement Pathway With Risk of Systemic Lupus Erythematosus and Primary Sjögren's Syndrome

OBJECTIVE

Complete genetic deficiency of the complement component C2 is a strong risk factor for monogenic systemic lupus erythematosus (SLE), but whether heterozygous C2 deficiency adds to the risk of SLE or primary Sjögren's syndrome (SS) has not been studied systematically. This study was undertaken to investigate potential associations of heterozygous C2 deficiency and C4 copy number variation with clinical manifestations in patients with SLE and patients with primary SS.

METHODS

The presence of the common 28-bp C2 deletion rs9332736 and C4 copy number variation was examined in Scandinavian patients who had received a diagnosis of SLE (n = 958) or primary SS (n = 911) and in 2,262 healthy controls through the use of DNA sequencing. The concentration of complement proteins in plasma and classical complement function were analyzed in a subgroup of SLE patients.

RESULTS

Heterozygous C2 deficiency-when present in combination with a low C4A copy number-substantially increased the risk of SLE (odds ratio [OR] 10.2 [95% confidence interval (95% CI) 3.5-37.0]) and the risk of primary SS (OR 13.0 [95% CI 4.5-48.4]) when compared to individuals with 2 C4A copies and normal C2. For patients heterozygous for rs9332736 with 1 C4A copy, the median age at diagnosis was 7 years earlier in patients with SLE and 12 years earlier in patients with primary SS when compared to patients with normal C2. Reduced C2 levels in plasma (P = 2 × 10-9 ) and impaired function of the classical complement pathway (P = 0.03) were detected in SLE patients with heterozygous C2 deficiency. Finally, in a primary SS patient homozygous for C2 deficiency, we observed low levels of anti-Scl-70, which suggests a risk of developing systemic sclerosis or potential overlap between primary SS and other systemic autoimmune diseases.

CONCLUSION

We demonstrate that a genetic pattern involving partial deficiencies of C2 and C4A in the classical complement pathway is a strong risk factor for SLE and for primary SS. Our results emphasize the central role of the complement system in the pathogenesis of both SLE and primary SS.

Complex genomic rearrangements: an underestimated cause of rare diseases

Complex genomic rearrangements (CGRs) are known contributors to disease but are often missed during routine genetic screening. Identifying CGRs requires (i) identifying copy number variants (CNVs) concurrently with inversions, (ii) phasing multiple breakpoint junctions incis, as well as (iii) detecting and resolving structural variants (SVs) within repeats. We demonstrate how combining cytogenetics and new sequencing methodologies is being successfully applied to gain insights into the genomic architecture of CGRs. In addition, we review CGR patterns and molecular features revealed by studying constitutional genomic disorders. These data offer invaluable lessons to individuals interested in investigating CGRs, evaluating their clinical relevance and frequency, as well as assessing their impact(s) on rare genetic diseases.

Genome interpretation using in silico predictors of variant impact

Estimating the effects of variants found in disease driver genes opens the door to personalized therapeutic opportunities. Clinical associations and laboratory experiments can only characterize a tiny fraction of all the available variants, leaving the majority as variants of unknown significance (VUS). In silico methods bridge this gap by providing instant estimates on a large scale, most often based on the numerous genetic differences between species. Despite concerns that these methods may lack reliability in individual subjects, their numerous practical applications over cohorts suggest they are already helpful and have a role to play in genome interpretation when used at the proper scale and context. In this review, we aim to gain insights into the training and validation of these variant effect predicting methods and illustrate representative types of experimental and clinical applications. Objective performance assessments using various datasets that are not yet published indicate the strengths and limitations of each method. These show that cautious use of in silico variant impact predictors is essential for addressing genome interpretation challenges.

Identification and functional characterization of a novel susceptibility locus for small vessel vasculitis with MPO-ANCA

OBJECTIVE

To identify and characterize genetic loci associated with the risk of developing ANCA-associated vasculitides (AAV).

METHODS

Genetic association analyses were performed after Illumina sequencing of 1853 genes and subsequent replication with genotyping of selected single nucleotide polymorphisms in a total cohort of 1110 Scandinavian cases with granulomatosis with polyangiitis or microscopic polyangiitis, and 1589 controls. A novel AAV-associated single nucleotide polymorphism was analysed for allele-specific effects on gene expression using luciferase reporter assay.

RESULTS

PR3-ANCA+ AAV was significantly associated with two independent loci in the HLA-DPB1/HLA-DPA1 region [rs1042335, P = 6.3 × 10-61, odds ratio (OR) 0.10; rs9277341, P = 1.5 × 10-44, OR 0.22] and with rs28929474 in the SERPINA1 gene (P = 2.7 × 10-10, OR 2.9). MPO-ANCA+ AAV was significantly associated with the HLA-DQB1/HLA-DQA2 locus (rs9274619, P = 5.4 × 10-25, OR 3.7) and with a rare variant in the BACH2 gene (rs78275221, P = 7.9 × 10-7, OR 3.0), the latter a novel susceptibility locus for MPO-ANCA+ granulomatosis with polyangiitis/microscopic polyangiitis. The rs78275221-A risk allele reduced luciferase gene expression in endothelial cells, specifically, as compared with the non-risk allele.

CONCLUSION

We identified a novel susceptibility locus for MPO-ANCA+ AAV and propose that the associated variant is of mechanistic importance, exerting a regulatory function on gene expression in specific cell types.

Unravelling genetic variants of a swedish family with high risk of prostate cancer

BACKGROUND

Prostate cancer is the most prevalent cancer in men worldwide. It is a polygenic disease with a substantial proportion of heritability. Identification of novel candidate biomarkers is crucial for clinical cancer prevention and the development of therapeutic strategies. Here, we describe the analysis of rare and common genetic variants that can predispose to the development of prostate cancer.

METHODS

Whole-genome sequencing was performed on germline DNA of five Swedish siblings which were diagnosed with prostate cancer. The high-risk variants were identified setting the minor allele frequency < 0.01, CADD > 10 and if tested in PRACTICAL, OR > 1.5, while the low-risk variants were identified minor allele frequency > 0.01, CADD > 10 and if tested in PRACTICAL, OR > 1.1.

RESULTS

We identified 38 candidate high-risk gene variants and 332 candidate low-risk gene variants, where 2 and 14 variants were in coding regions, respectively, that were shared by the brothers with prostate cancer.

CONCLUSIONS

This study expanded the knowledge of potential risk factor candidates involved in hereditary and familial prostate cancer. Our findings can be beneficial when applying targeted screening in families with a high risk of developing the disease.

Mosaic deletions of known genes explain skeletal dysplasias with high and low bone mass

Mosaicism, a state in which an individual has two or more genetically distinct populations of cells in the body, can be difficult to detect because of either mild or atypical clinical presentation and limitations in the commonly used detection methods. Knowledge of the role of mosaicism is limited in many skeletal disorders, including osteopathia striata with cranial sclerosis (OSCS) and cleidocranial dysplasia (CCD). We used whole‐genome sequencing (WGS) with coverage >40× to identify the genetic causes of disease in two clinically diagnosed patients. In a female patient with OSCS, we identified a mosaic 7‐nucleotide frameshift deletion in exon 2 of AMER1, NM_152424.4:c.855_861del:p.(His285Glnfs*7), affecting 8.3% of the WGS reads. In a male patient with CCD, approximately 34% of the WGS reads harbored a 3710‐basepair mosaic deletion, NC_000006.11:g.45514471_45518181del, starting in intron 8 of RUNX2 and terminating in the 3′ untranslated region. Droplet digital polymerase chain reaction was used to validate these deletions and quantify the absolute level of mosaicism in each patient. Although constitutional variants in AMER1 and RUNX2 are a known cause of OSCS and CCD, respectively, the mosaic changes here reported have not been described previously. Our study indicates that mosaicism should be considered in unsolved cases of skeletal dysplasia and should be investigated with comprehensive and sensitive detection methods. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Contribution of rare whole-genome sequencing variants to plasma protein levels and the missing heritability

Despite the success of genome-wide association studies, much of the genetic contribution to complex traits remains unexplained. Here, we analyse high coverage whole-genome sequencing data, to evaluate the contribution of rare genetic variants to 414 plasma proteins. The frequency distribution of genetic variants is skewed towards the rare spectrum, and damaging variants are more often rare. We estimate that less than 4.3% of the narrow-sense heritability is expected to be explained by rare variants in our cohort. Using a gene-based approach, we identify Cis-associations for 237 of the proteins, which is slightly more compared to a GWAS (N = 213), and we identify 34 associated loci in Trans. Several associations are driven by rare variants, which have larger effects, on average. We therefore conclude that rare variants could be of importance for precision medicine applications, but have a more limited contribution to the missing heritability of complex diseases.

Despite the success of genome-wide association studies, much of the genetic contribution to complex traits remains unexplained. Here, the authors identify effects by rare variants on plasma proteins, and estimate the contribution of rare variants to the heritability.

Whole-Exome Sequencing of Germline Variants in Non-BRCA Families with Hereditary Breast Cancer

Breast cancer is the most prevalent malignancy among women worldwide and hereditary breast cancer (HBC) accounts for about 5−10% of the cases. Today, the most recurrent genes known are BRCA1 and BRCA2, accounting for around 25% of familial cases. Although thousands of loss-of-function variants in more than twenty predisposing genes have been found, the majority of familial cases of HBC remain unexplained. The aim of this study was to identify new predisposing genes for HBC in three non-BRCA families with autosomal dominant inheritance pattern using whole-exome sequencing and functional prediction tools. No pathogenic variants in known hereditary cancer-related genes could explain the breast cancer susceptibility in these families. Among 2122 exonic variants with maximum minor allele frequency (MMAF) < 0.1%, between 17−35 variants with combined annotation-dependent depletion (CADD) > 20 segregated with disease in the three analyzed families. Selected candidate genes, i.e., UBASH3A, MYH13, UTP11L, and PAX7, were further evaluated using protein expression analysis but no alterations of cancer-related pathways were observed. In conclusion, identification of new high-risk cancer genes using whole-exome sequencing has been more challenging than initially anticipated, in spite of selected families with pronounced family history of breast cancer. A combination of low- and intermediate-genetic-risk variants may instead contribute the breast cancer susceptibility in these families.

National Genome Initiatives in Europe and the United Kingdom in the Era of Whole-Genome Sequencing: A Comprehensive Review

Identification of genomic variability in population plays an important role in the clinical diagnostics of human genetic diseases. Thanks to rapid technological development in the field of massive parallel sequencing technologies, also known as next-generation sequencing (NGS), complex genomic analyses are now easier and cheaper than ever before, which consequently leads to more effective utilization of these techniques in clinical practice. However, interpretation of data from NGS is still challenging due to several issues caused by natural variability of DNA sequences in human populations. Therefore, development and realization of projects focused on description of genetic variability of local population (often called "national or digital genome") with a NGS technique is one of the best approaches to address this problem. The next step of the process is to share such data via publicly available databases. Such databases are important for the interpretation of variants with unknown significance or (likely) pathogenic variants in rare diseases or cancer or generally for identification of pathological variants in a patient's genome. In this paper, we have compiled an overview of published results of local genome sequencing projects from United Kingdom and Europe together with future plans and perspectives for newly announced ones.