The impact of a fine-scale population stratification on rare variant association test results

Immune mechanisms in multiple sclerosis: CD3 levels on CD28+ CD4+ T cells link antibody responses to human herpesvirus 6

Compelling evidence suggests a significant association between antibody-mediated immune responses and multiple sclerosis (MS). However, the exact causal relationships between these immune responses and MS remain unclear. In this study, we conducted a comprehensive examination of the link between antibody-mediated immune responses and MS via Mendelian randomization (MR) analysis to identify specific infectious pathogens potentially involved in the onset and progression of MS. We compared immune cell infiltration between MS patients and control subjects. Furthermore, single-cell sequencing was employed to conduct a comparative analysis of the marker genes associated with each cell subtype between individuals diagnosed with MS and the control cohort. We revealed connections between antibody-mediated immune responses and immune cells, as well as the associations between these immune cells and MS. We discovered that CD3 levels on CD28+ CD4+ T cells significantly influence MS progression by altering the ratio of human herpesvirus 6 (HHV-6). These findings provide novel insights into the biological mechanisms underlying HHV-6-mediated MS.

GRM1 as a Candidate Gene for Buffalo Fertility: Insights from Genome-Wide Association Studies and Its Role in the FOXO Signaling Pathway

BACKGROUND

Water buffaloes represent a crucial genetic resource for the global dairy industry, yet enhancements in their production performance remain relatively constrained. The advent of advanced sequencing technologies, coupled with genome-wide association studies (GWASs), has significantly boosted the potential for breeding superior-quality water buffalo.

METHODS

An integrated genomic analysis was performed on sequencing data from 100 water buffaloes, utilizing the high-quality UOA_WB_1 genome assembly as a reference. This study particularly emphasized reproduction-related traits, with a focus on age at first calving (AFC).

RESULTS

Our analysis revealed two significant single-nucleotide polymorphisms (SNPs). Based on these genetic markers, the GRM1 gene was identified as a candidate gene. This gene shows substantial involvement in various reproduction-associated pathways, including the FOXO signaling pathway, calcium signaling pathway, and estrogen signaling pathway.

CONCLUSIONS

The identification of GRM1 as a candidate gene provides a robust theoretical basis for molecular breeding strategies aimed at enhancing fertility in water buffaloes. These findings offer critical scientific support for optimizing breeding programs, thereby improving overall production efficiency.

Genome-Wide Association Studies for Lactation Performance in Buffaloes

Background: Buffaloes are considered an indispensable genetic resource for dairy production. However, improvements in lactation performance have been relatively limited. Advances in sequencing technology, combined with genome-wide association studies, have facilitated the breeding of high-quality buffalo. Methods: We conducted an integrated analysis of genomic sequencing data from 120 water buffalo, the high-quality water buffalo genome assembly designated as UOA_WB_1, and milk production traits, including 305-day milk yield (MY), peak milk yield (PM), total protein yield (PY), protein percentage (PP), fat percentage (FP), and total milk fat yield (FY). Results: The results identified 56 significant SNPs, and based on these markers, 54 candidate genes were selected. These candidate genes were significantly enriched in lactation-related pathways, such as the cAMP signaling pathway (ABCC4), TGF-β signaling pathway (LEFTY2), Wnt signaling pathway (CAMK2D), and metabolic pathways (DGAT1). Conclusions: These candidate genes (e.g., ABCC4, LEFTY2, CAMK2D, DGAT1) provide a substantial theoretical foundation for molecular breeding to enhance milk production in buffaloes.

Appraisal of the causal effect of Chlamydia trachomatis infection on epithelial ovarian cancer risk: a two-sample Mendelian randomisation study

Background

History of Chlamydia trachomatis infection has previously been associated with epithelial ovarian cancer (EOC) in observational studies. We conducted a two-sample univariable Mendelian randomisation (MR) study to examine whether genetically predicted seropositivity to the C. trachomatis major outer membrane protein (momp) D is causally associated with EOC.

Methods

MR analyses employed genetic associations derived from UK Biobank as proxies for momp D seropositivity in 25 509 EOC cases and 40 941 controls that participated in the Ovarian Cancer Association Consortium. Findings were replicated using a GWAS meta-analyses of global biobanks including the UK Biobank, FinnGen and BioBank Japan.

Results

Genetically predicted momp D seropositivity was associated with overall and high-grade serous EOC risk in inverse-variance weighted (IVW) and MR-Egger univariable MR analysis (odds ratio (OR) 1.06; 95% confidence interval (CI) 1.02-1.10, and OR 1.08; 95%CI 1.01-1.16, respectively). Replication yielded similar results for overall EOC (OR 1.11; 95%CI 1.01-1.22).

Conclusion

This MR study supports a causative link between C. trachomatis infection and overall and high-grade serous EOC.

The importance of family-based sampling for biobanks

Biobanks aim to improve our understanding of health and disease by collecting and analysing diverse biological and phenotypic information in large samples. So far, biobanks have largely pursued a population-based sampling strategy, where the individual is the unit of sampling, and familial relatedness occurs sporadically and by chance. This strategy has been remarkably efficient and successful, leading to thousands of scientific discoveries across multiple research domains, and plans for the next wave of biobanks are underway. In this Perspective, we discuss the strengths and limitations of a complementary sampling strategy for future biobanks based on oversampling of close genetic relatives. Such family-based samples facilitate research that clarifies causal relationships between putative risk factors and outcomes, particularly in estimates of genetic effects, because they enable analyses that reduce or eliminate confounding due to familial and demographic factors. Family-based biobank samples would also shed new light on fundamental questions across multiple fields that are often difficult to explore in population-based samples. Despite the potential for higher costs and greater analytical complexity, the many advantages of family-based samples should often outweigh their potential challenges.

Integrative Analysis of Germline Rare Variants in Clear and Non-clear Cell Renal Cell Carcinoma

Background and objective

Previous germline studies on renal cell carcinoma (RCC) have usually pooled clear and non-clear cell RCCs and have not adequately accounted for population stratification, which might have led to an inaccurate estimation of genetic risk. Here, we aim to analyze the major germline drivers of RCC risk and clinically relevant but underexplored germline variant types.

Methods

We first characterized germline pathogenic variants (PVs), cryptic splice variants, and copy number variants (CNVs) in 1436 unselected RCC patients. To evaluate the enrichment of PVs in RCC, we conducted a case-control study of 1356 RCC patients ancestry matched with 16 512 cancer-free controls using approaches accounting for population stratification and histological subtypes, followed by characterization of secondary somatic events.

Key findings and limitations

Clear cell RCC patients (n = 976) exhibited a significant burden of PVs in VHL compared with controls (odds ratio [OR]: 39.1, p = 4.95e-05). Non-clear cell RCC patients (n = 380) carried enrichment of PVs in FH (OR: 77.9, p = 1.55e-08) and MET (OR: 1.98e11, p = 2.07e-05). In a CHEK2-focused analysis with European participants, clear cell RCC (n = 906) harbored nominal enrichment of low-penetrance CHEK2 variants-p.Ile157Thr (OR: 1.84, p = 0.049) and p.Ser428Phe (OR: 5.20, p = 0.045), while non-clear cell RCC (n = 295) exhibited nominal enrichment of CHEK2 loss of function PVs (OR: 3.51, p = 0.033). Patients with germline PVs in FH, MET, and VHL exhibited significantly earlier age of cancer onset than patients without germline PVs (mean: 46.0 vs 60.2 yr, p < 0.0001), and more than half had secondary somatic events affecting the same gene (n = 10/15, 66.7%). Conversely, CHEK2 PV carriers exhibited a similar age of onset to patients without germline PVs (mean: 60.1 vs 60.2 yr, p = 0.99), and only 30.4% carried somatic events in CHEK2 (n = 7/23). Finally, pathogenic germline cryptic splice variants were identified in SDHA and TSC1, and pathogenic germline CNVs were found in 18 patients, including CNVs in FH, SDHA, and VHL.

Conclusions and clinical implications

This analysis supports the existing link between several RCC risk genes and RCC risk manifesting in earlier age of onset. It calls for caution when assessing the role of CHEK2 due to the burden of founder variants with varying population frequency. It also broadens the definition of the RCC germline landscape of pathogenicity to incorporate previously understudied types of germline variants.

Patient summary

In this study, we carefully compared the frequency of rare inherited mutations with a focus on patients' genetic ancestry. We discovered that subtle variations in genetic background may confound a case-control analysis, especially in evaluating the cancer risk associated with specific genes, such as CHEK2. We also identified previously less explored forms of rare inherited mutations, which could potentially increase the risk of kidney cancer.

Risk and protection of different rare protein-coding variants of complement component C4A in age-related macular degeneration

Introduction: Age-related macular degeneration (AMD) is the leading cause of central vision loss in the elderly. One-third of the genetic contribution to this disease remains unexplained. Methods: We analyzed targeted sequencing data from two independent cohorts (4,245 cases, 1,668 controls) which included genomic regions of known AMD loci in 49 genes. Results: At a false discovery rate of <0.01, we identified 11 low-frequency AMD variants (minor allele frequency <0.05). Two of those variants were present in the complement C4A gene, including the replacement of the residues that contribute to the Rodgers-1/Chido-1 blood group antigens: [VDLL1207-1210ADLR (V1207A)] with discovery odds ratio (OR) = 1.7 (p = 3.2 × 10-5) which was replicated in the UK Biobank dataset (3,294 cases, 200,086 controls, OR = 1.52, p = 0.037). A novel variant associated with reduced risk for AMD in our discovery cohort was P1120T, one of the four C4A-isotypic residues. Gene-based tests yielded aggregate effects of nonsynonymous variants in 10 genes including C4A, which were associated with increased risk of AMD. In human eye tissues, immunostaining demonstrated C4A protein accumulation in and around endothelial cells of retinal and choroidal vasculature, and total C4 in soft drusen. Conclusion: Our results indicate that C4A protein in the complement activation pathways may play a role in the pathogenesis of AMD.

Severe adult hemophagocytic lymphohistiocytosis (HLHa) correlates with HLH-related gene variants

BACKGROUND

The contribution of genetic factors to the severity of adult hemophagocytic lymphohistiocytosis (HLHa) remains unclear.

OBJECTIVE

We sought to assess a potential link between HLHa outcomes and HLH-related gene variants.

METHODS

Clinical characteristics of 130 HLHa patients (age ≥ 18 years and HScore ≥ 169) and genotype of 8 HLH-related genes (LYST, PRF1, UNC13-D, STX11, STXBP2, RAB27A, XIAP, and SAP) were collected. A total of 34 variants found in only 6 genes were selected on the basis of their frequency and criteria predicted to impair protein function. Severity was defined by refractory disease to HLH treatment, death, or transfer to an intensive care unit.

RESULTS

HLHa-associated diseases (ADs) were neoplasia (n = 49 [37.7%]), autoimmune/inflammatory disease (n = 33 [25.4%]), or idiopathic when no AD was identified (n = 48 [36.9%]). Infectious events occurred in 76 (58.5%) patients and were equally distributed in all ADs. Severe and refractory HLHa were observed in 80 (61.5%) and 64 (49.2%) patients, respectively. HScore, age, sex ratio, AD, and infectious events showed no significant association with HLHa severity. Variants were identified in 71 alleles and were present in 56 (43.1%) patients. They were distributed as follows: 44 (34.4%), 9 (6.9%), and 3 (2.3%) patients carrying 1, 2, and 3 variant alleles, respectively. In a logistic regression model, only the number of variants was significantly associated with HLHa severity (1 vs 0: 3.86 [1.73-9.14], P = .0008; 2-3 vs 0: 29.4 [3.62-3810], P = .0002) and refractoriness (1 vs 0: 2.47 [1.17-5.34], P = .018; 2-3 vs 0: 13.2 [2.91-126.8], P = .0003).

CONCLUSIONS

HLH-related gene variants may be key components to the severity and refractoriness of HLHa.

Excalibur: A new ensemble method based on an optimal combination of aggregation tests for rare-variant association testing for sequencing data

The development of high-throughput next-generation sequencing technologies and large-scale genetic association studies produced numerous advances in the biostatistics field. Various aggregation tests, i.e. statistical methods that analyze associations of a trait with multiple markers within a genomic region, have produced a variety of novel discoveries. Notwithstanding their usefulness, there is no single test that fits all needs, each suffering from specific drawbacks. Selecting the right aggregation test, while considering an unknown underlying genetic model of the disease, remains an important challenge. Here we propose a new ensemble method, called Excalibur, based on an optimal combination of 36 aggregation tests created after an in-depth study of the limitations of each test and their impact on the quality of result. Our findings demonstrate the ability of our method to control type I error and illustrate that it offers the best average power across all scenarios. The proposed method allows for novel advances in Whole Exome/Genome sequencing association studies, able to handle a wide range of association models, providing researchers with an optimal aggregation analysis for the genetic regions of interest.

Integrative Analysis of Germline Rare Variants in Clear and Non-Clear Cell Renal Cell Carcinoma

IMPORTANCE

RCC encompasses a set of histologically distinct cancers with a high estimated genetic heritability, of which only a portion is currently explained. Previous rare germline variant studies in RCC have usually pooled clear and non-clear cell RCCs and have not adequately accounted for population stratification that may significantly impact the interpretation and discovery of certain candidate risk genes.

OBJECTIVE

To evaluate the enrichment of germline PVs in established cancer-predisposing genes (CPGs) in clear cell and non-clear cell RCC patients compared to cancer-free controls using approaches that account for population stratification and to identify unconventional types of germline RCC risk variants that confer an increased risk of developing RCC.

DESIGN SETTING AND PARTICIPANTS

In 1,436 unselected RCC patients with sufficient data quality, we systematically identified rare germline PVs, cryptic splice variants, and copy number variants (CNVs). From this unselected cohort, 1,356 patients were ancestry-matched with 16,512 cancer-free controls, and gene-level enrichment of rare germline PVs were assessed in 143 CPGs, followed by an investigation of somatic events in matching tumor samples.

MAIN OUTCOMES AND MEASURES

Gene-level burden of rare germline PVs, identification of secondary somatic events accompanying the germline PVs, and characterization of less-explored types of rare germline PVs in RCC patients.

RESULTS

In clear cell RCC (n = 976 patients), patients exhibited significantly higher prevalence of PVs in VHL compared to controls (OR: 39.1, 95% CI: 7.01-218.07, p-value:4.95e-05, q-value:0.00584). In non-clear cell RCC (n = 380 patients), patients carried enriched burden of PVs in FH (OR: 77.9, 95% CI: 18.68-324.97, p-value:1.55e-08, q-value: 1.83e-06) and MET (OR: 1.98e11, 95% CI: 0-inf, p-value: 2.07e-05, q-value: 3.50e-07). In a CHEK2-focused analysis with European cases and controls, clear cell RCC patients (n=906 European patients) harbored nominal enrichment of the previously reported low-penetrance CHEK2 variants, p.Ile157Thr (OR:1.84, 95% CI: 1.00-3.36, p-value:0.049) and p.Ser428Phe (OR:5.20, 95% CI: 1.00-26.40, p-value:0.045) while non-clear cell RCC patients (n=295 European patients) exhibited nominal enrichment of CHEK2 LOF germline PVs (OR: 3.51, 95% CI: 1.10-11.10, p-value: 0.033). RCC patients with germline PVs in FH, MET, and VHL exhibited significantly earlier age of cancer onset compared to patients without any germline PVs in CPGs (Mean: 46.0 vs 60.2 years old, Tukey adjusted p-value < 0.0001), and more than half had secondary somatic events affecting the same gene (n=10/15, 66.7%, 95% CI: 38.7-87.0%). Conversely, patients with rare germline PVs in CHEK2 exhibited a similar age of disease onset to patients without any identified germline PVs in CPGs (Mean: 60.1 vs 60.2 years old, Tukey adjusted p-value: 0.99), and only 30.4% of the patients carried secondary somatic events in CHEK2 (n=7/23, 95% CI: 14.1-53.0%). Finally, rare pathogenic germline cryptic splice variants underexplored in RCC were identified in SDHA and TSC1, and rare pathogenic germline CNVs were found in 18 patients, including CNVs in FH, SDHA, and VHL.

CONCLUSIONS AND RELEVANCE

This systematic analysis supports the existing link between several RCC risk genes and elevated RCC risk manifesting in earlier age of RCC onset. Our analysis calls for caution when assessing the role of germline PVs in CHEK2 due to the burden of founder variants with varying population frequency in different ancestry groups. It also broadens the definition of the RCC germline landscape of pathogenicity to incorporate previously understudied types of germline variants, such as cryptic splice variants and CNVs.

A comparison between similarity matrices for principal component analysis to assess population stratification in sequenced genetic data sets

Genetic similarity matrices are commonly used to assess population substructure (PS) in genetic studies. Through simulation studies and by the application to whole-genome sequencing (WGS) data, we evaluate the performance of three genetic similarity matrices: the unweighted and weighted Jaccard similarity matrices and the genetic relationship matrix. We describe different scenarios that can create numerical pitfalls and lead to incorrect conclusions in some instances. We consider scenarios in which PS is assessed based on loci that are located across the genome ('globally') and based on loci from a specific genomic region ('locally'). We also compare scenarios in which PS is evaluated based on loci from different minor allele frequency bins: common (>5%), low-frequency (5-0.5%) and rare (<0.5%) single-nucleotide variations (SNVs). Overall, we observe that all approaches provide the best clustering performance when computed based on rare SNVs. The performance of the similarity matrices is very similar for common and low-frequency variants, but for rare variants, the unweighted Jaccard matrix provides preferable clustering features. Based on visual inspection and in terms of standard clustering metrics, its clusters are the densest and the best separated in the principal component analysis of variants with rare SNVs compared with the other methods and different allele frequency cutoffs. In an application, we assessed the role of rare variants on local and global PS, using WGS data from multiethnic Alzheimer's disease data sets and European or East Asian populations from the 1000 Genome Project.

Recent advances and challenges of rare variant association analysis in the biobank sequencing era

Causal variants for rare genetic diseases are often rare in the general population. Rare variants may also contribute to common complex traits and can have much larger per-allele effect sizes than common variants, although power to detect these associations can be limited. Sequencing costs have steadily declined with technological advancements, making it feasible to adopt whole-exome and whole-genome profiling for large biobank-scale sample sizes. These large amounts of sequencing data provide both opportunities and challenges for rare-variant association analysis. Herein, we review the basic concepts of rare-variant analysis methods, the current state-of-the-art methods in utilizing variant annotations or external controls to improve the statistical power, and particular challenges facing rare variant analysis such as accounting for population structure, extremely unbalanced case-control design. We also review recent advances and challenges in rare variant analysis for familial sequencing data and for more complex phenotypes such as survival data. Finally, we discuss other potential directions for further methodology investigation.

Importance of Including Non-European Populations in Large Human Genetic Studies to Enhance Precision Medicine

One goal of genomic medicine is to uncover an individual's genetic risk for disease, which generally requires data connecting genotype to phenotype, as done in genome-wide association studies (GWAS). While there may be clinical promise to employing prediction tools such as polygenic risk scores (PRS), it currently stands that individuals of non-European ancestry may not reap the benefits of genomic medicine because of underrepresentation in large-scale genetics studies. Here, we discuss why this inequity poses a problem for genomic medicine and the reasons for the low transferability of PRS across populations. We also survey the ancestry representation of published GWAS and investigate how estimates of ancestry diversity in GWASparticipants might be biased. We highlight the importance of expanding genetic research in Africa, one of the most underrepresented regions in human genomics research, and discuss issues of ethics, resources, and technology for equitable advancement of genomic medicine.

Clarifying the causes of consistent and inconsistent findings in genetics

As research in genetics has advanced, some findings have been unexpected or shown to be inconsistent between studies or datasets. The reasons these inconsistencies arise are complex. Results from genetic studies can be affected by various factors including statistical power, linkage disequilibrium, quality control, confounding and selection bias, as well as real differences from interactions and effect modifiers, which may be informative about the mechanisms of traits and disease. Statistical artefacts can manifest as differences between results but they can also conceal underlying differences, which implies that their critical examination is important for understanding the underpinnings of traits. In this review, we examine these factors and outline how they can be identified and conceptualised with structural causal models. We explain the consequences they have on genetic estimates, such as genetic associations, polygenic scores, family‐ and genome‐wide heritability, and describe methods to address them to aid in the estimation of true effects of genetic variation. Clarifying these factors can help researchers anticipate when results are likely to diverge and aid researchers' understanding of causal relationships between genes and complex traits.

Controlling for human population stratification in rare variant association studies

Population stratification is a confounder of genetic association studies. In analyses of rare variants, corrections based on principal components (PCs) and linear mixed models (LMMs) yield conflicting conclusions. Studies evaluating these approaches generally focused on limited types of structure and large sample sizes. We investigated the properties of several correction methods through a large simulation study using real exome data, and several within- and between-continent stratification scenarios. We considered different sample sizes, with situations including as few as 50 cases, to account for the analysis of rare disorders. Large samples showed that accounting for stratification was more difficult with a continental than with a worldwide structure. When considering a sample of 50 cases, an inflation of type-I-errors was observed with PCs for small numbers of controls (≤ 100), and with LMMs for large numbers of controls (≥ 1000). We also tested a novel local permutation method (LocPerm), which maintained a correct type-I-error in all situations. Powers were equivalent for all approaches pointing out that the key issue is to properly control type-I-errors. Finally, we found that power of analyses including small numbers of cases can be increased, by adding a large panel of external controls, provided an appropriate stratification correction was used.

Taking population stratification into account by local permutations in rare-variant association studies on small samples

Many methods for rare variant association studies require permutations to assess the significance of tests. Standard permutations assume that all individuals are exchangeable and do not take population stratification (PS), a known confounding factor in genetic studies, into account. We propose a novel strategy, LocPerm, in which individual phenotypes are permuted only with their closest ancestry-based neighbors. We performed a simulation study, focusing on small samples, to evaluate and compare LocPerm with standard permutations and classical adjustment on first principal components. Under the null hypothesis, LocPerm was the only method providing an acceptable type I error, regardless of sample size and level of stratification. The power of LocPerm was similar to that of standard permutation in the absence of PS, and remained stable in different PS scenarios. We conclude that LocPerm is a method of choice for taking PS and/or small sample size into account in rare variant association studies.

Exploring the Causality Between Hypothyroidism and Non-alcoholic Fatty Liver: A Mendelian Randomization Study

The etiology of non-alcoholic fatty liver disease (NAFLD) involves complex interaction of genetic and environmental factors. A large number of observational studies have shown that hypothyroidism contributes to a high risk of NAFLD. However, the exact causality is still unknown. Due to the progress of genome-wide association study (GWAS) and the discovery of Mendelian randomization (MR), it is possible to explore the causality between the two diseases. In this study, in order to research into the influence of intermediate phenotypes on outcome, nine independent genetic variants of hypothyroidism obtained from the GWAS were used as instrumental variables (IVs) to perform MR analysis on NAFLD. Since there was no heterogeneity between IVs (P = 0.70), a fixed-effects model was used. The correlation between hypothyroidism and NAFLD was evaluated by using inverse-variance weighted (IVW) method and weighted median method. Then the sensitivity test was analyzed. The results showed that there was a high OR (1.7578; 95%CI 1.1897–2.5970; P = 0.0046) and a low intercept (−0.095; P = 0.431). None of the genetic variants drove the overall result (P < 0.01). Simply, we proved for the first time that the risk of NAFLD increases significantly on patients with hypothyroidism. Furthermore, we explained possible causes of NAFLD caused by hypothyroidism.

Exploring the midline soft tissue surface changes from 12 to 15 years of age in three distinct country population cohorts

INTRODUCTION

Several studies have highlighted differences in the facial features in a White European population. Genetics appear to have a major influence on normal facial variation, and environmental factors are likely to have minor influences on face shape directly or through epigenetic mechanisms.

AIM

The aim of this longitudinal cohort study is to determine the rate of change in midline facial landmarks in three distinct homogenous population groups (Finnish, Latvian, and Welsh) from 12.8 to 15.3 years of age. This age range covers the pubertal growth period for the majority of boys and girls.

METHODS

A cohort of children aged 12 were monitored for facial growth in three countries [Finland (n = 60), Latvia (n = 107), and Wales (n = 96)]. Three-dimensional facial surface images were acquired (using either laser or photogrammetric methods) at regular intervals (6-12 months) for 4 years. Ethical approval was granted in each country. Nine midline landmarks were identified and the relative spatial positions of these surface landmarks were measured relative to the mid-endocanthion (men) over a 4-year period.

RESULTS

This study reports the children who attended 95 per cent of all scanning sessions (Finland 48 out of 60; Latvia 104 out of 107; Wales 50 out of 96). Considerable facial variation is seen for all countries and sexes. There are clear patterns of growth that show different magnitudes at different age groups for the different country groups, sexes, and facial parameters. The greatest single yearly growth rate (5.4 mm) was seen for Welsh males for men-pogonion distance at 13.6 years of age. Males exhibit greater rates of growth compared to females. These variations in magnitude and timings are likely to be influenced by genetic ancestry as a result of population migration.

CONCLUSION

The midline points are a simple and valid method to assess the relative spatial positions of facial surface landmarks. This study confirms previous reports on the subtle differences in facial shapes and sizes of male and female children in different populations and also highlights the magnitudes and timings of growth for various midline landmark distances to the men point.

Autozygosity mapping and time-to-spontaneous delivery in Norwegian parent-offspring trios

Parental genetic relatedness may lead to adverse health and fitness outcomes in the offspring. However, the degree to which it affects human delivery timing is unknown. We use genotype data from ≃25 000 parent-offspring trios from the Norwegian Mother, Father and Child Cohort Study to optimize runs of homozygosity (ROH) calling by maximizing the correlation between parental genetic relatedness and offspring ROHs. We then estimate the effect of maternal, paternal and fetal autozygosity and that of autozygosity mapping (common segments and gene burden test) on the timing of spontaneous onset of delivery. The correlation between offspring ROH using a variety of parameters and parental genetic relatedness ranged between −0.2 and 0.6, revealing the importance of the minimum number of genetic variants included in an ROH and the use of genetic distance. The optimized compared to predefined parameters showed a ≃45% higher correlation between parental genetic relatedness and offspring ROH. We found no evidence of an effect of maternal, paternal nor fetal overall autozygosity on spontaneous delivery timing. Yet, through autozygosity mapping, we identified three maternal loci TBC1D1, SIGLECs and EDN1 gene regions reducing the median time-to-spontaneous onset of delivery by ≃2–5% (P-value < 2.3 × 10⁻⁶). We also found suggestive evidence of a fetal locus at 3q22.2, near the RYK gene region (P-value = 2.0 × 10⁻⁶). Autozygosity mapping may provide new insights on the genetic determinants of delivery timing beyond traditional genome-wide association studies, but particular and rigorous attention should be given to ROH calling parameter selection.

Genetic Determinants of Antibody-Mediated Immune Responses to Infectious Diseases Agents: A Genome-Wide and HLA Association Study

BACKGROUND

Infectious diseases are causally related to a large array of noncommunicable diseases (NCDs). Identifying genetic determinants of infections and antibody-mediated immune responses may shed light on this relationship and provide therapeutic targets for drug and vaccine development.

METHODS

We used the UK biobank cohort of up to 10 000 serological measurements of infectious diseases and genome-wide genotyping. We used data on 13 pathogens to define 46 phenotypes: 15 seropositivity case-control phenotypes and 31 quantitative antibody measurement phenotypes. For each of these, we performed genome-wide association studies (GWAS) using the fastGWA linear mixed model package and human leukocyte antigen (HLA) classical allele and amino acid residue associations analyses using Lasso regression for variable selection.

RESULTS

We included a total of 8735 individuals for case-control phenotypes, and an average (range) of 4286 (276-8555) samples per quantitative analysis. Fourteen of the GWAS yielded a genome-wide significant (P < 5 ×10-8) locus at the major histocompatibility complex (MHC) on chromosome 6. Outside the MHC, we found a total of 60 loci, multiple associated with Epstein-Barr virus (EBV)-related NCDs (eg, RASA3, MED12L, and IRF4). FUT2 was also identified as an important gene for polyomaviridae. HLA analysis highlighted the importance of DRB1*09:01, DQB1*02:01, DQA1*01:02, and DQA1*03:01 in EBV serologies and of DRB1*15:01 in polyomaviridae.

CONCLUSIONS

We have identified multiple genetic variants associated with antibody immune response to 13 infections, many of which are biologically plausible therapeutic or vaccine targets. This may help prioritize future research and drug development.

Inborn errors of type I IFN immunity in patients with life-threatening COVID-19

Clinical outcome upon infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) ranges from silent infection to lethal coronavirus disease 2019 (COVID-19). We have found an enrichment in rare variants predicted to be loss-of-function (LOF) at the 13 human loci known to govern Toll-like receptor 3 (TLR3)- and interferon regulatory factor 7 (IRF7)-dependent type I interferon (IFN) immunity to influenza virus in 659 patients with life-threatening COVID-19 pneumonia relative to 534 subjects with asymptomatic or benign infection. By testing these and other rare variants at these 13 loci, we experimentally defined LOF variants underlying autosomal-recessive or autosomal-dominant deficiencies in 23 patients (3.5%) 17 to 77 years of age. We show that human fibroblasts with mutations affecting this circuit are vulnerable to SARS-CoV-2. Inborn errors of TLR3- and IRF7-dependent type I IFN immunity can underlie life-threatening COVID-19 pneumonia in patients with no prior severe infection.

Multilevel principal components analysis of three-dimensional facial growth in adolescents

BACKGROUND AND OBJECTIVES

The study of age-related facial shape changes across different populations and sexes requires new multivariate tools to disentangle different sources of variations present in 3D facial images. Here we wish to use a multivariate technique called multilevel principal components analysis (mPCA) to study three-dimensional facial growth in adolescents.

METHODS

These facial shapes were captured for Welsh and Finnish subjects (both male and female) at multiple ages from 12 to 17 years old (i.e., repeated-measures data). 1000 "dense" 3D points were defined regularly for each shape by using a deformable template via "meshmonk" software. A three-level model was used here, namely: level 1 (sex/ethnicity); level 2, all "subject" variations excluding sex, ethnicity, and age; and level 3, age. The technicalities underpinning the mPCA method are presented in Appendices.

RESULTS

Eigenvalues via mPCA predicted that: level 1 (ethnicity/sex) contained 7.9% of variation; level 2 contained 71.5%; and level 3 (age) contained 20.6%. The results for the eigenvalues via mPCA followed a similar pattern to those results of single-level PCA. Results for modes of variation made sense, where effects due to ethnicity, sex, and age were reflected in modes at appropriate levels of the model. Standardised scores at level 1 via mPCA showed much stronger differentiation between sex and ethnicity groups than results of single-level PCA. Results for standardised scores from both single-level PCA and mPCA at level 3 indicated that females had different average "trajectories" with respect to these scores than males, which suggests that facial shape matures in different ways for males and females. No strong evidence of differences in growth patterns between Finnish and Welsh subjects was observed.

CONCLUSIONS

mPCA results agree with existing research relating to the general process of facial changes in adolescents with respect to age quoted in the literature. They support previous evidence that suggests that males demonstrate larger changes and for a longer period of time compared to females, especially in the lower third of the face. These calculations are therefore an excellent initial test that multivariate multilevel methods such as mPCA can be used to describe such age-related changes for "dense" 3D point data.

Genome-wide Association Study for Vitamin D Levels Reveals 69 Independent Loci

We aimed to increase our understanding of the genetic determinants of vitamin D levels by undertaking a large-scale genome-wide association study (GWAS) of serum 25 hydroxyvitamin D (25OHD). To do so, we used imputed genotypes from 401,460 white British UK Biobank participants with available 25OHD levels, retaining single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) > 0.1% and imputation quality score > 0.3. We performed a linear mixed model GWAS on standardized log-transformed 25OHD, adjusting for age, sex, season of measurement, and vitamin D supplementation. These results were combined with those from a previous GWAS including 42,274 Europeans. In silico functional follow-up of the GWAS results was undertaken to identify enrichment in gene sets, pathways, and expression in tissues, and to investigate the partitioned heritability of 25OHD and its shared heritability with other traits. Using this approach, the SNP heritability of 25OHD was estimated to 16.1%. 138 conditionally independent SNPs were detected (p value < 6.6 × 10-9) among which 53 had MAF < 5%. Single variant association signals mapped to 69 distinct loci, among which 63 were previously unreported. We identified enrichment in hepatic and lipid metabolism gene pathways and enriched expression of the 25OHD genes in liver, skin, and gastrointestinal tissues. We observed partially shared heritability between 25OHD and socio-economic traits, a feature which may be mediated through time spent outdoors. Therefore, through a large 25OHD GWAS, we identified 63 loci that underline the contribution of genes outside the vitamin D canonical metabolic pathway to the genetic architecture of 25OHD.

The genetic history of France

The study of the genetic structure of different countries within Europe has provided significant insights into their demographic history and population structure. Although France occupies a particular location at the western part of Europe and at the crossroads of migration routes, few population genetic studies have been conducted so far with genome-wide data. In this study, we analyzed SNP-chip genetic data from 2184 individuals born in France who were enrolled in two independent population cohorts. Using FineSTRUCTURE, six different genetic clusters of individuals were found that were very consistent between the two cohorts. These clusters correspond closely to geographic, historical, and linguistic divisions of France, and contain different proportions of ancestry from Stone and Bronze Age populations. By modeling the relationship between genetics and geography using EEMS, we were able to detect gene flow barriers that are similar across the two cohorts and correspond to major rivers and mountain ranges. Estimations of effective population sizes also revealed very similar patterns in both cohorts with a rapid increase of effective population sizes over the last 150 generations similar to other European countries. A marked bottleneck is also consistently seen in the two datasets starting in the 14^th century when the Black Death raged in Europe. In conclusion, by performing the first exhaustive study of the genetic structure of France, we fill a gap in genetic studies of Europe that will be useful to medical geneticists, historians, and archeologists.

Population and evolutionary genetics of the PAH locus to uncover overdominance and adaptive mechanisms in phenylketonuria: Results from a multiethnic study

BACKGROUND

Phenylketonuria (PKU) is the most common inborn error of amino acid metabolism in Europe. The reasons underlying the high prevalence of heterozygous carriers are not clearly understood. We aimed to look for pathogenic PAH variant enrichment according to geographical areas and patients' ethnicity using a multiethnic nationwide cohort of patients with PKU in France. We subsequently appraised the population differentiation, balancing selection and the molecular evolutionary history of the PAH locus.

METHODS

The French nationwide PKU study included patients who have been referred at the national level to the University Hospital of Nancy, and for whom a molecular diagnosis of phenylketonuria was made by Sanger sequencing. We performed enrichment analyses by comparing alternative allele frequencies using Fisher's exact test with Bonferroni adjustment. We estimated the amount of genetic differentiation among populations using Wright's fixation index (Fst). To estimate the molecular evolutionary history of the PAH gene, we performed phylogenetic and evolutionary analyses using whole-genome and exome-sequencing data from healthy individuals and non-PKU patients, respectively. Finally, we used exome-wide association study to decipher potential genetic loci associated with population divergence on PAH.

FINDINGS

The study included 696 patients and revealed 132 pathogenic PAH variants. Three geographical areas showed significant enrichment for a pathogenic PAH variant: North of France (p.Arg243Leu), North-West of France (p.Leu348Val), and Mediterranean coast (p.Ala403Val). One PAH variant (p.Glu280Gln) was significantly enriched among North-Africans (OR = 23·23; 95% CI: 9·75-55·38). PAH variants exhibiting a strong genetic differentiation were significantly enriched in the 'Biopterin_H' domain (OR = 6·45; 95% CI: 1·99-20·84), suggesting a balancing selection pressure on the biopterin function of PAH. Phylogenetic and timetree analyses were consistent with population differentiation events on European-, African-, and Asian-ancestry populations. The five PAH variants most strongly associated with a high selection pressure were phylogenetically close and were located within the biopterin domain coding region of PAH or in its vicinity. Among the non-PAH loci potentially associated with population divergence, two reached exome-wide significance: SSPO (SCO-spondin) and DBH (dopamine beta-hydroxylase), involved in neuroprotection and metabolic adaptation, respectively.

INTERPRETATION

Our data provide evidence on the combination of evolutionary and adaptive events in populations with distinct ancestries, which may explain the overdominance of some genetic variants on PAH.

FUNDING

French National Institute of Health and Medical Research (INSERM) UMR_S 1256.

Genome-wide Association Studies in Ancestrally Diverse Populations: Opportunities, Methods, Pitfalls, and Recommendations

Genome-wide association studies (GWASs) have focused primarily on populations of European descent, but it is essential that diverse populations become better represented. Increasing diversity among study participants will advance our understanding of genetic architecture in all populations and ensure that genetic research is broadly applicable. To facilitate and promote research in multi-ancestry and admixed cohorts, we outline key methodological considerations and highlight opportunities, challenges, solutions, and areas in need of development. Despite the perception that analyzing genetic data from diverse populations is difficult, it is scientifically and ethically imperative, and there is an expanding analytical toolbox to do it well.

A Cautionary Note on the Effects of Population Stratification Under an Extreme Phenotype Sampling Design

Extreme phenotype sampling (EPS) is a popular study design used to reduce genotyping or sequencing costs. Assuming continuous phenotype data are available on a large cohort, EPS involves genotyping or sequencing only those individuals with extreme phenotypic values. Although this design has been shown to have high power to detect genetic effects even at smaller sample sizes, little attention has been paid to the effects of confounding variables, and in particular population stratification. Using extensive simulations, we demonstrate that the false positive rate under the EPS design is greatly inflated relative to a random sample of equal size or a “case-control”-like design where the cases are from one phenotypic extreme and the controls randomly sampled. The inflated false positive rate is observed even with allele frequency and phenotype mean differences taken from European population data. We show that the effects of confounding are not reduced by increasing the sample size. We also show that including the top principal components in a logistic regression model is sufficient for controlling the type 1 error rate using data simulated with a population genetics model and using 1,000 Genomes genotype data. Our results suggest that when an EPS study is conducted, it is crucial to adjust for all confounding variables. For genetic association studies this requires genotyping a sufficient number of markers to allow for ancestry estimation. Unfortunately, this could increase the costs of a study if sequencing or genotyping was only planned for candidate genes or pathways; the available genetic data would not be suitable for ancestry correction as many of the variants could have a true association with the trait.