Structural diversity and African origin of the 17q21.31 inversion polymorphism

Genetic diversity of North African populations in the 17q21 genomic region

The demographic history of human populations in North Africa has been characterized by complex migration processes that have determined the current genetic structure of these populations. We examined the autosomal markers of eight sampled populations in northern Africa (Tunisia and Libya) to explore their genetic structure and to place them in a global context. We genotyped a set of 30 autosomal single-nucleotide polymorphisms (SNPs) extending 9.5 Mb and encompassing the 17q21 inversion region. Our data include 403 individuals from Tunisia and Libya. To put our populations in the global context, we analyzed our data in comparison with other populations, including those of the 1000 Genomes Project. To evaluate the data, we conducted genetic diversity, principal component, STRUCTURE, and haplotype analyses. The analysis of genetic composition revealed the genetic heterogeneity of North African populations. The principal component and STRUCTURE analyses converged and revealed the intermediate position of North Africans between Europeans and Asians. Haplotypic analysis demonstrated that the normal (H1) and inverted (H2) polymorphisms in the chromosome 17q21 region occur in North Africa at frequencies similar to those found in European and Southwest Asian populations. The results highlight the complex demographic history of North Africa, reflecting the influence of genetic flow from Europe and the Near East that dates to the prehistoric period. These gene flows added to demographic factors (inbreeding, endogamy), natural factors (topography, Sahara), and cultural factors that play a role in the emergence of the diverse and heterogeneous genetic structures of North African populations. This study contributes to a better understanding of the complex structure of North African populations.

Exploring gene content with pangene graphs

MOTIVATION

The gene content regulates the biology of an organism. It varies between species and between individuals of the same species. Although tools have been developed to identify gene content changes in bacterial genomes, none is applicable to collections of large eukaryotic genomes such as the human pangenome.

RESULTS

We developed pangene, a computational tool to identify gene orientation, gene order and gene copy-number changes in a collection of genomes. Pangene aligns a set of input protein sequences to the genomes, resolves redundancies between protein sequences and constructs a gene graph with each genome represented as a walk in the graph. It additionally finds subgraphs, which we call bibubbles, that capture gene content changes. Applied to the human pangenome, pangene identifies known gene-level variations and reveals complex haplotypes that are not well studied before. Pangene also works with high-quality bacterial pangenome and reports similar numbers of core and accessory genes in comparison to existing tools.

AVAILABILITY AND IMPLEMENTATION

Source code at https://github.com/lh3/pangene; pre-built pangene graphs can be downloaded from https://zenodo.org/records/8118576 and visualized at https://pangene.bioinweb.org.

Structural polymorphism and diversity of human segmental duplications

Segmental duplications (SDs) contribute significantly to human disease, evolution, and diversity yet have been difficult to resolve at the sequence level. We present a population genetics survey of SDs by analyzing 170 human genome assemblies where the majority of SDs are fully resolved using long-read sequence assembly. Excluding the acrocentric short arms, we identify 173.2 Mbp of duplicated sequence (47.4 Mbp not present in the telomere-to-telomere reference) distinguishing fixed from structurally polymorphic events. We find that intrachromosomal SDs are among the most variable with rare events mapping near their progenitor sequences. African genomes harbor significantly more intrachromosomal SDs and are more likely to have recently duplicated gene families with higher copy number when compared to non-African samples. A comparison to a resource of 563 million full-length Iso-Seq reads identifies 201 novel, potentially protein-coding genes corresponding to these copy number polymorphic SDs.

Exploring gene content with pangene graphs

Motivation

The gene content regulates the biology of an organism. It varies between species and between individuals of the same species. Although tools have been developed to identify gene content changes in bacterial genomes, none is applicable to collections of large eukaryotic genomes such as the human pangenome.

Results

We developed pangene, a computational tool to identify gene orientation, gene order and gene copy-number changes in a collection of genomes. Pangene aligns a set of input protein sequences to the genomes, resolves redundancies between protein sequences and constructs a gene graph with each genome represented as a walk in the graph. It additionally finds subgraphs, which we call bibubbles, that capture gene content changes. Applied to the human pangenome, pangene identifies known gene-level variations and reveals complex haplotypes that are not well studied before. Pangene also works with high-quality bacterial pangenome and reports similar numbers of core and accessory genes in comparison to existing tools.

Availability and implementation

Source code at https://github.com/lh3/pangene; pre-built pangene graphs can be downloaded from https://zenodo.org/records/8118576 and visualized at https://pangene.bioinweb.org.

Unraveling the complex role of MAPT-containing H1 and H2 haplotypes in neurodegenerative diseases

A ~ 1 Mb inversion polymorphism exists within the 17q21.31 locus of the human genome as direct (H1) and inverted (H2) haplotype clades. This inversion region demonstrates high linkage disequilibrium, but the frequency of each haplotype differs across ancestries. While the H1 haplotype exists in all populations and shows a normal pattern of genetic variability and recombination, the H2 haplotype is enriched in European ancestry populations, is less frequent in African ancestry populations, and nearly absent in East Asian ancestry populations. H1 is a known risk factor for several neurodegenerative diseases, and has been associated with many other traits, suggesting its importance in cellular phenotypes of the brain and entire body. Conversely, H2 is protective for these diseases, but is associated with predisposition to recurrent microdeletion syndromes and neurodevelopmental disorders such as autism. Many single nucleotide variants and copy number variants define H1/H2 haplotypes and sub-haplotypes, but identifying the causal variant(s) for specific diseases and phenotypes is complex due to the extended linkage equilibrium. In this review, we assess the current knowledge of this inversion region regarding genomic structure, gene expression, cellular phenotypes, and disease association. We discuss recent discoveries and challenges, evaluate gaps in knowledge, and highlight the importance of understanding the effect of the 17q21.31 haplotypes to promote advances in precision medicine and drug discovery for several diseases.

Multiple paralogues and recombination mechanisms drive the high incidence of 22q11.2 Deletion Syndrome

The 22q11.2 deletion syndrome (22q11.2DS) is the most common microdeletion disorder. Why the incidence of 22q11.2DS is much greater than that of other genomic disorders remains unknown. Short read sequencing cannot resolve the complex segmental duplicon structure to provide direct confirmation of the hypothesis that the rearrangements are caused by non-allelic homologous recombination between the low copy repeats on chromosome 22 (LCR22s). To enable haplotype-specific assembly and rearrangement mapping in LCR22 regions, we combined fiber-FISH optical mapping with whole genome (ultra-)long read sequencing or rearrangement-specific long-range PCR on 24 duos (22q11.2DS patient and parent-of-origin) comprising several different LCR22-mediated rearrangements. Unexpectedly, we demonstrate that not only different paralogous segmental duplicon but also palindromic AT-rich repeats (PATRR) are driving 22q11.2 rearrangements. In addition, we show the existence of two different inversion polymorphisms preceding rearrangement, and somatic mosaicism. The existence of different recombination sites and mechanisms in paralogues and PATRRs which are copy number expanding in the human population are a likely explanation for the high 22q11.2DS incidence.

A new blood DNA methylation signature for Koolen-de Vries syndrome: Classification of missense KANSL1 variants and comparison to fibroblast cells

Pathogenic variants in KANSL1 and 17q21.31 microdeletions are causative of Koolen-de Vries syndrome (KdVS), a neurodevelopmental syndrome with characteristic facial dysmorphia. Our previous work has shown that syndromic conditions caused by pathogenic variants in epigenetic regulatory genes have identifiable patterns of DNA methylation (DNAm) change: DNAm signatures or episignatures. Given the role of KANSL1 in histone acetylation, we tested whether variants underlying KdVS are associated with a DNAm signature. We profiled whole-blood DNAm for 13 individuals with KANSL1 variants, four individuals with 17q21.31 microdeletions, and 21 typically developing individuals, using Illumina's Infinium EPIC array. In this study, we identified a robust DNAm signature of 456 significant CpG sites in 8 individuals with KdVS, a pattern independently validated in an additional 7 individuals with KdVS. We also demonstrate the diagnostic utility of the signature and classify two KANSL1 VUS as well as four variants in individuals with atypical clinical presentation. Lastly, we investigated tissue-specific DNAm changes in fibroblast cells from individuals with KdVS. Collectively, our findings contribute to the understanding of the epigenetic landscape related to KdVS and aid in the diagnosis and classification of variants in this structurally complex genomic region.

Association of Structural Forms of 17q21.31 with the Risk of Progressive Supranuclear Palsy and MAPT Sub-haplotypes

Importance

The chromosome 17q21.31 region, containing a 900 Kb inversion that defines H1 and H2 haplotypes, represents the strongest genetic risk locus in progressive supranuclear palsy (PSP). In addition to H1 and H2, various structural forms of 17q21.31, characterized by the copy number of α, β, and γ duplications, have been identified. However, the specific effect of each structural form on the risk of PSP has never been evaluated in a large cohort study.

Objective

To assess the association of different structural forms of 17q.21.31, defined by the copy numbers of α, β, and γ duplications, with the risk of PSP and MAPT sub-haplotypes.

Design setting and participants

Utilizing whole genome sequencing data of 1,684 (1,386 autopsy confirmed) individuals with PSP and 2,392 control subjects, a case-control study was conducted to investigate the association of copy numbers of α, β, and γ duplications and structural forms of 17q21.31 with the risk of PSP. All study subjects were selected from the Alzheimer's Disease Sequencing Project (ADSP) Umbrella NG00067.v7. Data were analyzed between March 2022 and November 2023.

Main outcomes and measures

The main outcomes were the risk (odds ratios [ORs]) for PSP with 95% CIs. Risks for PSP were evaluated by logistic regression models.

Results

The copy numbers of α and β were associated with the risk of PSP only due to their correlation with H1 and H2, while the copy number of γ was independently associated with the increased risk of PSP. Each additional duplication of γ was associated with 1.10 (95% CI, 1.04-1.17; P = 0.0018) fold of increased risk of PSP when conditioning H1 and H2. For the H1 haplotype, addition γ duplications displayed a higher odds ratio for PSP: the odds ratio increases from 1.21 (95%CI 1.10-1.33, P = 5.47 × 10-5) for H1β1γ1 to 1.29 (95%CI 1.16-1.43, P = 1.35 × 10-6) for H1β1γ2, 1.45 (95%CI 1.27-1.65, P = 3.94 × 10-8) for H1β1γ3, and 1.57 (95%CI 1.10-2.26, P = 1.35 × 10-2) for H1β1γ4. Moreover, H1β1γ3 is in linkage disequilibrium with H1c (R2 = 0.31), a widely recognized MAPT sub-haplotype associated with increased risk of PSP. The proportion of MAPT sub-haplotypes associated with increased risk of PSP (i.e., H1c, H1d, H1g, H1o, and H1h) increased from 34% in H1β1γ1 to 77% in H1β1γ4.

Conclusions and relevance

This study revealed that the copy number of γ was associated with the risk of PSP independently from H1 and H2. The H1 haplotype with more γ duplications showed a higher odds ratio for PSP and were associated with MAPT sub-haplotypes with increased risk of PSP. These findings expand our understanding of how the complex structure at 17q21.31 affect the risk of PSP.

The selection landscape and genetic legacy of ancient Eurasians

The Holocene (beginning around 12,000 years ago) encompassed some of the most significant changes in human evolution, with far-reaching consequences for the dietary, physical and mental health of present-day populations. Using a dataset of more than 1,600 imputed ancient genomes1, we modelled the selection landscape during the transition from hunting and gathering, to farming and pastoralism across West Eurasia. We identify key selection signals related to metabolism, including that selection at the FADS cluster began earlier than previously reported and that selection near the LCT locus predates the emergence of the lactase persistence allele by thousands of years. We also find strong selection in the HLA region, possibly due to increased exposure to pathogens during the Bronze Age. Using ancient individuals to infer local ancestry tracts in over 400,000 samples from the UK Biobank, we identify widespread differences in the distribution of Mesolithic, Neolithic and Bronze Age ancestries across Eurasia. By calculating ancestry-specific polygenic risk scores, we show that height differences between Northern and Southern Europe are associated with differential Steppe ancestry, rather than selection, and that risk alleles for mood-related phenotypes are enriched for Neolithic farmer ancestry, whereas risk alleles for diabetes and Alzheimer's disease are enriched for Western hunter-gatherer ancestry. Our results indicate that ancient selection and migration were large contributors to the distribution of phenotypic diversity in present-day Europeans.

Structural Variation Evolution at the 15q11-q13 Disease-Associated Locus

The impact of segmental duplications on human evolution and disease is only just starting to unfold, thanks to advancements in sequencing technologies that allow for their discovery and precise genotyping. The 15q11-q13 locus is a hotspot of recurrent copy number variation associated with Prader-Willi/Angelman syndromes, developmental delay, autism, and epilepsy and is mediated by complex segmental duplications, many of which arose recently during evolution. To gain insight into the instability of this region, we characterized its architecture in human and nonhuman primates, reconstructing the evolutionary history of five different inversions that rearranged the region in different species primarily by accumulation of segmental duplications. Comparative analysis of human and nonhuman primate duplication structures suggests a human-specific gain of directly oriented duplications in the regions flanking the GOLGA cores and HERC segmental duplications, representing potential genomic drivers for the human-specific expansions. The increasing complexity of segmental duplication organization over the course of evolution underlies its association with human susceptibility to recurrent disease-associated rearrangements.

Chromosomal inversion polymorphisms shape human brain morphology

The impact of chromosomal inversions on human brain morphology remains underexplored. We studied 35 common inversions classified from genotypes of 33,018 adults with European ancestry. The inversions at 2p22.3, 16p11.2, and 17q21.31 reach genome-wide significance, followed by 8p23.1 and 6p21.33, in their association with cortical and subcortical morphology. The 17q21.31, 8p23.1, and 16p11.2 regions comprise the LRRC37, OR7E, and NPIP duplicated gene families. We find the 17q21.31 MAPT inversion region, known for harboring neurological risk, to be the most salient locus among common variants for shaping and patterning the cortex. Overall, we observe the inverted orientations decreasing brain size, with the exception that the 2p22.3 inversion is associated with increased subcortical volume and the 8p23.1 inversion is associated with increased motor cortex. These significant inversions are in the genomic hotspots of neuropsychiatric loci. Our findings are generalizable to 3,472 children and demonstrate inversions as essential genetic variation to understand human brain phenotypes.

MAPT allele and haplotype frequencies in Nigerian Africans: Population distribution and association with Parkinson's disease risk and age at onset

INTRODUCTION

The association between MAPT and PD risk may be subject to ethnic variability even within populations of similar geographical origin. Data on MAPT haplotype frequencies, and its association with PD risk in black Africans are lacking. We aimed to determine the frequencies of MAPT haplotypes and their role as risk factors for PD and age at onset in Nigerians.

METHODS

The haplotype and genotype frequencies of MAPT rs1052553 were analysed in 907 individuals with PD and 1022 age-matched healthy controls from the Nigeria Parkinson's Disease Research network cohort. Clinical data related to PD included age at study, age at onset (AAO), and disease duration.

RESULTS

The frequency of the H1 haplotype was 98.7% in PD, and 99.1% in controls (p = 0.19). The H2 haplotype was present in - 1.3% of PD and 0.9% of controls (p = 0.24). The most frequent MAPT genotype was H1H1 (PD - 97.5%, controls - 98.2%). The H1 haplotype was not associated with PD risk after accounting for gender and AAO (Odds ratio for H1/H1 vs H1/H2 and H2/H2: 0.68 (95% CI:0.39-1.28); p = 0.23).

CONCLUSIONS

Our findings support previous studies that report a low frequency of the MAPT H2 haplotype in black ancestry Africans but document its occurrence in Nigerians. The MAPT H1 haplotype was not associated with an increased risk or age at onset of PD in this cohort.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

Cognitive Performance as a Function of MAPT Haplotype: A Prospective Longitudinal Study of an Essential Tremor Cohort

Background

Cognitive impairment is a feature of essential tremor (ET). There are no studies of the genetic drivers of this association. We examined whether the microtubule-associated protein tau (MAPT) H1 haplotype is associated with cognitive performance in ET.

Methods

ET cases genotyped for the MAPT H1 and H2 haplotypes completed a battery of neuropsychological tests at baseline and four follow-up evaluations. Chi-square, t-tests, and analyses of covariance examined associations between the presence of the MAPT H1 haplotype, cognitive diagnoses of normal, mild cognitive impairment (MCI), and dementia, and performance in specific cognitive domains.

Results

We observed no evidence of cognitive differences as a function of the presence of the MAPT H1 haplotype. Specifically, cases with (n = 57) and without (n = 42) this haplotype did not differ with respect to the prevalence of diagnoses of MCI or dementia, p ≥ 0.87. Moreover, cases with vs without this haplotype did not differ in either the age or point in the disease course at which observed conversions to MCI or dementia occurred, p's ≥ 0.51. Finally, no haplotype-related differences were observed in performance in the cognitive domains of attention, executive function, language, memory, visuospatial or global ability, p's ≥ 0.21, or in changes in performance in these domains across time, p's ≥ 0.08.

Discussion

The study in an ET cohort revealed no influence of MAPT haplotypes on cognitive performance. This study serves as a valuable foundation for future studies to expand our understanding of the genetic drivers of cognitive impairment in ET.

Highlights

This study found no evidence of cognitive differences between individuals with and without the MAPT H1 haplotype. Our work provides a valuable foundation for future work to expand our knowledge of the genetic drivers of cognitive impairment in ET.

Inversion polymorphism in a complete human genome assembly

The telomere-to-telomere (T2T) complete human reference has significantly improved our ability to characterize genome structural variation. To understand its impact on inversion polymorphisms, we remapped data from 41 genomes against the T2T reference genome and compared it to the GRCh38 reference. We find a ~ 21% increase in sensitivity improving mapping of 63 inversions on the T2T reference. We identify 26 misorientations within GRCh38 and show that the T2T reference is three times more likely to represent the correct orientation of the major human allele. Analysis of 10 additional samples reveals novel rare inversions at chromosomes 15q25.2, 16p11.2, 16q22.1-23.1, and 22q11.21.

MAPT allele and haplotype frequencies in Nigerian Africans: population distribution and association with Parkinson's disease risk and age at onset

BACKGROUND

The microtubule-associated protein tau ( MAPT ) gene is critical because of its putative role in the causal pathway of neurodegenerative diseases including Parkinson's disease (PD). However, there is a lack of clarity regarding the link between the main H1 haplotype and risk of PD. Inconsistencies in reported association may be driven by genetic variability in the populations studied to date. Data on MAPT haplotype frequencies in the general population and association studies exploring the role of MAPT haplotypes in conferring PD risk in black Africans are lacking.

OBJECTIVES

To determine the frequencies of MAPT haplotypes and explore the role of the H1 haplotype as a risk factor for PD risk and age at onset in Nigerian Africans.

METHODS

The haplotype and genotype frequencies of MAPT rs1052553 were analysed using PCR-based KASP™ in 907 individuals with PD and 1,022 age-matched neurologically normal controls from the Nigeria Parkinson's Disease Research (NPDR) network cohort. Clinical data related to PD included age at study, age at onset, and disease duration.

RESULTS

The frequency of the main MAPT H1 haplotype in this cohort was 98.7% in individuals with PD, and 99.1% in healthy controls (p=0.19). The H2 haplotype was present in 41/1929 (2.1%) of the cohort (PD - 1.3%; Controls - 0.9%; p=0.24). The most frequent MAPT genotype was H1H1 (PD - 97.5%, controls - 98.2%). The H1 haplotype was not associated with PD risk after accounting for gender and age at onset (Odds ratio for H1/H1 vs H1/H2 and H2/H2: 0.68 (95% CI:0.39-1.28); p=0.23).

CONCLUSIONS

Our findings support previous studies that report a low frequency of the MAPT H2 haplotype in black ancestry Africans, but document its occurrence in the Nigerian population (2.1%). In this cohort of black Africans with PD, the MAPT H1 haplotype was not associated with an increased risk or age at onset of PD.

Inversions on human chromosomes

Human chromosome inversions are types of balanced structural variations, making them difficult to analyze. Thanks to PEM (paired-end sequencing and mapping), there has been tremendous progress in studying inversions. Inversions play an important role as an evolutionary factor, contributing to the formation of gonosomes, speciation of chimpanzees and humans, and inv17q21.3 or inv8p23.1 exhibit the features of natural selection. Both inversions have been related to pathogenic phenotype by directly affecting a gene structure (e.g., inv5p15.1q14.1), regulating gene expression (e.g., inv7q21.3q35) and by predisposing to other secondary arrangements (e.g., inv7q11.23). A polymorphism of human inversions is documented by the InvFEST database (a database that stores information about clinical predictions, validations, frequency of inversions, etc.), but only a small fraction of these inversions is validated, and a detailed analysis is complicated by the frequent location of breakpoints within regions of repetitive sequences.

Neolithic expansion and the 17q21.31 inversion in Iberia: an evolutionary approach to H2 haplotype distribution in the Near East and Europe

The chromosomal region 17q21.31 harbors a 900 kb inversion polymorphism named after the microtubule-associated protein tau (MAPT) gene. Since no recombination occurs, two haplotypes are recognized: a directly oriented variant (H1) and an inverted variant (H2). The H2 haplotype features a distribution pattern with high frequencies in the Near East and Europe, medium levels in South Asia and North Africa, and low levels elsewhere. Studies of this genomic region are relevant owing to its likely association with numerous neurodegenerative diseases. However, the causes underlying the geographic distribution of the haplotype frequencies remain a bone of contention among researchers. With this work, we have intended to outline a plausible hypothesis on the origin of the high European H2 frequencies. To that end, we have analyzed an extensive population database (including three new Iberian populations) to explore potential clinal variations of H2 frequencies. We found a sigmoidal frequency cline with an upward trend from South Asia to Europe. The maximum value was detected in the Basques from Gipuzkoa province (0.494) with the curve's inflection point in the Near East. From our results, we suggest that the most likely scenario for high H2 frequencies in Europe would be a founding event in the Near East during the late Paleolithic or early Neolithic. Subsequently, such H2 overrepresentation would have reached Europe with the arrival of the first Neolithic farmers. The current frequencies and geographic distribution of the 17q21.31 inversion suggest that the founding events mainly affected the H2D subhaplotype.

Association studies between chromosomal regions 1q21.3, 5q21.3, 14q21.2 and 17q21.31 and numbers of children in Poland

Number of children is an important human trait and studies have indicated associations with single-nucleotide polymorphisms (SNPs). Aim: to give further evidence for four associations using a large sample of Polish subjects. Data from the POPULOUS genetic database was provided from anonymous, healthy, unrelated, Polish volunteers of both sexes (N = 5760). SNPs (n = 173) studied: (a) 69 from the chromosome 17 H1/H2 inversion; (b) six from 1q21.3, 5q21.3 and 14q21.2; and (c) 98 random negative controls. Zero-inflated negative-binomial regression (z.i.) was performed (0-3 numbers of children per individual (NCI) set as non-events; adjustors: year of birth, sex). Significance level p = 0.05 with Bonferroni correction. Statistically-significant differences (with data from both sexes combined) were obtained from highly-linked inversion SNPs: representative rs12373123 gave means: homozygotes TT: 2.31 NCI (n = 1418); heterozygotes CT: 2.35 NCI (n = 554); homozygotes CC: 2.44 NCI (n = 43) (genotype p = 0.01; TTvs.CC p = 0.004; CTvs.CC p = 0.009). (Male data alone gave similar results.) Recessive modeling indicated that H2-homozygotes had 0.118 more children than H1-homozygotes + heterozygotes (z.i.-count estimates ± standard errors: CT, - 0.508 ± 0.194; TT, - 0.557 ± 0.191). The non-over-dispersed count model detected no interactions: of importance there was no significant interaction with age. No positive results were obtained from negative-control SNPs or (b). Conclusions: association between the H1/H2 inversion and numbers of children (previously reported in Iceland) has been confirmed, albeit using a different statistical model. One limitation is the small amount of data, despite initially ~ 6000 subjects. Causal studies require further investigation.

The influence of 17q21.31 and APOE genetic ancestry on neurodegenerative disease risk

Advances in genomic research over the last two decades have greatly enhanced our knowledge concerning the genetic landscape and pathophysiological processes involved in multiple neurodegenerative diseases. However, current insights arise almost exclusively from studies on individuals of European ancestry. Despite this, studies have revealed that genetic variation differentially impacts risk for, and clinical presentation of neurodegenerative disease in non-European populations, conveying the importance of ancestry in predicting disease risk and understanding the biological mechanisms contributing to neurodegeneration. We review the genetic influence of two important disease-associated loci, 17q21.31 (the "MAPT locus") and APOE, to neurodegenerative disease risk in non-European populations, touching on global population differences and evolutionary genetics by ancestry that may underlie some of these differences. We conclude there is a need to increase representation of non-European ancestry individuals in genome-wide association studies (GWAS) and biomarker analyses in order to help resolve existing disparities in understanding risk for, diagnosis of, and treatment for neurodegenerative diseases in diverse populations.

Genome-wide profiling of 5-hydroxymethylcytosines in circulating cell-free DNA reveals population-specific pathways in the development of multiple myeloma

Multiple myeloma (MM) and its precursors monoclonal gammopathy of undetermined significance (MGUS) and smoldering myeloma (SMM) are 2–3 times more common in African Americans (AA) than European Americans (EA). Although epigenetic changes are well recognized in the context of myeloma cell biology, the contribution of 5-hydroxymethylcytosines (5hmC) to racial disparities in MM is unknown. Using the 5hmC-Seal and next-generation sequencing, we profiled genome-wide 5hmC in circulating cell-free DNA (cfDNA) from 342 newly diagnosed patients with MM (n = 294), SMM (n = 18), and MGUS (n = 30). We compared differential 5hmC modifications between MM and its precursors among 227 EA and 115 AA patients. The captured 5hmC modifications in cfDNA were found to be enriched in B-cell and T-cell-derived histone modifications marking enhancers. Of the top 500 gene bodies with differential 5hmC levels between MM and SMM/MGUS, the majority (94.8%) were distinct between EA and AA and enriched with population-specific pathways, including amino acid metabolism in AA and mainly cancer-related signaling pathways in EA. These findings improved our understanding of the epigenetic contribution to racial disparities in MM and suggest epigenetic pathways that could be exploited as novel preventive strategies in high-risk populations.

Genomic architecture and functional effects of potential human inversion supergenes

Supergenes are involved in adaptation in multiple organisms, but they are little known in humans. Genomic inversions are the most common mechanism of supergene generation and maintenance. Here, we review the information about two large inversions that are the best examples of potential human supergenes. In addition, we do an integrative analysis of the newest data to understand better their functional effects and underlying genetic changes. We have found that the highly divergent haplotypes of the 17q21.31 inversion of approximately 1.5 Mb have multiple phenotypic associations, with consistent effects in brain-related traits, red and white blood cells, lung function, male and female characteristics and disease risk. By combining gene expression and nucleotide variation data, we also analysed the molecular differences between haplotypes, including gene duplications, amino acid substitutions and regulatory changes, and identify CRHR1, KANLS1 and MAPT as good candidates to be responsible for these phenotypes. The situation is more complex for the 8p23.1 inversion, where there is no clear genetic differentiation. However, the inversion is associated with several related phenotypes and gene expression differences that could be linked to haplotypes specific of one orientation. Our work, therefore, contributes to the characterization of both exceptional variants and illustrates the important role of inversions. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Detailed stratified GWAS analysis for severe COVID-19 in four European populations

Given the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), a deeper analysis of the host genetic contribution to severe COVID-19 is important to improve our understanding of underlying disease mechanisms. Here, we describe an extended genome-wide association meta-analysis of a well-characterized cohort of 3255 COVID-19 patients with respiratory failure and 12 488 population controls from Italy, Spain, Norway and Germany/Austria, including stratified analyses based on age, sex and disease severity, as well as targeted analyses of chromosome Y haplotypes, the human leukocyte antigen region and the SARS-CoV-2 peptidome. By inversion imputation, we traced a reported association at 17q21.31 to a ~0.9-Mb inversion polymorphism that creates two highly differentiated haplotypes and characterized the potential effects of the inversion in detail. Our data, together with the 5th release of summary statistics from the COVID-19 Host Genetics Initiative including non-Caucasian individuals, also identified a new locus at 19q13.33, including NAPSA, a gene which is expressed primarily in alveolar cells responsible for gas exchange in the lung.

Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1 retrotransposition; 80% of the larger inversions are balanced and affect twice as many nucleotides as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or retrotransposons. Since flanking repeats promote non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7 × 10-4 per locus per generation. Recurrent inversions exhibit a sex-chromosomal bias and co-localize with genomic disorder critical regions. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes specific haplotypes to disease-causing CNVs.

The early-life exposome modulates the effect of polymorphic inversions on DNA methylation

Polymorphic genomic inversions are chromosomal variants with intrinsic variability that play important roles in evolution, environmental adaptation, and complex traits. We investigated the DNA methylation patterns of three common human inversions, at 8p23.1, 16p11.2, and 17q21.31 in 1,009 blood samples from children from the Human Early Life Exposome (HELIX) project and in 39 prenatal heart tissue samples. We found inversion-state specific methylation patterns within and nearby flanking each inversion region in both datasets. Additionally, numerous inversion-exposure interactions on methylation levels were identified from early-life exposome data comprising 64 exposures. For instance, children homozygous at inv-8p23.1 and higher meat intake were more susceptible to TDH hypermethylation (P = 3.8 × 10⁻²²); being the inversion, exposure, and gene known risk factors for adult obesity. Inv-8p23.1 associated hypermethylation of GATA4 was also detected across numerous exposures. Our data suggests that the pleiotropic influence of inversions during development and lifetime could be substantially mediated by allele-specific methylation patterns which can be modulated by the exposome.

Analysis of the relationship between presence of common DNA sequence inversions and DNA methylation patterns suggests a role for environmental exposures (such as food intake) in mediating inversion state-specific methylation patterns.

Genetics of common cerebral small vessel disease

Cerebral small vessel disease (cSVD) is a leading cause of ischaemic and haemorrhagic stroke and a major contributor to dementia. Covert cSVD, which is detectable with brain MRI but does not manifest as clinical stroke, is highly prevalent in the general population, particularly with increasing age. Advances in technologies and collaborative work have led to substantial progress in the identification of common genetic variants that are associated with cSVD-related stroke (ischaemic and haemorrhagic) and MRI-defined covert cSVD. In this Review, we provide an overview of collaborative studies - mostly genome-wide association studies (GWAS) - that have identified >50 independent genetic loci associated with the risk of cSVD. We describe how these associations have provided novel insights into the biological mechanisms involved in cSVD, revealed patterns of shared genetic variation across cSVD traits, and shed new light on the continuum between rare, monogenic and common, multifactorial cSVD. We consider how GWAS summary statistics have been leveraged for Mendelian randomization studies to explore causal pathways in cSVD and provide genetic evidence for drug effects, and how the combination of findings from GWAS with gene expression resources and drug target databases has enabled identification of putative causal genes and provided proof-of-concept for drug repositioning potential. We also discuss opportunities for polygenic risk prediction, multi-ancestry approaches and integration with other omics data.

Genome-wide association study and functional validation implicates JADE1 in tauopathy

Primary age-related tauopathy (PART) is a neurodegenerative pathology with features distinct from but also overlapping with Alzheimer disease (AD). While both exhibit Alzheimer-type temporal lobe neurofibrillary degeneration alongside amnestic cognitive impairment, PART develops independently of amyloid-β (Aβ) plaques. The pathogenesis of PART is not known, but evidence suggests an association with genes that promote tau pathology and others that protect from Aβ toxicity. Here, we performed a genetic association study in an autopsy cohort of individuals with PART (n = 647) using Braak neurofibrillary tangle stage as a quantitative trait. We found some significant associations with candidate loci associated with AD (SLC24A4, MS4A6A, HS3ST1) and progressive supranuclear palsy (MAPT and EIF2AK3). Genome-wide association analysis revealed a novel significant association with a single nucleotide polymorphism on chromosome 4 (rs56405341) in a locus containing three genes, including JADE1 which was significantly upregulated in tangle-bearing neurons by single-soma RNA-seq. Immunohistochemical studies using antisera targeting JADE1 protein revealed localization within tau aggregates in autopsy brains with four microtubule-binding domain repeats (4R) isoforms and mixed 3R/4R, but not with 3R exclusively. Co-immunoprecipitation in post-mortem human PART brain tissue revealed a specific binding of JADE1 protein to four repeat tau lacking N-terminal inserts (0N4R). Finally, knockdown of the Drosophila JADE1 homolog rhinoceros (rno) enhanced tau-induced toxicity and apoptosis in vivo in a humanized 0N4R mutant tau knock-in model, as quantified by rough eye phenotype and terminal deoxynucleotidyl transferase dUTP nick end-labeling (TUNEL) in the fly brain. Together, these findings indicate that PART has a genetic architecture that partially overlaps with AD and other tauopathies and suggests a novel role for JADE1 as a modifier of neurofibrillary degeneration.

An approach to gene-based testing accounting for dependence of tests among nearby genes

In genome-wide association studies (GWAS), it has become commonplace to test millions of single-nucleotide polymorphisms (SNPs) for phenotypic association. Gene-based testing can improve power to detect weak signal by reducing multiple testing and pooling signal strength. While such tests account for linkage disequilibrium (LD) structure of SNP alleles within each gene, current approaches do not capture LD of SNPs falling in different nearby genes, which can induce correlation of gene-based test statistics. We introduce an algorithm to account for this correlation. When a gene's test statistic is independent of others, it is assessed separately; when test statistics for nearby genes are strongly correlated, their SNPs are agglomerated and tested as a locus. To provide insight into SNPs and genes driving association within loci, we develop an interactive visualization tool to explore localized signal. We demonstrate our approach in the context of weakly powered GWAS for autism spectrum disorder, which is contrasted to more highly powered GWAS for schizophrenia and educational attainment. To increase power for these analyses, especially those for autism, we use adaptive $P$-value thresholding, guided by high-dimensional metadata modeled with gradient boosted trees, highlighting when and how it can be most useful. Notably our workflow is based on summary statistics.

Tau and MAPT genetics in tauopathies and synucleinopathies

MAPT encodes the microtubule-associated protein tau, which is the main component of neurofibrillary tangles (NFTs) and found in other protein aggregates. These aggregates are among the pathological hallmarks of primary tauopathies such as frontotemporal dementia (FTD). Abnormal tau can also be observed in secondary tauopathies such as Alzheimer's disease (AD) and synucleinopathies such as Parkinson's disease (PD). On top of pathological findings, genetic data also links MAPT to these disorders. MAPT variations are a cause or risk factors for many tauopathies and synucleinopathies and are associated with certain clinical and pathological features in affected individuals. In addition to clinical, pathological, and genetic overlap, evidence also suggests that tau and alpha-synuclein may interact on the molecular level, and thus might collaborate in the neurodegenerative process. Understanding the role of MAPT variations in tauopathies and synucleinopathies is therefore essential to elucidate the role of tau in the pathogenesis and phenotype of those disorders, and ultimately to develop targeted therapies. In this review, we describe the role of MAPT genetic variations in tauopathies and synucleinopathies, several genotype-phenotype and pathological features, and discuss their implications for the classification and treatment of those disorders.

Evidence for opposing selective forces operating on human-specific duplicated TCAF genes in Neanderthals and humans

TRP channel-associated factor 1/2 (TCAF1/TCAF2) proteins antagonistically regulate the cold-sensor protein TRPM8 in multiple human tissues. Understanding their significance has been complicated given the locus spans a gap-ridden region with complex segmental duplications in GRCh38. Using long-read sequencing, we sequence-resolve the locus, annotate full-length TCAF models in primate genomes, and show substantial human-specific TCAF copy number variation. We identify two human super haplogroups, H4 and H5, and establish that TCAF duplications originated ~1.7 million years ago but diversified only in Homo sapiens by recurrent structural mutations. Conversely, in all archaic-hominin samples the fixation for a specific H4 haplotype without duplication is likely due to positive selection. Here, our results of TCAF copy number expansion, selection signals in hominins, and differential TCAF2 expression between haplogroups and high TCAF2 and TRPM8 expression in liver and prostate in modern-day humans imply TCAF diversification among hominins potentially in response to cold or dietary adaptations.

Duplications of gene segments can allow novel physiological adaptations to evolve. A detailed analysis of the TCAF gene family in primates and archaic humans suggest rapid duplication and diversification in this gene family is associated with cold or dietary adaptations.

Insights into genetic factors contributing to variability in SARS-CoV-2 susceptibility and COVID-19 disease severity

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we applied colocalization to compare summary statistics for 16 GWASs from the COVID-19 Host Genetics Initiative to investigate similarities and differences in their genetic signals. We identified 9 loci associated with susceptibility (one with two independent GWAS signals; one with an ethnicity-specific signal), 14 associated with severity (one with two independent GWAS signals; two with ethnicity-specific signals) and one harboring two discrepant GWAS signals (one for susceptibility; one for severity). Utilizing colocalization we also identified 45 GTEx tissues that had eQTL(s) for 18 genes strongly associated with GWAS signals in eleven loci (1-4 genes per locus). Some of these genes showed tissue-specific altered expression and others showed altered expression in up to 41 different tissue types. Our study provides insights into the complex molecular mechanisms underlying inherited predispositions to COVID-19-disease phenotypes.

Diverse Molecular Mechanisms Contribute to Differential Expression of Human Duplicated Genes

Emerging evidence links genes within human-specific segmental duplications (HSDs) to traits and diseases unique to our species. Strikingly, despite being nearly identical by sequence (>98.5%), paralogous HSD genes are differentially expressed across human cell and tissue types, though the underlying mechanisms have not been examined. We compared cross-tissue mRNA levels of 75 HSD genes from 30 families between humans and chimpanzees and found expression patterns consistent with relaxed selection on or neofunctionalization of derived paralogs. In general, ancestral paralogs exhibited greatest expression conservation with chimpanzee orthologs, though exceptions suggest certain derived paralogs may retain or supplant ancestral functions. Concordantly, analysis of long-read isoform sequencing data sets from diverse human tissues and cell lines found that about half of derived paralogs exhibited globally lower expression. To understand mechanisms underlying these differences, we leveraged data from human lymphoblastoid cell lines (LCLs) and found no relationship between paralogous expression divergence and post-transcriptional regulation, sequence divergence, or copy-number variation. Considering cis-regulation, we reanalyzed ENCODE data and recovered hundreds of previously unidentified candidate CREs in HSDs. We also generated large-insert ChIP-sequencing data for active chromatin features in an LCL to better distinguish paralogous regions. Some duplicated CREs were sufficient to drive differential reporter activity, suggesting they may contribute to divergent cis-regulation of paralogous genes. This work provides evidence that cis-regulatory divergence contributes to novel expression patterns of recent gene duplicates in humans.

Identifying chromosomal subpopulations based on their recombination histories advances the study of the genetic basis of phenotypic traits

Recombination is a main source of genetic variability. However, the potential role of the variation generated by recombination in phenotypic traits, including diseases, remains unexplored because there is currently no method to infer chromosomal subpopulations based on recombination pattern differences. We developed recombClust, a method that uses SNP-phased data to detect differences in historic recombination in a chromosome population. We validated our method by performing simulations and by using real data to accurately predict the alleles of well-known recombination modifiers, including common inversions in Drosophila melanogaster and human, and the chromosomes under selective pressure at the lactase locus in humans. We then applied recombClust to the complex human 1q21.1 region, where nonallelic homologous recombination produces deleterious phenotypes. We discovered and validated the presence of two different recombination histories in these regions that significantly associated with the differential expression of ANKRD35 in whole blood and that were in high linkage with variants previously associated with hypertension. By detecting differences in historic recombination, our method opens a way to assess the influence of recombination variation in phenotypic traits.

Clinical Genetics Can Solve the Pitfalls of Genome-Wide Investigations: Lesson from Mismapping a Loss-of-Function Variant in KANSL1

Massive parallel sequencing of 70 genes in a girl with a suspicion of chromatinopathy detected the (NM_015443.4:)c.985_986delTT variant in exon 2 of KANSL1, which led to a diagnostic consideration of Koolen De Vries syndrome. The same variant was present in the healthy mother, consistent with either incomplete penetrance or variant mismapping. A network of second opinion was implemented among clinical geneticists first, and a diagnosis of Koolen De Vries syndrome was considered unlikely. By MLPA, a duplication spanning exons 1-3 of KANSL1 was detected in both the mother and the daughter. On cDNA sequencing, biallelic wild type mRNA was observed. We concluded that the variant affects the noncoding duplicated gene region in our family, and we finally classified it as benign. Parallel wide genomic sequencing is increasingly the first genetic investigation in individuals with intellectual disability. The c.985_986delTT variant in KANSL1 was described both in individuals with typical KdVS and in a limited number of healthy subjects. This report highlights the role of clinical genetics to correctly classify variants and to define proper clinical and diagnostic correlations.

Determining the impact of uncharacterized inversions in the human genome by droplet digital PCR

Despite the interest in characterizing genomic variation, the presence of large repeats at the breakpoints hinders the analysis of many structural variants. This is especially problematic for inversions, since there is typically no gain or loss of DNA. Here, we tested novel linkage-based droplet digital PCR (ddPCR) assays to study 20 inversions ranging from 3.1 to 742 kb flanked by inverted repeats (IRs) up to 134 kb long. Of those, we validated 13 inversions predicted by different genome-wide techniques. In addition, we obtained new experimental human population information across 95 African, European, and East Asian individuals for 16 inversions, including four already validated variants without high-throughput genotyping methods. Through comparison with previous data, independent replicates and both inversion breakpoints, we demonstrate that the technique is highly accurate and reproducible. Most studied inversions are widespread across continents, and their frequency is negatively correlated with genetic length. Moreover, all except two show clear signs of being recurrent, and we could better define the factors affecting recurrence levels and estimate the inversion rate across the genome. Finally, the generated genotypes have allowed us to check inversion functional effects, validating gene expression differences reported before for two inversions and finding new candidate associations. Therefore, the developed methodology makes it possible to screen these and other complex genomic variants quickly in a large number of samples for the first time, highlighting the importance of direct genotyping to assess their potential consequences and clinical implications.

Adaptive archaic introgression of copy number variants and the discovery of previously unknown human genes

Copy number variants (CNVs) are subject to stronger selective pressure than single-nucleotide variants, but their roles in archaic introgression and adaptation have not been systematically investigated. We show that stratified CNVs are significantly associated with signatures of positive selection in Melanesians and provide evidence for adaptive introgression of large CNVs at chromosomes 16p11.2 and 8p21.3 from Denisovans and Neanderthals, respectively. Using long-read sequence data, we reconstruct the structure and complex evolutionary history of these polymorphisms and show that both encode positively selected genes absent from most human populations. Our results collectively suggest that large CNVs originating in archaic hominins and introgressed into modern humans have played an important role in local population adaptation and represent an insufficiently studied source of large-scale genetic variation.

The 22q11 low copy repeats are characterized by unprecedented size and structural variability

Low copy repeats (LCRs) are recognized as a significant source of genomic instability, driving genome variability and evolution. The Chromosome 22 LCRs (LCR22s) mediate nonallelic homologous recombination (NAHR) leading to the 22q11 deletion syndrome (22q11DS). However, LCR22s are among the most complex regions in the genome, and their structure remains unresolved. The difficulty in generating accurate maps of LCR22s has also hindered localization of the deletion end points in 22q11DS patients. Using fiber FISH and Bionano optical mapping, we assembled LCR22 alleles in 187 cell lines. Our analysis uncovered an unprecedented level of variation in LCR22s, including LCR22A alleles ranging in size from 250 to 2000 kb. Further, the incidence of various LCR22 alleles varied within different populations. Additionally, the analysis of LCR22s in 22q11DS patients and their parents enabled further refinement of the rearrangement site within LCR22A and -D, which flank the 22q11 deletion. The NAHR site was localized to a 160-kb paralog shared between the LCR22A and -D in seven 22q11DS patients. Thus, we present the most comprehensive map of LCR22 variation to date. This will greatly facilitate the investigation of the role of LCR variation as a driver of 22q11 rearrangements and the phenotypic variability among 22q11DS patients.

scoreInvHap: Inversion genotyping for genome-wide association studies

Polymorphic inversions contribute to adaptation and phenotypic variation. However, large multi-centric association studies of inversions remain challenging. We present scoreInvHap, a method to genotype inversions from SNP data for genome-wide association studies (GWASs), overcoming important limitations of current methods and outperforming them in accuracy and applicability. scoreInvHap calls individual inversion-genotypes from a similarity score to the SNPs of experimentally validated references. It can be used on different sources of SNP data, including those with low SNP coverage such as exome sequencing, and is easily adaptable to genotype new inversions, either in humans or in other species. We present 20 human inversions that can be reliably and easily genotyped with scoreInvHap to discover their role in complex human traits, and illustrate a first genome-wide association study of experimentally-validated human inversions. scoreInvHap is implemented in R and it is freely available from Bioconductor.

Chromosomal inversions are structural variants consisting on an orientation change of a chromosome segment. Inversions have been linked to some phenotypic differences between individuals and to genetic divergence. However, their overall contribution to complex diseases is largely underdetermined as there are no high-throughput methods to call inversion-genotypes in large cohort studies. Here, we propose a new method, scoreInvHap, to call individual inversion genotypes from their haplotype similarity. We show that scoreInvHap has a high performance when analyzing heterogeneous sources of SNP data. Our current implementation contains 20 human inversions that can be readily genotyped in existing GWAS datasets. We exemplify the utility of scoreInvHap by running the first-genome wide association of experimentally validated inversions and a multi-centric inversion association study. All in all, scoreInvHap can substantially contribute to increase our knowledge of the role of chromosomal inversions in complex diseases by re-analyzing data from existing genetic association studies.

Bayesian meta-analysis across genome-wide association studies of diverse phenotypes

Genome-wide association studies (GWAS) are a powerful tool for understanding the genetic basis of diseases and traits, but most studies have been conducted in isolation, with a focus on either a single or a set of closely related phenotypes. We describe MetABF, a simple Bayesian framework for performing integrative meta-analysis across multiple GWAS using summary statistics. The approach is applicable across a wide range of study designs and can increase the power by 50% compared with standard frequentist tests when only a subset of studies have a true effect. We demonstrate its utility in a meta-analysis of 20 diverse GWAS which were part of the Wellcome Trust Case Control Consortium 2. The novelty of the approach is its ability to explore, and assess the evidence for a range of possible true patterns of association across studies in a computationally efficient framework.

Genomic inversions and GOLGA core duplicons underlie disease instability at the 15q25 locus

Human chromosome 15q25 is involved in several disease-associated structural rearrangements, including microdeletions and chromosomal markers with inverted duplications. Using comparative fluorescence in situ hybridization, strand-sequencing, single-molecule, real-time sequencing and Bionano optical mapping analyses, we investigated the organization of the 15q25 region in human and nonhuman primates. We found that two independent inversions occurred in this region after the fission event that gave rise to phylogenetic chromosomes XIV and XV in humans and great apes. One of these inversions is still polymorphic in the human population today and may confer differential susceptibility to 15q25 microdeletions and inverted duplications. The inversion breakpoints map within segmental duplications containing core duplicons of the GOLGA gene family and correspond to the site of an ancestral centromere, which became inactivated about 25 million years ago. The inactivation of this centromere likely released segmental duplications from recombination repression typical of centromeric regions. We hypothesize that this increased the frequency of ectopic recombination creating a hotspot of hominid inversions where dispersed GOLGA core elements now predispose this region to recurrent genomic rearrangements associated with disease.

The intramembrane protease SPPL2c promotes male germ cell development by cleaving phospholamban

Signal peptide peptidase (SPP) and the four homologous SPP-like (SPPL) proteases constitute a family of intramembrane aspartyl proteases with selectivity for type II-oriented transmembrane segments. Here, we analyse the physiological function of the orphan protease SPPL2c, previously considered to represent a non-expressed pseudogene. We demonstrate proteolytic activity of SPPL2c towards selected tail-anchored proteins. Despite shared ER localisation, SPPL2c and SPP exhibit distinct, though partially overlapping substrate spectra and inhibitory profiles, and are organised in different high molecular weight complexes. Interestingly, SPPL2c is specifically expressed in murine and human testis where it is primarily localised in spermatids. In mice, SPPL2c deficiency leads to a partial loss of elongated spermatids and reduced motility of mature spermatozoa, but preserved fertility. However, matings of male and female SPPL2c-/- mice exhibit reduced litter sizes. Using proteomics we identify the sarco/endoplasmic reticulum Ca2+-ATPase (SERCA2)-regulating protein phospholamban (PLN) as a physiological SPPL2c substrate. Accumulation of PLN correlates with a decrease in intracellular Ca2+ levels in elongated spermatids that likely contribute to the compromised male germ cell differentiation and function of SPPL2c-/- mice.

Lung cancer health disparities

Compared with all other racial and ethnic groups in the United States, African Americans are disproportionally affected by lung cancer, both in terms of incidence and survival. It is likely that smoking, as the main etiological factor associated with lung cancer, contributes to these disparities, but the precise mechanism is still unclear. This paper seeks to explore the history of lung cancer disparities and review to the literature regarding the various factors that contribute to them.

Long-read sequence and assembly of segmental duplications

We have developed a computational method based on polyploid phasing of long sequence reads to resolve collapsed regions of segmental duplications within genome assemblies. Segmental Duplication Assembler (SDA; https://github.com/mvollger/SDA ) constructs graphs in which paralogous sequence variants define the nodes and long-read sequences provide attraction and repulsion edges, enabling the partition and assembly of long reads corresponding to distinct paralogs. We apply it to single-molecule, real-time sequence data from three human genomes and recover 33-79 megabase pairs (Mb) of duplications in which approximately half of the loci are diverged (<99.8%) compared to the reference genome. We show that the corresponding sequence is highly accurate (>99.9%) and that the diverged sequence corresponds to copy-number-variable paralogs that are absent from the human reference genome. Our method can be applied to other complex genomes to resolve the last gene-rich gaps, improve duplicate gene annotation, and better understand copy-number-variant genetic diversity at the base-pair level.

The role of monogenic genes in idiopathic Parkinson's disease

In the past two decades, mutations in multiple genes have been linked to autosomal dominant or recessive forms of monogenic Parkinson's disease (PD). Collectively, these monogenic (often familial) cases account for less than 5% of all PD, the majority being apparently sporadic cases. More recently, large-scale genome-wide association studies have identified over 40 loci that increase risk of PD. Importantly, there is overlap between monogenic and sporadic PD genes, particularly for the loci that contain the genes SNCA and LRRK2, which are mutated in monogenic dominant PD. There have also been reports of idiopathic PD cases with heterozygous variants in autosomal recessive genes suggesting that these mutations may increase risk of PD. These observations suggest that monogenic and idiopathic PD may have shared pathogenic mechanisms. Here, we focus mainly on the role of monogenic PD genes that represent pleomorphic risk loci for idiopathic PD. We also discuss the functional mechanisms that may play a role in increasing risk of disease in both monogenic and idiopathic forms.

Divergent brain gene expression patterns associate with distinct cell-specific tau neuropathology traits in progressive supranuclear palsy

Progressive supranuclear palsy (PSP) is a neurodegenerative parkinsonian disorder characterized by tau pathology in neurons and glial cells. Transcriptional regulation has been implicated as a potential mechanism in conferring disease risk and neuropathology for some PSP genetic risk variants. However, the role of transcriptional changes as potential drivers of distinct cell-specific tau lesions has not been explored. In this study, we integrated brain gene expression measurements, quantitative neuropathology traits and genome-wide genotypes from 268 autopsy-confirmed PSP patients to identify transcriptional associations with unique cell-specific tau pathologies. We provide individual transcript and transcriptional network associations for quantitative oligodendroglial (coiled bodies = CB), neuronal (neurofibrillary tangles = NFT), astrocytic (tufted astrocytes = TA) tau pathology, and tau threads and genomic annotations of these findings. We identified divergent patterns of transcriptional associations for the distinct tau lesions, with the neuronal and astrocytic neuropathologies being the most different. We determined that NFT are positively associated with a brain co-expression network enriched for synaptic and PSP candidate risk genes, whereas TA are positively associated with a microglial gene-enriched immune network. In contrast, TA is negatively associated with synaptic and NFT with immune system transcripts. Our findings have implications for the diverse molecular mechanisms that underlie cell-specific vulnerability and disease risk in PSP.

Comparative Transcriptome Profiling Reveals Coding and Noncoding RNA Differences in NSCLC from African Americans and European Americans

Purpose: To determine whether racial differences in gene and miRNA expression translates to differences in lung tumor biology with clinical relevance in African Americans (AAs) and European Americans (EAs).Experimental Design: The NCI-Maryland Case Control Study includes seven Baltimore City hospitals and is overrepresented with AA patients (∼40%). Patients that underwent curative NSCLC surgery between 1998 and 2014 were enrolled. Comparative molecular profiling used mRNA (n = 22 AAs and 19 EAs) and miRNA (n = 42 AAs and 55 EAs) expression arrays to track differences in paired fresh frozen normal tissues and lung tumor specimens from AAs and EAs. Pathway enrichment, predicted drug response, tumor microenvironment infiltration, cancer immunotherapy antigen profiling, and miRNA target enrichment were assessed.Results: AA-enriched differential gene expression was characterized by stem cell and invasion pathways. Differential gene expression in lung tumors from EAs was primarily characterized by cell proliferation pathways. Population-specific gene expression was partly driven by population-specific miRNA expression profiles. Drug susceptibility predictions revealed a strong inverse correlation between AA resistance and EA sensitivity to the same panel of drugs. Statistically significant differences in M1 and M2 macrophage infiltration were observed in AAs (P < 0.05); however, PD-L1, PD-L2 expression was similar between both.Conclusions: Comparative transcriptomic profiling revealed clear differences in lung tumor biology between AAs and EAs. Increased participation by AAs in lung cancer clinical trials are needed to integrate, and leverage, transcriptomic differences with other clinical information to maximize therapeutic benefit in both AAs and EAs. Clin Cancer Res; 23(23); 7412-25. ©2017 AACR.

Inversion variants in human and primate genomes

For many years, inversions have been proposed to be a direct driving force in speciation since they suppress recombination when heterozygous. Inversions are the most common large-scale differences among humans and great apes. Nevertheless, they represent large events easily distinguishable by classical cytogenetics, whose resolution, however, is limited. Here, we performed a genome-wide comparison between human, great ape, and macaque genomes using the net alignments for the most recent releases of genome assemblies. We identified a total of 156 putative inversions, between 103 kb and 91 Mb, corresponding to 136 human loci. Combining literature, sequence, and experimental analyses, we analyzed 109 of these loci and found 67 regions inverted in one or multiple primates, including 28 newly identified inversions. These events overlap with 81 human genes at their breakpoints, and seven correspond to sites of recurrent rearrangements associated with human disease. This work doubles the number of validated primate inversions larger than 100 kb, beyond what was previously documented. We identified 74 sites of errors, where the sequence has been assembled in the wrong orientation, in the reference genomes analyzed. Our data serve two purposes: First, we generated a map of evolutionary inversions in these genomes representing a resource for interrogating differences among these species at a functional level; second, we provide a list of misassembled regions in these primate genomes, involving over 300 Mb of DNA and 1978 human genes. Accurately annotating these regions in the genome references has immediate applications for evolutionary and biomedical studies on primates.

Tau haplotypes support the Asian ancestry of the Roma population settled in the Basque Country

We examined tau haplotype frequencies in two different ethnical groups from the Basque Country (BC): Roma people and residents of European ancestry (general population). In addition, we analyzed the spatial distribution of tau haplotypes in Eurasian populations to explore the genetic affinities of the Romani groups living in Europe in a broader scope. The 17q21.31 genomic region was characterized through the genotyping of two diagnostic single nucleotide polymorphisms, SNPs (rs10514879 and rs199451), which allow the identification of H1 and H2 haplotypes. A significant heterozygous deficit was detected in the Romani for rs10514879. The H2 haplotype frequency proved to be more than twice in the BC general population (0.283) than in the Roma people (0.127). In contrast, H2 frequency proved to be very similar between Basque and Hungarian Romani, and similar to the H2 frequencies found in northwestern India and Pakistan as well. Several statistical analyses unveiled genetic structuring for the MAPT diversity, mirrored in a significant association between geography and genetic distances, with an upward trend of H2 haplotype frequencies from Asia to Europe. Yet, Roma samples did not fit into this general spatial patterning because of their discrepancy between geographical position and H2 frequency. Despite the long spatial coexistence in the Basque region between the residents of European ancestry and the Roma, the latter have preserved their Asian genetic ancestry. Bearing in mind the lack of geographical barriers between both ethnical groups, these findings support the notion that sociocultural mores might promote assortative matings in human populations.

Detailed analysis of inversions predicted between two human genomes: errors, real polymorphisms, and their origin and population distribution

The growing catalogue of structural variants in humans often overlooks inversions as one of the most difficult types of variation to study, even though they affect phenotypic traits in diverse organisms. Here, we have analysed in detail 90 inversions predicted from the comparison of two independently assembled human genomes: the reference genome (NCBI36/HG18) and HuRef. Surprisingly, we found that two thirds of these predictions (62) represent errors either in assembly comparison or in one of the assemblies, including 27 misassembled regions in HG18. Next, we validated 22 of the remaining 28 potential polymorphic inversions using different PCR techniques and characterized their breakpoints and ancestral state. In addition, we determined experimentally the derived allele frequency in Europeans for 17 inversions (DAF = 0.01-0.80), as well as the distribution in 14 worldwide populations for 12 of them based on the 1000 Genomes Project data. Among the validated inversions, nine have inverted repeats (IRs) at their breakpoints, and two show nucleotide variation patterns consistent with a recurrent origin. Conversely, inversions without IRs have a unique origin and almost all of them show deletions or insertions at the breakpoints in the derived allele mediated by microhomology sequences, which highlights the importance of mechanisms like FoSTeS/MMBIR in the generation of complex rearrangements in the human genome. Finally, we found several inversions located within genes and at least one candidate to be positively selected in Africa. Thus, our study emphasizes the importance of careful analysis and validation of large-scale genomic predictions to extract reliable biological conclusions.