A general approach for haplotype phasing across the full spectrum of relatedness

Reconstructing complex admixture history using a hierarchical model

Various methods have been proposed to reconstruct admixture histories by analyzing the length of ancestral chromosomal tracts, such as estimating the admixture time and number of admixture events. However, available methods do not explicitly consider the complex admixture structure, which characterizes the joining and mixing patterns of different ancestral populations during the admixture process, and instead assume a simplified one-by-one sequential admixture model. In this study, we proposed a novel approach that considers the non-sequential admixture structure to reconstruct admixture histories. Specifically, we introduced a hierarchical admixture model that incorporated four ancestral populations and developed a new method, called HierarchyMix, which uses the length of ancestral tracts and the number of ancestry switches along genomes to reconstruct the four-way admixture history. By automatically selecting the optimal admixture model using the Bayesian information criterion principles, HierarchyMix effectively estimates the corresponding admixture parameters. Simulation studies confirmed the effectiveness and robustness of HierarchyMix. We also applied HierarchyMix to Uyghurs and Kazakhs, enabling us to reconstruct the admixture histories of Central Asians. Our results highlight the importance of considering complex admixture structures and demonstrate that HierarchyMix is a useful tool for analyzing complex admixture events.

Immune-response 3'UTR alternative polyadenylation quantitative trait loci contribute to variation in human complex traits and diseases

Genome-wide association studies (GWASs) have identified thousands of non-coding variants that are associated with human complex traits and diseases. The analysis of such GWAS variants in different contexts and physiological states is essential for deciphering the regulatory mechanisms underlying human disease. Alternative polyadenylation (APA) is a key post-transcriptional modification for most human genes that substantially impacts upon cell behavior. Here, we mapped 9,493 3'-untranslated region APA quantitative trait loci in 18 human immune baseline cell types and 8 stimulation conditions (immune 3'aQTLs). Through the comparison between baseline and stimulation data, we observed the high responsiveness of 3'aQTLs to immune stimulation (response 3'aQTLs). Co-localization and mendelian randomization analyses of immune 3'aQTLs identified 678 genes where 3'aQTL are associated with variation in complex traits, 27.3% of which were derived from response 3'aQTLs. Overall, these analyses reveal the role of immune 3'aQTLs in the determination of complex traits, providing new insights into the regulatory mechanisms underlying disease etiologies.

A purebred South American breed showing high effective population size and independent breed ancestry: The Chilean Terrier

The Chilean Terrier is a known breed in Chile that has not been genetically assessed despite its distinctive color patterns, agility, and hardiness across the diversity of climates encountered within the Chilean landscape. The population structure and its relatedness with other breeds, as well as the actual origin of the breed, remain unknown. We estimated several population parameters using samples from individuals representing the distribution of the Chilean Terrier across the country. By utilizing the Illumina HD canine genotyping array, we computed the effective population size (Ne ), individual inbreeding, and relatedness to evaluate the genetic diversity of the breed. The results show that linkage disequilibrium was relatively low and decayed rapidly; in fact, Ne was very high when compared to other breeds, and similar to other American indigenous breeds (such as the Chihuahua with values of Ne near 500). These results are in line with the low estimates of genomic inbreeding and relatedness and the relatively large number of effective chromosome segments (Me  = 2467) obtained using the properties of the genomic relationship matrix. Between population analysis (cross-population extended haplotype homozygosity, di ) with other breeds such as the Jack Russell Terrier, the Peruvian-Inca Orchid, and the Chihuahua suggested that candidate regions harboring FGF5, PAX3, and ASIP, probably explained some morphological traits, such as the distinctive color pattern characteristic of the breed. When considering Admixture estimates and phylogenetic analysis, together with other breeds of American and European origin, the Chilean Terrier does not have a recent European ancestry. Overall, the results suggest that the breed has evolved independently in Chile from other terrier breeds, from an unknown European terrier ancestor.

Identifying genetic variants that influence the abundance of cell states in single-cell data

To understand genetic mechanisms driving disease, it is essential but difficult to map how risk alleles affect the composition of cells present in the body. Single-cell profiling quantifies granular information about tissues, but variant-associated cell states may reflect diverse combinations of the profiled cell features that are challenging to predefine. We introduce GeNA (Genotype-Neighborhood Associations), a statistical tool to identify cell state abundance quantitative trait loci (csaQTLs) in high-dimensional single-cell datasets. Instead of testing associations to predefined cell states, GeNA flexibly identifies the cell states whose abundance is most associated with genetic variants. In a genome-wide survey of scRNA-seq peripheral blood profiling from 969 individuals,1GeNA identifies five independent loci associated with shifts in the relative abundance of immune cell states. For example, rs3003-T (p=1.96×10-11) associates with increased abundance of NK cells expressing TNF-α response programs. This csaQTL colocalizes with increased risk for psoriasis, an autoimmune disease that responds to anti-TNF treatments. Flexibly characterizing csaQTLs for granular cell states may help illuminate how genetic background alters cellular composition to confer disease risk.

Eurasian back-migration into Northeast Africa was a complex and multifaceted process

Recent studies have identified Northeast Africa as an important area for human movements during the Holocene. Eurasian populations have moved back into Northeastern Africa and contributed to the genetic composition of its people. By gathering the largest reference dataset to date of Northeast, North, and East African as well as Middle Eastern populations, we give new depth to our knowledge of Northeast African demographic history. By employing local ancestry methods, we isolated the Non-African parts of modern-day Northeast African genomes and identified the best putative source populations. Egyptians and Sudanese Copts bore most similarities to Levantine populations whilst other populations in the region generally had predominantly genetic contributions from the Arabian peninsula rather than Levantine populations for their Non-African genetic component. We also date admixture events and investigated which factors influenced the date of admixture and find that major linguistic families were associated with the date of Eurasian admixture. Taken as a whole we detect complex patterns of admixture and diverse origins of Eurasian admixture in Northeast African populations of today.

Prediction of breast cancer risk for sisters of women attending screening

BACKGROUND

Risk assessment is important for breast cancer prevention and early detection. We aimed to examine whether common risk factors, mammographic features, and breast cancer risk prediction scores of a woman were associated with breast cancer risk for her sisters.

METHODS

We included 53 051 women from the Karolinska Mammography Project for Risk Prediction of Breast Cancer (KARMA) study. Established risk factors were derived using self-reported questionnaires, mammograms, and single nucleotide polymorphism genotyping. Using the Swedish Multi-Generation Register, we identified 32 198 sisters of the KARMA women (including 5352 KARMA participants and 26 846 nonparticipants). Cox models were used to estimate the hazard ratios of breast cancer for both women and their sisters, respectively.

RESULTS

A higher breast cancer polygenic risk score, a history of benign breast disease, and higher breast density in women were associated with an increased risk of breast cancer for both women and their sisters. No statistically significant association was observed between breast microcalcifications and masses in women and breast cancer risk for their sisters. Furthermore, higher breast cancer risk scores in women were associated with an increased risk of breast cancer for their sisters. Specifically, the hazard ratios for breast cancer per 1 standard deviation increase in age-adjusted KARMA, Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm (BOADICEA), and Tyrer-Cuzick risk scores were 1.16 (95% confidence interval [CI] = 1.07 to 1.27), 1.23 (95% CI = 1.12 to 1.35), and 1.21 (95% CI = 1.11 to 1.32), respectively.

CONCLUSION

A woman's breast cancer risk factors are associated with her sister's breast cancer risk. However, the clinical utility of these findings requires further investigation.

The parasitic nematode Strongyloides ratti exists predominantly as populations of long-lived asexual lineages

Nematodes are important parasites of people and animals, and in natural ecosystems they are a major ecological force. Strongyloides ratti is a common parasitic nematode of wild rats and we have investigated its population genetics using single-worm, whole-genome sequencing. We find that S. ratti populations in the UK consist of mixtures of mainly asexual lineages that are widely dispersed across a host population. These parasite lineages are likely very old and may have originated in Asia from where rats originated. Genes that underly the parasitic phase of the parasite's life cycle are hyperdiverse compared with the rest of the genome, and this may allow the parasites to maximise their fitness in a diverse host population. These patterns of parasitic nematode population genetics have not been found before and may also apply to Strongyloides spp. that infect people, which will affect how we should approach their control.

Genetic Origins and Adaptive Evolution of the Deng People on the Tibetan Plateau

The Tibetan Plateau is populated by diverse ethnic groups, but most of them are underrepresented in genomics studies compared with the Tibetans (TIB). Here, to gain further insight into the genetic diversity and evolutionary history of the people living in the Tibetan Plateau, we sequenced 54 whole genomes of the Deng people with high coverage (30-60×) and analyzed the data together with that of TIB and Sherpas, as well as 968 ancient Asian genomes and available archaic and modern human data. We identified 17.74 million novel single-nucleotide variants from the newly sequenced genomes, although the Deng people showed reduced genomic diversity and a relatively small effective population size. Compared with the other Tibetan highlander groups which are highly admixed, the Deng people are dominated by a sole ancestry that could be traced to some ancient northern East Asian populations. The divergence between Deng and Tibetan people (∼4,700-7,200 years) was more recent than that between highlanders and the Han Chinese (Deng-HAN, ∼9,000-14,000 years; TIB-HAN, 7,200-10,000 years). Adaptive genetic variants (AGVs) identified in the Deng are only partially shared with those previously reported in the TIB like HLA-DQB1, whereas others like KLHL12 were not reported in TIB. In contrast, the top candidate genes harboring AGVs as previously identified in TIB, like EPAS1 and EGLN1, do not show strong positive selection signals in Deng. Interestingly, Deng also showed a different archaic introgression scenario from that observed in the TIB. Our results suggest that convergent adaptation might be prevalent on the Tibetan Plateau.

Identifying crossovers and shared genetic material in whole genome sequencing data from families

Large, whole-genome sequencing (WGS) data sets containing families provide an important opportunity to identify crossovers and shared genetic material in siblings. However, the high variant calling error rates of WGS in some areas of the genome can result in spurious crossover calls, and the special inheritance status of the X Chromosome presents challenges. We have developed a hidden Markov model that addresses these issues by modeling the inheritance of variants in families in the presence of error-prone regions and inherited deletions. We call our method PhasingFamilies. We validate PhasingFamilies using the platinum genome family NA1281 (precision: 0.81; recall: 0.97), as well as simulated genomes with known crossover positions (precision: 0.93; recall: 0.92). Using 1925 quads from the Simons Simplex Collection, we found that PhasingFamilies resolves crossovers to a median resolution of 3527.5 bp. These crossovers recapitulate existing recombination rate maps, including for the X Chromosome; produce sibling pair IBD that matches expected distributions; and are validated by the haplotype estimation tool SHAPEIT. We provide an efficient, open-source implementation of PhasingFamilies that can be used to identify crossovers from family sequencing data.

Interaction of Mitochondrial Polygenic Score and Lifestyle Factors in LRRK2 p.Gly2019Ser Parkinsonism

BACKGROUND

A mitochondrial polygenic score (MGS) is composed of genes related to mitochondrial function and found to be associated with Parkinson's disease (PD) risk.

OBJECTIVE

To investigate the impact of the MGS and lifestyle/environment on age at onset (AAO) in LRRK2 p.Gly2019Ser parkinsonism (LRRK2-PD) and idiopathic PD (iPD).

METHODS

We included N = 486 patients with LRRK2-PD and N = 9259 with iPD from the Accelerating Medicines Partnership® Parkinson's Disease Knowledge Platform (AMP-PD), Fox Insight, and a Tunisian Arab-Berber founder population. Genotyping data were used to perform the MGS analysis. Additionally, lifestyle/environmental data were obtained from the PD Risk Factor Questionnaire (PD-RFQ). Linear regression models were used to assess the relationship between MGS, lifestyle/environment, and AAO.

RESULTS

Our derived MGS was significantly higher in PD cases compared with controls (P = 1.1 × 10-8 ). We observed that higher MGS was significantly associated with earlier AAO in LRRK2-PD (P = 0.047, β = -1.40) and there was the same trend with a smaller effect size in iPD (P = 0.231, β = 0.22). There was a correlation between MGS and AAO in LRRK2-PD patients of European descent (P = 0.049, r = -0.12) that was visibly less pronounced in Tunisians (P = 0.449, r = -0.05). We found that the MGS interacted with caffeinated soda consumption (P = 0.003, β = -5.65) in LRRK2-PD and with tobacco use (P = 0.010, β = 1.32) in iPD. Thus, patients with a high MGS had an earlier AAO only if they consumed caffeinated soda or were non-smokers.

CONCLUSIONS

The MGS was more strongly associated with earlier AAO in LRRK2-PD compared with iPD. Caffeinated soda consumption or tobacco use interacted with MGS to predict AAO. Our study suggests gene-environment interactions as modifiers of AAO in LRRK2-PD. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

A novel antisense lncRNA, ARBAG harboring an RNA destabilizing GWAS variant for C-peptide dictates the transcript isoforms of GABRA6 in cerebellum

Human disease-associated genetic variations often map to long non-coding RNA (lncRNA) genes; however, elucidation of their functional impact is challenging. We previously identified a new genetic variant rs4454083 (A/G) residing in exon of an uncharacterized lncRNA ARBAG that strongly associates with plasma levels of C-peptide, a hormone that regulates insulin bioavailability. On the opposite strand, rs4454083 also corresponds to an intron of a cerebellum-specific GABA receptor subunit gene GABRA6 that mediates strengthening of inhibitory synapses by insulin. Here, we show that alleles of rs4454083 modulate transcript levels of the antisense gene, ARBAG, which then controls the expression of the sense gene, GABRA6. Predisposing to low C-peptide, GG (a minor allele genotype across ethnicities) stabilizes ARBAG lncRNA causing higher transcript levels in cerebellum. ARBAG lncRNA abundance leads to cleavage of GABRA6 mRNA at the complementary region, resulting in a dysfunctional GABRA6 protein that would not be recruited for synapse strengthening. Together, our findings in human cerebellar cell-line and induced Pluripotent Stem Cells (iPSCs) demonstrate biological role of a novel lncRNA in determining the ratio of mRNA isoforms of a protein-coding gene and the ability of an embedded variant in modulating lncRNA stability leading to inter-individual differences in protein expression.

A Breast Cancer Polygenic Risk Score Is Feasible for Risk Stratification in the Norwegian Population

BACKGROUND

Statistical associations of numerous single nucleotide polymorphisms with breast cancer (BC) have been identified in genome-wide association studies (GWAS). Recent evidence suggests that a Polygenic Risk Score (PRS) can be a useful risk stratification instrument for a BC screening strategy, and a PRS test has been developed for clinical use. The performance of the PRS is yet unknown in the Norwegian population.

AIM

To evaluate the performance of PRS models for BC in a Norwegian dataset.

METHODS

We investigated a sample of 1053 BC cases and 7094 controls from different regions of Norway. PRS values were calculated using four PRS models, and their performance was evaluated by the area under the curve (AUC) and the odds ratio (OR). The effect of the PRS on the age of onset of BC was determined by a Cox regression model, and the lifetime absolute risk of developing BC was calculated using the iCare tool.

RESULTS

The best performing PRS model included 3820 SNPs, which yielded an AUC = 0.625 and an OR = 1.567 per one standard deviation increase. The PRS values of the samples correlate with an increased risk of BC, with a hazard ratio of 1.494 per one standard deviation increase (95% confidence interval of 1.406-1.588). The individuals in the highest decile of the PRS have at least twice the risk of developing BC compared to the individuals with a median PRS. The results in this study with Norwegian samples are coherent with the findings in the study conducted using Estonian and UK Biobank samples.

CONCLUSION

The previously validated PRS models have a similar observed accuracy in the Norwegian data as in the UK and Estonian populations. A PRS provides a meaningful association with the age of onset of BC and lifetime risk. Therefore, as suggested in Estonia, a PRS may also be integrated into the screening strategy for BC in Norway.

Genetic control of mRNA splicing as a potential mechanism for incomplete penetrance of rare coding variants

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-sequencing data in GTEx V8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased whole genome sequencing data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

Microevolutionary processes analysis in the Lithuanian genome

Differences in the relative fitness of genomic variants are foundational, without these, neither natural selection nor adaption can exist. This research analyzed two microevolutionary forces, mutations, and positive selection, using whole genome sequencing data from Lithuanians across three generations: newborns (generation I), their parents (generation II), 60 years old Lithuanians, and the root ancestors (generation III). The main objective was to determine the frequency of mutations under selection in modern humans and how allele frequencies change across generations. Our results show that going through all the landscapes of the relative fitness on each chromosome, the general relative fitness background pattern remains the same in analysed generations. However, the tendency of relative fitness to decrease, in general, is noted. We hypothesize that the de novo genome variants or genome variants with a very low frequency that formed in the previous generation did not have time to be as affected by natural selection, thus, in the following generation, the force of natural selection acting on them is greater and their cumulative relative fitness also decreases. The strong natural selection pressure on the genetic regions that encode the NEGR1 and PTPN1/PTNP21 genes were also identified, highlighting the evolution of the Lithuanian population's genome over generations, and possible genomic "deficiencies" for better adaptation.

Time-varying effects are common in genetic control of gestational duration

Preterm birth is a major burden to neonatal health worldwide, determined in part by genetics. Recently, studies discovered several genes associated with this trait or its continuous equivalent-gestational duration. However, their effect timing, and thus clinical importance, is still unclear. Here, we use genotyping data of 31 000 births from the Norwegian Mother, Father and Child cohort (MoBa) to investigate different models of the genetic pregnancy 'clock'. We conduct genome-wide association studies using gestational duration or preterm birth, replicating known maternal associations and finding one new fetal variant. We illustrate how the interpretation of these results is complicated by the loss of power when dichotomizing. Using flexible survival models, we resolve this complexity and find that many of the known loci have time-varying effects, often stronger early in pregnancy. The overall polygenic control of birth timing appears to be shared in the term and preterm, but not very preterm, periods and exploratory results suggest involvement of the major histocompatibility complex genes in the latter. These findings show that the known gestational duration loci are clinically relevant and should help design further experimental studies.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

History of tuberculosis disease is associated with genetic regulatory variation in Peruvians

A quarter of humanity is estimated to be latently infected with Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n=63) or did not progress to TB (controls, n=63). Transcriptomic profiling of monocyte-derived dendritic cells (DCs) and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Five genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

Invited review: Good practices in genome-wide association studies to identify candidate sequence variants in dairy cattle

Genotype data from dairy cattle selection programs have greatly facilitated GWAS to identify variants related to economic traits. Results can enhance the accuracy of genomic prediction, analyze more complex models that go beyond additive effects, elucidate the genetic architecture of a trait, and finally, decipher the underlying biology of traits. The entire process, comprising data generation, quality control, statistical analyses, interpretation of association results, and linking results to biology should be designed and executed to minimize the generation of false-positive and false-negative associations and misleading links to biological processes. This review aims to provide general guidelines for data analysis that address data quality control, association tests, adjustment for population stratification, and significance evaluation to improve the reliability of conclusions. We also provide guidance on post-GWAS strategy and the interpretation of results. These guidelines are tailored to dairy cattle, which are characterized by long-range linkage disequilibrium, large half-sib families, and routinely collected phenotypes, requiring different approaches than those applied in human GWAS. We discuss common limitations and challenges that have been overlooked in the analysis and interpretation of GWAS to identify candidate sequence variants in dairy cattle.

Germline de novo mutations in families with Mendelian cancer syndromes caused by defects in DNA repair

DNA repair defects underlie many cancer syndromes. We tested whether de novo germline mutations (DNMs) are increased in families with germline defects in polymerase proofreading or base excision repair. A parent with a single germline POLE or POLD1 mutation, or biallelic MUTYH mutations, had 3-4 fold increased DNMs over sex-matched controls. POLE had the largest effect. The DNMs carried mutational signatures of the appropriate DNA repair deficiency. No DNM increase occurred in offspring of MUTYH heterozygous parents. Parental DNA repair defects caused about 20-150 DNMs per child, additional to the ~60 found in controls, but almost all extra DNMs occurred in non-coding regions. No increase in post-zygotic mutations was detected, excepting a child with bi-allelic MUTYH mutations who was excluded from the main analysis; she had received chemotherapy and may have undergone oligoclonal haematopoiesis. Inherited DNA repair defects associated with base pair-level mutations increase DNMs, but phenotypic consequences appear unlikely.

Multi-ethnic imputation system (MI-System): a genotype imputation server for high-dimensional data

OBJECTIVE

Genotype imputation is a commonly used technique that infers un-typed variants into a study's genotype data, allowing better identification of causal variants in disease studies. However, due to overrepresentation of Caucasian studies, there's a lack of understanding of genetic basis of health-outcomes in other ethnic populations. Therefore, facilitating imputation of missing key-predictor-variants that can potentially improve a risk health-outcome prediction model, specifically for Asian ancestry, is of utmost relevance.

METHODS

We aimed to construct an imputation and analysis web-platform, that primarily facilitates, but is not limited to, genotype imputation on East-Asians. The goal is to provide a collaborative imputation platform for researchers in the public domain, towards rapidly and efficiently conducting accurate genotype imputation.

RESULTS

We present an online genotype imputation platform, Multi-ethnic Imputation System (MI-System) (https://misystem.cgm.ntu.edu.tw/), that offers users 3 established pipelines, SHAPEIT2-IMPUTE2, SHAPEIT4-IMPUTE5, and Beagle5.1 for conducting imputation analyses. In addition to 1000 Genomes and Hapmap3, a new customized Taiwan Biobank (TWB) reference panel, specifically created for Taiwanese-Chinese ancestry is provided. MI-System further offers functions to create customized reference panels to be used for imputation, conduct quality control, split whole genome data into chromosomes, and convert genome builds.

CONCLUSION

Users can upload their genotype data and perform imputation with minimum effort and resources. The utility functions further can be utilized to preprocess user uploaded data with easy clicks. MI-System potentially contributes to Asian-population genetics research, while eliminating the requirement for high performing computational resources and bioinformatics expertise. It will enable an increased pace of research and provide a knowledge-base for genetic carriers of complex diseases, therefore greatly enhancing patient-driven research.

STATEMENT OF SIGNIFICANCE

Multi-ethnic Imputation System (MI-System), primarily facilitates, but is not limited to, imputation on East-Asians, through 3 established prephasing-imputation pipelines, SHAPEIT2-IMPUTE2, SHAPEIT4-IMPUTE5, and Beagle5.1, where users can upload their genotype data and perform imputation and other utility functions with minimum effort and resources. A new customized Taiwan Biobank (TWB) reference panel, specifically created for Taiwanese-Chinese ancestry is provided. Utility functions include (a) create customized reference panels, (b) conduct quality control, (c) split whole genome data into chromosomes, and (d) convert genome builds. Users can also combine 2 reference panels using the system and use combined panels as reference to conduct imputation using MI-System.

Genome-wide aggregated trans-effects on risk of type 1 diabetes: A test of the "omnigenic" sparse effector hypothesis of complex trait genetics

The "omnigenic" hypothesis postulates that the polygenic effects of common SNPs on a typical complex trait are mediated through trans-effects on expression of a relatively sparse set of effector ("core") genes. We tested this hypothesis in a study of 4,964 cases of type 1 diabetes (T1D) and 7,497 controls by using summary statistics to calculate aggregated (excluding the HLA region) trans-scores for gene expression in blood. From associations of T1D with aggregated trans-scores, nine putative core genes were identified, of which three-STAT1, CTLA4 and FOXP3-are genes in which variants cause monogenic forms of autoimmune diabetes. Seven of these genes affect the activity of regulatory T cells, and two are involved in immune responses to microbial lipids. Four T1D-associated genomic regions could be identified as master regulators via trans-effects on gene expression. These results support the sparse effector hypothesis and reshape our understanding of the genetic architecture of T1D.

The making of the oral microbiome in Agta hunter-gatherers

Ecological and genetic factors have influenced the composition of the human microbiome during our evolutionary history. We analysed the oral microbiota of the Agta, a hunter-gatherer population where some members have adopted an agricultural diet. We show that age is the strongest factor modulating the microbiome, probably through immunosenescence since we identified an increase in the number of species classified as pathogens with age. We also characterised biological and cultural processes generating sexual dimorphism in the oral microbiome. A small subset of oral bacteria is influenced by the host genome, linking host collagen genes to bacterial biofilm formation. Our data also suggest that shifting from a fish/meat diet to a rice-rich diet transforms their microbiome, mirroring the Neolithic transition. All of these factors have implications in the epidemiology of oral diseases. Thus, the human oral microbiome is multifactorial and shaped by various ecological and social factors that modify the oral environment.

Automatic landmarking identifies new loci associated with face morphology and implicates Neanderthal introgression in human nasal shape

We report a genome-wide association study of facial features in >6000 Latin Americans based on automatic landmarking of 2D portraits and testing for association with inter-landmark distances. We detected significant associations (P-value <5 × 10-8) at 42 genome regions, nine of which have been previously reported. In follow-up analyses, 26 of the 33 novel regions replicate in East Asians, Europeans, or Africans, and one mouse homologous region influences craniofacial morphology in mice. The novel region in 1q32.3 shows introgression from Neanderthals and we find that the introgressed tract increases nasal height (consistent with the differentiation between Neanderthals and modern humans). Novel regions include candidate genes and genome regulatory elements previously implicated in craniofacial development, and show preferential transcription in cranial neural crest cells. The automated approach used here should simplify the collection of large study samples from across the world, facilitating a cosmopolitan characterization of the genetics of facial features.

Clinical case study meets population cohort: identification of a BRCA1 pathogenic founder variant in Orcadians

We multiply ascertained the BRCA1 pathogenic missense variant c.5207T > C; p.Val1736Ala (V1736A) in clinical investigation of breast and ovarian cancer families from Orkney in the Northern Isles of Scotland, UK. We sought to investigate the frequency and clinical relevance of this variant in those of Orcadian ancestry as an exemplar of the value of population cohorts in clinical care, especially in isolated populations. Oral history and birth, marriage and death registrations indicated genealogical linkage of the clinical cases to ancestors from the Isle of Westray, Orkney. Further clinical cases were identified through targeted testing for V1736A in women of Orcadian ancestry attending National Health Service (NHS) genetic clinics for breast and ovarian cancer family risk assessments. The variant segregates with female breast and ovarian cancer in clinically ascertained cases. Separately, exome sequence data from 2088 volunteer participants with three or more Orcadian grandparents, in the ORCADES research cohort, was interrogated to estimate the population prevalence of V1736A in Orcadians. The effects of the variant were assessed using Electronic Health Record (EHR) linkage. Twenty out of 2088 ORCADES research volunteers (~1%) carry V1736A, with a common haplotype around the variant. This allele frequency is ~480-fold higher than in UK Biobank participants. Cost-effectiveness of population screening for BRCA1 founder pathogenic variants has been demonstrated at a carrier frequency below the ~1% observed here. Thus we suggest that Orcadian women should be offered testing for the BRCA1 V1736A founder pathogenic variant, starting with those with known Westray ancestry.

Genome-Wide Meta-Analysis Identifies Multiple Novel Rare Variants to Predict Common Human Infectious Diseases Risk

Infectious diseases still threaten global human health, and host genetic factors have been indicated as determining risk factors for observed variations in disease susceptibility, severity, and outcome. We performed a genome-wide meta-analysis on 4624 subjects from the 10,001 Dalmatians cohort, with 14 infection-related traits. Despite a rather small number of cases in some instances, we detected 29 infection-related genetic associations, mostly belonging to rare variants. Notably, the list included the genes CD28, INPP5D, ITPKB, MACROD2, and RSF1, all of which have known roles in the immune response. Expanding our knowledge on rare variants could contribute to the development of genetic panels that could assist in predicting an individual's life-long susceptibility to major infectious diseases. In addition, longitudinal biobanks are an interesting source of information for identifying the host genetic variants involved in infectious disease susceptibility and severity. Since infectious diseases continue to act as a selective pressure on our genomes, there is a constant need for a large consortium of biobanks with access to genetic and environmental data to further elucidate the complex mechanisms behind host-pathogen interactions and infectious disease susceptibility.

Genetic effects on the timing of parturition and links to fetal birth weight

The timing of parturition is crucial for neonatal survival and infant health. Yet, its genetic basis remains largely unresolved. We present a maternal genome-wide meta-analysis of gestational duration (n = 195,555), identifying 22 associated loci (24 independent variants) and an enrichment in genes differentially expressed during labor. A meta-analysis of preterm delivery (18,797 cases, 260,246 controls) revealed six associated loci and large genetic similarities with gestational duration. Analysis of the parental transmitted and nontransmitted alleles (n = 136,833) shows that 15 of the gestational duration genetic variants act through the maternal genome, whereas 7 act both through the maternal and fetal genomes and 2 act only via the fetal genome. Finally, the maternal effects on gestational duration show signs of antagonistic pleiotropy with the fetal effects on birth weight: maternal alleles that increase gestational duration have negative fetal effects on birth weight. The present study provides insights into the genetic effects on the timing of parturition and the complex maternal-fetal relationship between gestational duration and birth weight.

Whole-genome re-sequencing provides key genomic insights in farmed Arctic charr (Salvelinus alpinus) populations of anadromous and landlocked origin from Scandinavia

Arctic charr (Salvelinus alpinus) is a niche-market high-value species for Nordic aquaculture. Similar to other salmonids, both anadromous and landlocked populations are encountered. Whole-genome re-sequencing (22X coverage) was performed on two farmed populations of anadromous (Sigerfjord; n = 24) and landlocked (Arctic Superior; n = 24) origin from Norway and Sweden respectively. More than 5 million SNPs were used to study their genetic diversity and to scan for selection signatures. The two populations were clearly distinguished through principal component analysis, with the mean fixation index being ~0.12. Furthermore, the levels of genomic inbreeding estimated from runs of homozygosity were 6.23% and 8.66% for the Norwegian and the Swedish population respectively. Biological processes that could be linked to selection pressure associated primarily with the anadromous background and/or secondarily with domestication were suggested. Overall, our study provided insights regarding the genetic composition of two main strains of farmed Arctic charr from Scandinavia. At the same time, ample genomic resources were produced in the magnitude of millions of SNPs that could assist the transition of Nordic Arctic charr farming in the genomics era.

Cold-induced vasodilation response in a Japanese cohort: insights from cold-water immersion and genome-wide association studies

BACKGROUND

Cold-induced vasodilation (CIVD) occurs after blood vessels in the skin are constricted due to local cold exposure. Although many CIVD studies have been conducted, the underlying molecular mechanisms are yet to be clarified. Therefore, we explored genetic variants associated with CIVD response using the largest-scale dataset reported to date in a CIVD study involving wavelet analysis; thus, the findings improve our understanding of the molecular mechanisms that regulate the CIVD response.

METHODS

We performed wavelet analysis of three skin blood flow signals [endothelial nitric oxide (eNO)-independent, eNO-dependent, and neurogenic activities] during finger cold-water immersion at 5 °C in 94 Japanese young adults. Additionally, we conducted genome-wide association studies of CIVD using saliva samples collected from the participants.

RESULTS

We found that the mean wavelet amplitudes of eNO-independent and neurogenic activities significantly increased and decreased prior to CIVD, respectively. Our results also implied that as many as ~ 10% of the Japanese subjects did not show an apparent CIVD response. Our genome-wide association studies of CIVD using ~ 4,040,000 imputed data found no apparent CIVD-related genetic variants; however, we identified 10 genetic variants, including 2 functional genes (COL4A2 and PRLR) that are associated with notable blunted eNO-independent and neurogenic activity responses in individuals without CIVD response during local cold exposure.

CONCLUSIONS

Our findings indicate that individuals without CIVD response differentiated by genotypes with COL4A2 and PRLR genetic variants exhibited notable blunted eNO-independent and neurogenic activity responses during local cold exposure.

Haplotype-aware analysis of somatic copy number variations from single-cell transcriptomes

Genome instability and aberrant alterations of transcriptional programs both play important roles in cancer. Single-cell RNA sequencing (scRNA-seq) has the potential to investigate both genetic and nongenetic sources of tumor heterogeneity in a single assay. Here we present a computational method, Numbat, that integrates haplotype information obtained from population-based phasing with allele and expression signals to enhance detection of copy number variations from scRNA-seq. Numbat exploits the evolutionary relationships between subclones to iteratively infer single-cell copy number profiles and tumor clonal phylogeny. Analysis of 22 tumor samples, including multiple myeloma, gastric, breast and thyroid cancers, shows that Numbat can reconstruct the tumor copy number profile and precisely identify malignant cells in the tumor microenvironment. We identify genetic subpopulations with transcriptional signatures relevant to tumor progression and therapy resistance. Numbat requires neither sample-matched DNA data nor a priori genotyping, and is applicable to a wide range of experimental settings and cancer types.

Identification and characterization of the long non-coding RNA NFIA-AS2 as a novel locus for body mass index in American Indians

BACKGROUND

Genome-wide association studies have shown that body mass index (BMI), an estimate of obesity, is highly polygenic. Individual variants typically have small effect sizes, making it challenging to identify unique loci in under-represented ethnic groups which lack statistical power due to their small sample size. Yet obesity is a major health disparity and is particularly prevalent in southwestern American Indians. Here, we identify and characterize a new locus for BMI that was detected by analyzing moderate associations with BMI obtained in a population-based sample of southwestern American Indians together with the well-powered GIANT dataset.

METHODS

Genotypes for 10.5 million variants were tested for association with BMI in 5870 American Indians and 2600 variants that showed an association P < 10^-3 in the American Indian sample were combined in a meta-analysis with the BMI data reported in GIANT (N = 240,608). The newly identified gene, NFIA-AS2 was functionally characterized, and the impact of its lead associated variant rs1777538 was studied both in-silico and in-vitro.

RESULTS

Rs1777538 (T/C; C allele frequency = 0.16 in American Indians and 0.04 in GIANT, meta-analysis P = 5.0 × 10^-7) exhibited a large effect in American Indians (1 kg/m² decrease in BMI per copy of C allele). NFIA-AS2 was found to be a nuclear localized long non-coding RNA expressed in tissues pertinent to human obesity. Analysis of this variant in human brown preadipocytes showed that NFIA-AS2 transcripts carrying the C allele had increased RNA degradation compared to the T allele transcripts (half-lives = 9 h, 13 h respectively). During brown adipogenesis, NFIA-AS2 featured a stage-specific regulation of nearby gene expression where rs1777538 demonstrated an allelic difference in regulation in the mature adipocytes (the strongest difference was observed for L1TD1, P = 0.007).

CONCLUSION

Our findings support a role for NFIA-AS2 in regulating pathways that impact BMI.

Time-varying effects are common in genetic control of gestational duration

Preterm birth is a major burden to neonatal health worldwide, determined in part by genetics. Recently, studies discovered several genes associated with this trait or its continuous equivalent - gestational duration. However, their effect timing, and thus clinical importance, is still unclear. Here, we use genotyping data of 31,000 births from the Norwegian Mother, Father and Child cohort (MoBa) to investigate different models of the genetic pregnancy "clock". We conduct genome-wide association studies using gestational duration or preterm birth, replicating known maternal associations and finding one new foetal variant. We illustrate how the interpretation of these results is complicated by the loss of power when dichotomizing. Using flexible survival models, we resolve this complexity and find that many of the known loci have time-varying effects, often stronger early in pregnancy. The overall polygenic control of birth timing appears to be shared in the term and preterm, but not very preterm periods, and exploratory results suggest involvement of the major histocompatibility complex genes in the latter. These findings show that the known gestational duration loci are clinically relevant, and should help design further experimental studies.

Genetic control of mRNA splicing as a potential mechanism for incomplete penetrance of rare coding variants

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-seq data in GTEx v8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased WGS data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

Constructing a polygenic risk score for childhood obesity using functional data analysis

Obesity is a highly heritable condition that affects increasing numbers of adults and, concerningly, of children. However, only a small fraction of its heritability has been attributed to specific genetic variants. These variants are traditionally ascertained from genome-wide association studies (GWAS), which utilize samples with tens or hundreds of thousands of individuals for whom a single summary measurement (e.g., BMI) is collected. An alternative approach is to focus on a smaller, more deeply characterized sample in conjunction with advanced statistical models that leverage longitudinal phenotypes. Novel functional data analysis (FDA) techniques are used to capitalize on longitudinal growth information from a cohort of children between birth and three years of age. In an ultra-high dimensional setting, hundreds of thousands of single nucleotide polymorphisms (SNPs) are screened, and selected SNPs are used to construct two polygenic risk scores (PRS) for childhood obesity using a weighting approach that incorporates the dynamic and joint nature of SNP effects. These scores are significantly higher in children with (vs. without) rapid infant weight gain-a predictor of obesity later in life. Using two independent cohorts, it is shown that the genetic variants identified in very young children are also informative in older children and in adults, consistent with early childhood obesity being predictive of obesity later in life. In contrast, PRSs based on SNPs identified by adult obesity GWAS are not predictive of weight gain in the cohort of young children. This provides an example of a successful application of FDA to GWAS. This application is complemented with simulations establishing that a deeply characterized sample can be just as, if not more, effective than a comparable study with a cross-sectional response. Overall, it is demonstrated that a deep, statistically sophisticated characterization of a longitudinal phenotype can provide increased statistical power to studies with relatively small sample sizes; and shows how FDA approaches can be used as an alternative to the traditional GWAS.

Haplotyping-Assisted Diploid Assembly and Variant Detection with Linked Reads

Phasing is essential for determining the origins of each set of alleles in the whole-genome sequencing data of individuals. As such, it provides essential information for the causes of hereditary diseases and the sources of individual variability. Recent technical breakthroughs in linked-read (referred to as co-barcoding in other chapters of the book) and long-read sequencing and downstream analysis have brought the goal of accurate and complete phasing within reach. Here we review recent progress related to the assembly and phasing of personal genomes based on linked-reads and related applications. Motivated by current limitations in generating high-quality diploid assemblies and detecting variants, a new suite of software tools, Aquila, was developed to fully take advantage of linked-read sequencing technology. The overarching goal of Aquila is to exploit the strengths of linked-read technology including long-range connectivity and inherent phasing of variants for reference-assisted local de novo assembly at the whole-genome scale. The diploid nature of the assemblies facilitates detection and phasing of genetic variation, including single nucleotide variations (SNVs), small insertions and deletions (indels), and structural variants (SVs). An extension of Aquila, Aquila_stLFR, focuses on another newly developed linked-reads sequencing technology, single-tube long-fragment read (stLFR). AquilaSV, a region-based diploid assembly approach, is used to characterize structural variants and can achieve diploid assembly in one target region at a time. Lastly, we introduce HAPDeNovo, a program that exploits phasing information from linked-read sequencing to improve detection of de novo mutations. Use of these tools is expected to harness the advantages of linked-reads technology, improve phasing, and advance variant discovery.

Pancan-MNVQTLdb: systematic identification of multi-nucleotide variant quantitative trait loci in 33 cancer types

Multi-nucleotide variants (MNVs) are defined as clusters of two or more nearby variants existing on the same haplotype in an individual. Recent studies have identified millions of MNVs in human populations, but their functions remain largely unknown. Numerous studies have demonstrated that single-nucleotide variants could serve as quantitative trait loci (QTLs) by affecting molecular phenotypes. Therefore, we propose that MNVs can also affect molecular phenotypes by influencing regulatory elements. Using the genotype data from The Cancer Genome Atlas (TCGA), we first identified 223 759 unique MNVs in 33 cancer types. Then, to decipher the functions of these MNVs, we investigated the associations between MNVs and six molecular phenotypes, including coding gene expression, miRNA expression, lncRNA expression, alternative splicing, DNA methylation and alternative polyadenylation. As a result, we identified 1 397 821 cis-MNVQTLs and 402 381 trans-MNVQTLs. We further performed survival analysis and identified 46 173 MNVQTLs associated with patient overall survival. We also linked the MNVQTLs to genome-wide association studies (GWAS) data and identified 119 762 MNVQTLs that overlap with existing GWAS loci. Finally, we developed Pancan-MNVQTLdb (http://gong_lab.hzau.edu.cn/mnvQTLdb/) for data retrieval and download. Pancan-MNVQTLdb will help decipher the functions of MNVs in different cancer types and be an important resource for genetic and cancer research.

Domestication and selection footprints in Persian walnuts (Juglans regia)

Walnut (Juglans) species are economically important hardwood trees cultivated worldwide for both edible nuts and high-quality wood. Broad-scale assessments of species diversity, evolutionary history, and domestication are needed to improve walnut breeding. In this study, we sequenced 309 walnut accessions from around the world, including 55 Juglans relatives, 98 wild Persian walnuts (J. regia), 70 J. regia landraces, and 86 J. regia cultivars. The phylogenetic tree indicated that J. regia samples (section Dioscaryon) were monophyletic within Juglans. The core areas of genetic diversity of J. regia germplasm were southwestern China and southern Asia near the Qinghai-Tibet Plateau and the Himalayas, and the uplift of the Himalayas was speculated to be the main factor leading to the current population dynamics of Persian walnut. The pattern of genomic variation in terms of nucleotide diversity, linkage disequilibrium, single nucleotide polymorphisms, and insertions/deletions revealed the domestication and selection footprints in Persian walnut. Selective sweep analysis, GWAS, and expression analysis further identified two transcription factors, JrbHLH and JrMYB6, that influence the thickness of the nut diaphragm as loci under selection during domestication. Our results elucidate the domestication and selection footprints in Persian walnuts and provide a valuable resource for the genomics-assisted breeding of this important crop.

A Pipeline for Phasing and Genotype Imputation on Mixed Human Data (Parents-Offspring Trios and Unrelated Subjects) by Reviewing Current Methods and Software

Genotype imputation has become an essential prerequisite when performing association analysis. It is a computational technique that allows us to infer genetic markers that have not been directly genotyped, thereby increasing statistical power in subsequent association studies, which consequently has a crucial impact on the identification of causal variants. Many features need to be considered when choosing the proper algorithm for imputation, including the target sample on which it is performed, i.e., related individuals, unrelated individuals, or both. Problems could arise when dealing with a target sample made up of mixed data, composed of both related and unrelated individuals, especially since the scientific literature on this topic is not sufficiently clear. To shed light on this issue, we examined existing algorithms and software for performing phasing and imputation on mixed human data from SNP arrays, specifically when related subjects belong to trios. By discussing the advantages and limitations of the current algorithms, we identified LD-based methods as being the most suitable for reconstruction of haplotypes in this specific context, and we proposed a feasible pipeline that can be used for imputing genotypes in both phased and unphased human data.

Identity-by-descent analysis of a large Tourette's syndrome pedigree from Costa Rica implicates genes involved in neuronal development and signal transduction

Tourette Syndrome (TS) is a heritable, early-onset neuropsychiatric disorder that typically begins in early childhood. Identifying rare genetic variants that make a significant contribution to risk in affected families may provide important insights into the molecular aetiology of this complex and heterogeneous syndrome. Here we present a whole-genome sequencing (WGS) analysis from the 11-generation pedigree (>500 individuals) of a densely affected Costa Rican family which shares ancestry from six founder pairs. By conducting an identity-by-descent (IBD) analysis using WGS data from 19 individuals from the extended pedigree we have identified putative risk haplotypes that were not seen in controls, and can be linked with four of the six founder pairs. Rare coding and non-coding variants present on the haplotypes and only seen in haplotype carriers show an enrichment in pathways such as regulation of locomotion and signal transduction, suggesting common mechanisms by which the haplotype-specific variants may be contributing to TS-risk in this pedigree. In particular we have identified a rare deleterious missense variation in RAPGEF1 on a chromosome 9 haplotype and two ultra-rare deleterious intronic variants in ERBB4 and IKZF2 on the same chromosome 2 haplotype. All three genes play a role in neurodevelopment. This study, using WGS data in a pedigree-based approach, shows the importance of investigating both coding and non-coding variants to identify genes that may contribute to disease risk. Together, the genes and variants identified on the IBD haplotypes represent biologically relevant targets for investigation in other pedigree and population-based TS data.

Long-range chromosomal interactions increase and mark repressed gene expression during adipogenesis

Obesity perturbs central functions of human adipose tissue, centred on differentiation of preadipocytes to adipocytes, i.e., adipogenesis. The large environmental component of obesity makes it important to elucidate epigenetic regulatory factors impacting adipogenesis. Promoter Capture Hi-C (pCHi-C) has been used to identify chromosomal interactions between promoters and associated regulatory elements. However, long range interactions (LRIs) greater than 1 Mb are often filtered out of pCHi-C datasets, due to technical challenges and their low prevalence. To elucidate the unknown role of LRIs in adipogenesis, we investigated preadipocyte differentiation to adipocytes using pCHi-C and bulk and single nucleus RNA-seq data. We first show that LRIs are reproducible between biological replicates, and they increase >2-fold in frequency across adipogenesis. We further demonstrate that genomic loci containing LRIs are more epigenetically repressed than regions without LRIs, corresponding to lower gene expression in the LRI regions. Accordingly, as preadipocytes differentiate into adipocytes, LRI regions are more likely to contain repressed preadipocyte marker genes; whereas these same LRI regions are depleted of actively expressed adipocyte marker genes. Finally, we show that LRIs can be used to restrict multiple testing of the long-range cis-eQTL analysis to identify variants that regulate genes via LRIs. We exemplify this by identifying a putative long range cis regulatory mechanism at the LYPLAL1/TGFB2 obesity locus. In summary, we identify LRIs that mark repressed regions of the genome, and these interactions increase across adipogenesis, pinpointing developmental regions that need to be repressed in a cell-type specific way for adipogenesis to proceed.

Genomic prediction of carcass traits using different haplotype block partitioning methods in beef cattle

Genomic prediction (GP) based on haplotype alleles can capture quantitative trait loci (QTL) effects and increase predictive ability because the haplotypes are expected to be in linkage disequilibrium (LD) with QTL. In this study, we constructed haploblocks using LD-based and the fixed number of single nucleotide polymorphisms (fixed-SNP) methods with Illumina BovineHD chip in beef cattle. To evaluate the performance of different haplotype block partitioning methods, we constructed haploblocks based on LD thresholds (from r 2 > 0.2 to r 2 > 0.8) and the number of fixed-SNPs (5, 10, 20). The performance of predictive methods for three carcass traits including liveweight (LW), dressing percentage (DP), and longissimus dorsi muscle weight (LDMW) was evaluated using three approaches (GBLUP and BayesB model based on the SNP, GHBLUP, and BayesBH models based on the haploblock, and GHBLUP+GBLUP and BayesBH+BayesB models based on the combined haploblock and the nonblocked SNPs, which were located between blocks). In this study, we found the accuracies of LD-based and fixed-SNP haplotype Bayesian methods outperformed the Bayesian models (up to 8.54 ± 7.44% and 5.74 ± 2.95%, respectively). GHBLUP showed a high improvement (up to 11.29 ± 9.87%) compared with GBLUP. The Bayesian models have higher accuracies than BLUP models in most scenarios. The average computing time of the BayesBH+BayesB model can reduce by 29.3% compared with the BayesB model. The prediction accuracies using the LD-based haplotype method showed higher improvements than the fixed-SNP haplotype method. In addition, to avoid the influence of rare haplotypes generated from haplotype construction, we compared the performance of GP by filtering four types of minor haplotype allele frequency (MHAF) (0.01, 0.025, 0.05, and 0.1) under different conditions (LD levels were set at r 2 > 0.3, and the fixed number of SNPs was 5). We found the optimal MHAF threshold for LW was 0.01, and the optimal MHAF threshold for DP and LDMW was 0.025.

ARHap: Association Rule Haplotype Phasing

This article proposes a novel approach for Individual Human phasing through discovery of interesting hidden relations among single variant sites. The proposed framework, called ARHap, learns strong association rules among variant loci on the genome and develops a combinatorial approach for fast and accurate haplotype phasing based on the discovered associations. ARHap is composed of two main modules or processing phases. In the first phase, called association rule learning, ARHap identifies quantitative association rules from a collection of DNA reads of the organism under study, resulting in a set of strong rules that reveal the inter-dependency of alleles. In the next phase, called haplotype reconstruction, we develop algorithms to utilize the learned rules to construct highly reliable haplotypes at individual single nucleotide polymorphism (SNP) sites. ARHap has several features that lead to both fast and accurate haplotyping. It uses an incremental haplotype reconstruction approach that enables us to generate association rules according to the unreconstructed SNP sites during each round of the algorithm. During each round, the association rule learning module generates rules while constraining the length of the rules and limiting the rules to those that contribute to reconstruction of unreconstructed sites only. The framework begins by generating rules of small size and highly strong. The rule length can increase and/or criteria about strongness of the rule are adjusted gradually, during subsequent rounds, if some SNP sites have remained unreconstructed. This adaptive approach, which uses feedback from haplotype reconstruction module, eliminates generation of rules that do not contribute to haplotype reconstruction as well as weak rules that may introduce error in the final haplotypes. Extensive experimental analyses on datasets representing diploid organisms demonstrate superiority of ARHap in diploid haplotyping compared to the state-of-the-art algorithms. In particular, we show that this novel approach to haplotype phasing not only is fast but also achieves significantly better accuracy performance compared to other read-based computational approaches.

Ancestry-related distribution of Runs of homozygosity and functional variants in Qatari population

Background

Describing how genetic history shapes the pattern of medically relevant variants could improve the understanding of how specific loci interact with each other and affect diseases and traits prevalence. The Qatari population is characterized by a complex history of admixture and substructure, and the study of its population genomic features would provide valuable insights into the genetic landscape of functional variants. Here, we analyzed the genomic variation of 186 newly-genotyped healthy individuals from the Qatari peninsula.

Results

We discovered an intricate genetic structure using ancestry related analyses. In particular, the presence of three different clusters, Cluster 1, Cluster 2 and Cluster 3 (with Near Eastern, South Asian and African ancestry, respectively), was detected with an additional fourth one (Cluster 4) with East Asian ancestry. These subpopulations show differences in the distribution of runs of homozygosity (ROH) and admixture events in the past, ranging from 40 to 5 generations ago. This complex genetic history led to a peculiar pattern of functional markers under positive selection, differentiated in shared signals and private signals. Interestingly we found several signatures of shared selection on SNPs in the FADS2 gene, hinting at a possible common evolutionary link to dietary intake. Among the private signals, we found enrichment for markers associated with HDL and LDL for Cluster 1(Near Eastern ancestry) and Cluster 3 (South Asian ancestry) and height and blood traits for Cluster 2 (African ancestry).

The differences in genetic history among these populations also resulted in the different frequency distribution of putative loss of function variants. For example, homozygous carriers for rs2884737, a variant linked to an anticoagulant drug (warfarin) response, are mainly represented by individuals with predominant Bedouin ancestry (risk allele frequency G at 0.48).

Conclusions

We provided a detailed catalogue of the different ancestral pattern in the Qatari population highlighting differences and similarities in the distribution of selected variants and putative loss of functions. Finally, these results would provide useful guidance for assessing genetic risk factors linked to consanguinity and genetic ancestry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01087-1.

Variants influencing age at diagnosis of HNF1A-MODY

BACKGROUND

HNF1A-MODY is a monogenic form of diabetes caused by variants in the HNF1A gene. Different HNF1A variants are associated with differences in age of disease onset, but other factors are postulated to influence this trait. Here, we searched for genetic variants influencing age of HNF1A-MODY onset.

METHODS

Blood samples from 843 HNF1A-MODY patients from Czech Republic, France, Poland, Slovakia, the UK and the US were collected. A validation set consisted of 121 patients from the US. We conducted a genome-wide association study in 843 HNF1A-MODY patients. Samples were genotyped using Illumina Human Core arrays. The core analysis was performed using the GENESIS package in R statistical software. Kinship coefficients were estimated with the KING and PC-Relate algorithms. In the linear mixed model, we accounted for year of birth, sex, and location of the HNF1A causative variant.

RESULTS

A suggestive association with age of disease onset was observed for rs2305198 (p = 2.09E-07) and rs7079157 (p = 3.96E-06) in the HK1 gene, rs2637248 in the LRMDA gene (p = 2.44E-05), and intergenic variant rs2825115 (p = 2.04E-05). Variant rs2637248 reached nominal significance (p = 0.019), while rs7079157 (p = 0.058) and rs2825115 (p = 0.068) showed suggestive association with age at diabetes onset in the validation set.

CONCLUSIONS

rs2637248 in the LRMDA gene is associated with age at diabetes onset in HNF1A-MODY patients.

High-coverage whole-genome sequencing of the expanded 1000 Genomes Project cohort including 602 trios

The 1000 Genomes Project (1kGP) is the largest fully open resource of whole-genome sequencing (WGS) data consented for public distribution without access or use restrictions. The final, phase 3 release of the 1kGP included 2,504 unrelated samples from 26 populations and was based primarily on low-coverage WGS. Here, we present a high-coverage 3,202-sample WGS 1kGP resource, which now includes 602 complete trios, sequenced to a depth of 30X using Illumina. We performed single-nucleotide variant (SNV) and short insertion and deletion (INDEL) discovery and generated a comprehensive set of structural variants (SVs) by integrating multiple analytic methods through a machine learning model. We show gains in sensitivity and precision of variant calls compared to phase 3, especially among rare SNVs as well as INDELs and SVs spanning frequency spectrum. We also generated an improved reference imputation panel, making variants discovered here accessible for association studies.

Sodium-calcium exchanger-3 regulates pain "wind-up": From human psychophysics to spinal mechanisms

Repeated application of noxious stimuli leads to a progressively increased pain perception; this temporal summation is enhanced in and predictive of clinical pain disorders. Its electrophysiological correlate is "wind-up," in which dorsal horn spinal neurons increase their response to repeated nociceptor stimulation. To understand the genetic basis of temporal summation, we undertook a GWAS of wind-up in healthy human volunteers and found significant association with SLC8A3 encoding sodium-calcium exchanger type 3 (NCX3). NCX3 was expressed in mouse dorsal horn neurons, and mice lacking NCX3 showed normal, acute pain but hypersensitivity to the second phase of the formalin test and chronic constriction injury. Dorsal horn neurons lacking NCX3 showed increased intracellular calcium following repetitive stimulation, slowed calcium clearance, and increased wind-up. Moreover, virally mediated enhanced spinal expression of NCX3 reduced central sensitization. Our study highlights Ca2+ efflux as a pathway underlying temporal summation and persistent pain, which may be amenable to therapeutic targeting.

The recombination landscape of the Khoe-San likely represents the upper limits of recombination divergence in humans

BACKGROUND

Recombination maps are  important resources for epidemiological and evolutionary analyses; however, there are currently no recombination maps representing any African population outside of those with West African ancestry. We infer the demographic history for the Nama, an indigenous Khoe-San population of southern Africa, and derive a novel, population-specific recombination map from the whole genome sequencing of 54 Nama individuals. We hypothesise that there are no publicly available recombination maps representative of the Nama, considering the deep population divergence and subsequent isolation of the Khoe-San from other African groups.

RESULTS

We show that the recombination landscape of the Nama does not cluster with any continental groups with publicly available representative recombination maps. Finally, we use selection scans as an example of how fine-scale differences between the Nama recombination map and the combined Phase II HapMap recombination map can impact the outcome of selection scans.

CONCLUSIONS

Fine-scale differences in recombination can meaningfully alter the results of a selection scan. The recombination map we infer likely represents an upper bound on the extent of divergence we expect to see for a recombination map in humans and would be of interest to any researcher that wants to test the sensitivity of population genetic or GWAS analysis to recombination map input.

A common deletion at BAK1 reduces enhancer activity and confers risk of intracranial germ cell tumors

Intracranial germ cell tumors (IGCTs) are rare brain neoplasms that mainly occur in children and adolescents with a particularly high incidence in East Asian populations. Here, we conduct a genome-wide association study (GWAS) of 133 patients with IGCTs and 762 controls of Japanese ancestry. A common 4-bp deletion polymorphism in an enhancer adjacent to BAK1 is significantly associated with the disease risk (rs3831846; P = 2.4 × 10⁻⁹, odds ratio = 2.46 [95% CI: 1.83–3.31], minor allele frequency = 0.43). Rs3831846 is in strong linkage disequilibrium with a testicular GCTs susceptibility variant rs210138. In-vitro reporter assays reveal rs3831846 to be a functional variant attenuating the enhancer activity, suggesting its contribution to IGCTs predisposition through altering BAK1 expression. Risk alleles of testicular GCTs derived from the European GWAS show significant positive correlations in the effect sizes with the Japanese IGCTs GWAS (P = 1.3 × 10⁻⁴, Spearman’s ρ = 0.48). These results suggest the shared genetic susceptibility of GCTs beyond ethnicity and primary sites.

Intracranial germ cell tumors (IGCTs) are rare brain tumors mainly diagnosed in children and young adults. Here, the authors conduct a genome-wide association study for IGCTs, identify a risk locus at BAK1, and characterize its functional consequences.

Demographic history differences between Hispanics and Brazilians imprint haplotype features

Admixture is known to greatly impact the genetic landscape of a population and, while genetic variation underlying human phenotypes has been shown to differ among populations, studies on admixed subjects are still scarce. Latin American populations are the result of complex demographic history, such as 2 or 3-way admixing events, bottlenecks and/or expansions, and adaptive events unique to the American continent. To explore the impact of these events on the genetic structure of Latino populations, we evaluated the following haplotype features: linkage disequilibrium, shared identity by descent segments, runs of homozygosity, and extended haplotype homozygosity (integrated haplotype score) in Latinos represented in the 1000 Genome Project along with array data from 171 Brazilians sampled in the South and Southeast regions of Brazil. We found that linkage disequilibrium decay relates to the amount of American and African ancestry. The extent of identity by descent sharing positively correlates with historical effective population sizes, which we found to be steady or growing, except for Puerto Ricans and Colombians. Long runs of homozygosity, a particular instance of autozygosity, was only enriched in Peruvians and Native Americans. We used simulations to account for random sampling and linkage disequilibrium to filter positive selection indexes and found 244 unique markers under selection, 26 of which are common to 2 or more populations. Some markers exhibiting positive selection signals had estimated time to the most recent common ancestor consistent with human adaptation to the American continent. In conclusion, Latino populations present highly divergent haplotype characteristics that impact genetic architecture and underlie complex phenotypes.

Genomic Insights Into the Demographic History of the Southern Chinese

Southern China is the birthplace of rice-cultivating agriculture and different language families and has also witnessed various human migrations that facilitated cultural diffusions. The fine-scale demographic history in situ that forms present-day local populations, however, remains unclear. To comprehensively cover the genetic diversity in East and Southeast Asia, we generated genome-wide SNP data from 211 present-day Southern Chinese and co-analyzed them with ∼1,200 ancient and modern genomes. In Southern China, language classification is significantly associated with genetic variation but with a different extent of predictability, and there is strong evidence for recent shared genetic history particularly in Hmong–Mien and Austronesian speakers. A geography-related genetic sub-structure that represents the major genetic variation in Southern East Asians is established pre-Holocene and its extremes are represented by Neolithic Fujianese and First Farmers in Mainland Southeast Asia. This sub-structure is largely reduced by admixture in ancient Southern Chinese since &gt; ∼2,000 BP, which forms a “Southern Chinese Cluster” with a high level of genetic homogeneity. Further admixture characterizes the demographic history of the majority of Hmong–Mien speakers and some Kra-Dai speakers in Southwest China happened ∼1,500–1,000 BP, coeval to the reigns of local chiefdoms. In Yellow River Basin, we identify a connection of local populations to genetic sub-structure in Southern China with geographical correspondence appearing &gt; ∼9,000 BP, while the gene flow likely closely related to “Southern Chinese Cluster” since the Longshan period (∼5,000–4,000 BP) forms ancestry profile of Han Chinese Cline.