Genome sequencing elucidates Sardinian genetic architecture and augments association analyses for lipid and blood inflammatory markers

Genetic aspects in schizophrenia. Receptoral theories. Metabolic theories

Ties between schizophrenia (SCZ) and genetics are undeniably significant issue prone to be discussed in the nowadays psychology. Recent research on this domain focuses more on specific genes and heredity (specifically monozygotic pairs of twins) for diagnosing SCZ, than on environmental influences. SCZ is considered a multifactorial disease, thought to convert from a merger of risk and biological genes and environmental factors that could alter and reshape the trajectory of brain development. On this regard, this review sums up recent and innovative methods of distinguishing schizophrenic features from other mental illnesses in patients, based on chromosomal and genes changes. The term “reverse genetics” is no longer up to date, being replaced with “genome scanning” and “positional cloning”. For many researchers, genome scanning is continuing the reverse of the sensible strategy for detecting various important biological disorders, which may start from the discovery of a protein or any other molecule involved in a biological process, being followed by its gene cloning. Genes being discovered in this manner could become candidate genes for the disease. However, genome scanning occurs through testing each chromosomal segment (or mitochondrial genome) for the counter transmission of the disease.

Genetic variant within CDK6 regulates immune response to palbociclib treatment

Everyone carries a set of genetic variants that contribute to regulation of the levels of blood cells, with unknown clinical impact. One of them, rs445 within the cell-cycle checkpoint gene CDK6, reduces the levels of myeloid cell types including granulocytes. We treated CD3+ T cells and whole blood with palbociclib in 41 individuals, who were stratified by genotype for analyses. In T cells we assessed cell cycle and apoptosis, whereas in whole blood, apoptosis in activated (CD11b+), unactivated (CD11b-) granulocytes, cytotoxic (CD8 + CD4-), and helper (CD8-CD4+) T cells. We find that rs445 modulates the immune response of CD8+ T cells. It also increases the level of apoptotic CD11b + activated granulocytes after palbociclib treatment, which, in synergy with neutropenia, may affect drug related adverse events. These results suggest that the effect of palbociclib treatment may depend on underlying genetically encoded individual immune response as well as the direct response to the drug.

Into the Deep: New Data on the Lipid and Fatty Acid Profile of Redfish Sebastes mentella Inhabiting Different Depths in the Irminger Sea

New data on lipid and fatty acid profiles are presented, and the dynamics of the studied components in muscles in the males and females of the beaked redfish, Sebastes mentella, in the depth gradient of the Irminger Sea (North Atlantic) is discussed. The contents of the total lipids (TLs), total phospholipids (PLs), monoacylglycerols (MAGs), diacylglycerols (DAGs), triacylglycerols (TAGs), cholesterol (Chol), Chol esters, non-esterified fatty acids (NEFAs), and wax esters were determined by HPTLC; the phosphatidylserine (PS), phosphatidylethanolamine (PE), phosphatidylinositol (PI), phosphatidylcholine (PC), and lysophosphatidylcholine (LPC) were determined by HPLC; and fatty acids of total lipids were determined using GC. The Chol esters prevailed in muscles over the storage TAGs, and the wax ester content was high, which is a characteristic trait of vertically migrating species. Specific dynamics in certain PL in redfish were found to be depended on depth, suggesting that PLs are involved in the re-arrangement of the membrane physicochemical state and the maintenance of motor activity under high hydrostatic pressure. The high contents of DHA and EPA were observed in beaked redfish muscles is the species' characteristic trait. The MUFAs in muscles include dietary markers of zooplankton (copepods)-20:1(n-9) and 22:1(n-11), whose content was found to be lower in fish sampled from greater depths.

Genome-wide search for genes affecting the age at diagnosis of type 1 diabetes

BACKGROUND

Type 1 diabetes (T1D) is an autoimmune disease affecting individuals in the early years of life. Although previous studies have identified genetic loci influencing T1D diagnosis age, these studies did not investigate the genome with high resolution.

OBJECTIVE AND METHODS

We performed a genome-wide meta-analysis for age at diagnosis with cohorts from Finland (Finnish Diabetic Nephropathy Study), the United Kingdom (UK Genetic Resource Investigating Diabetes) and Sardinia. Through SNP associations, transcriptome-wide association analysis linked T1D diagnosis age and gene expression.

RESULTS

We identified two chromosomal regions associated with T1D diagnosis age: multiple independent variants in the HLA region on chromosome 6 and a locus on chromosome 17q12. We performed gene-level association tests with transcriptome prediction models from two whole blood datasets, lymphocyte cell line, spleen, pancreas and small intestine tissues. Of the non-HLA genes, lower PNMT expression in whole blood, and higher IKZF3 and ZPBP2, and lower ORMDL3 and GSDMB transcription levels in multiple tissues were associated with lower T1D diagnosis age (FDR = 0.05). These genes lie on chr17q12 which is associated with T1D, other autoimmune diseases, and childhood asthma. Additionally, higher expression of PHF20L1, a gene not previously implicated in T1D, was associated with lower diagnosis age in lymphocytes, pancreas, and spleen. Altogether, the non-HLA associations were enriched in open chromatin in various blood cells, blood vessel tissues and foetal thymus tissue.

CONCLUSION

Multiple genes on chr17q12 and PHF20L1 on chr8 were associated with T1D diagnosis age and only further studies may elucidate the role of these genes for immunity and T1D onset.

Genome-wide association study of circulating interleukin 6 levels identifies novel loci

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Geographic variation in the polygenic score of height in Japan

A geographical gradient of height has existed in Japan for approximately 100 years. People in northern Japan tend to be taller than those in southern Japan. The differences in annual temperature and day length between the northern and southern prefectures of Japan have been suggested as possible causes of the height gradient. Although height is well known to be a polygenic trait with high heritability, the genetic contributions to the gradient have not yet been explored. Polygenic score (PS) is calculated by aggregating the effects of genetic variants identified by genome-wide association studies (GWASs) to predict the traits of individual subjects. Here, we calculated the PS of height for 10,840 Japanese individuals from all 47 prefectures in Japan. The median height PS for each prefecture was significantly correlated with the mean height of females and males obtained from another independent Japanese nation-wide height dataset, suggesting genetic contribution to the observed height gradient. We also found that individuals and prefectures genetically closer to continental East Asian ancestry tended to have a higher PS; modern Japanese people are considered to have originated as result of admixture between indigenous Jomon people and immigrants from continental East Asia. Another PS analysis based on the GWAS using only the mainland Japanese was conducted to evaluate the effect of population stratification on PS. The result also supported genetic contribution to height, and indicated that the PS might be affected by a bias due to population stratification even in a relatively homogenous population like Japanese.

The forensic landscape and the population genetic analyses of Hainan Li based on massively parallel sequencing DNA profiling

Due to the formation of the Qiongzhou Strait by climate change and marine transition, Hainan island was isolated from the mainland southern China during the Last Glacial Maximum. Hainan island, located at the southernmost part of China and separated from the Leizhou Peninsula by the Qiongzhou Strait, laid on one of the modern human northward migration routes from Southeast Asia to East Asia. The Hlai language-speaking Li minority, the second largest population after Han Chinese in Hainan island, is the direct descendants of the initial migrants in Hainan island and has unique ethnic properties and derived characteristics; however, the forensic-associated studies on Hainan Li population are still insufficient. Hence, 136 Hainan Li individuals were genotyped in this study using the MPS-based ForenSeq™ DNA Signature Prep Kit (DNA Primer Set A, DPMA) to characterize the forensic genetic polymorphism landscape, and DNA profiles were obtained from 152 different molecular genetic markers (27 autosomal STRs, 24 Y-STRs, 7 X-STRs, and 94 iiSNPs). A total of 419 distinct length variants and 586 repeat sequence sub-variants, with 31 novel alleles (at 17 loci), were identified across the 58 STR loci from the DNA profiles of Hainan Li population. We evaluated the forensic characteristics and efficiencies of DPMA, demonstrating that the STRs and iiSNPs in DPMA were highly polymorphic in Hainan Li population and could be employed in forensic applications. In addition, we set up three datasets, which included the genetic data of (i) iiSNPs (27 populations, 2640 individuals), (ii) Y-STRs (42 populations, 8281 individuals), and (iii) Y haplogroups (123 populations, 4837 individuals) along with the population ancestries and language families, to perform population genetic analyses separately from different perspectives. In conclusion, the phylogenetic analyses indicated that Hainan Li, with a southern East Asia origin and Tai-Kadai language-speaking language, is an isolated population relatively. But the genetic pool of Hainan Li influenced by the limited gene flows from other Tai-Kadai populations and Hainan populations. Furthermore, the establishment of isolated population models will be beneficial to clarify the exquisite population structures and develop specific genetic markers for subpopulations in forensic genetic fields.

PRF1 mutation alters immune system activation, inflammation, and risk of autoimmunity

BACKGROUND

Defective alleles within the PRF1 gene, encoding the pore-forming protein perforin, in combination with environmental factors, cause familial type 2 hemophagocytic lymphohistiocytosis (FHL2), a rare, severe autosomal recessive childhood disorder characterized by massive release of cytokines-cytokine storm.

OBJECTIVE

The aim of this study was to determine the function of hypomorph PRF1:p.A91V g.72360387 G > A on multiple sclerosis (MS) and type 1 diabetes (T1D).

METHODS

We cross-compare the association data for PRF1:p.A91V mutation derived from GWAS on adult MS and pediatric T1D in Sardinians. The novel association with T1D was replicated in metanalysis in 12,584 cases and 17,692 controls from Sardinia, the United Kingdom, and Scotland. To dissect this mutation function, we searched through the coincident association immunophenotypes in additional set of general population Sardinians.

RESULTS

We report that PRF1:p.A91V, is associated with increase of lymphocyte levels, especially within the cytotoxic memory T-cells, at general population level with reduced interleukin 7 receptor expression on these cells. The minor allele increased risk of MS, in 2903 cases and 2880 controls from Sardinia p = 2.06 × 10^-4, odds ratio OR = 1.29, replicating a previous finding, whereas it protects from T1D p = 1.04 × 10^-5, OR = 0.82.

CONCLUSION

Our results indicate opposing contributions of the cytotoxic T-cell compartment to MS and T1D pathogenesis.

Functional and clinical implications of genetic structure in 1686 Italian exomes

To reconstruct the phenotypical and clinical implications of the Italian genetic structure, we thoroughly analyzed a whole-exome sequencing data set comprised of 1686 healthy Italian individuals. We found six previously unreported variants with remarkable frequency differences between Northern and Southern Italy in the HERC2, OR52R1, ADH1B, and THBS4 genes. We reported 36 clinically relevant variants (submitted as pathogenic, risk factors, or drug response in ClinVar) with significant frequency differences between Italy and Europe. We then explored putatively pathogenic variants in the Italian exome. On average, our Italian individuals carried 16.6 protein-truncating variants (PTVs), with 2.5% of the population having a PTV in one of the 59 American College of Medical Genetics (ACMG) actionable genes. Lastly, we looked for PTVs that are likely to cause Mendelian diseases. We found four heterozygous PTVs in haploinsufficient genes (KAT6A, PTCH1, and STXBP1) and three homozygous PTVs in genes causing recessive diseases (DPYD, FLG, and PYGM). Comparing frequencies from our data set to other public databases, like gnomAD, we showed the importance of population-specific databases for a more accurate assessment of variant pathogenicity. For this reason, we made aggregated frequencies from our data set publicly available as a tool for both clinicians and researchers (http://nigdb.cineca.it; NIG-ExIT).

Population structure of indigenous inhabitants of Arabia

Modern day Saudi Arabia occupies the majority of historical Arabia, which may have contributed to ancient waves of migration out of Africa. This ancient history has left a lasting imprint in the genetics of the region, including the diverse set of tribes that call Saudi Arabia their home. How these tribes relate to each other and to the world's major populations remains an unanswered question. In an attempt to improve our understanding of the population structure of Saudi Arabia, we conducted genomic profiling of 957 unrelated individuals who self-identify with 28 large tribes in Saudi Arabia. Consistent with the tradition of intra-tribal unions, the subjects showed strong clustering along tribal lines with the distance between clusters correlating with their geographical proximities in Arabia. However, these individuals form a unique cluster when compared to the world's major populations. The ancient origin of these tribal affiliations is supported by analyses that revealed little evidence of ancestral origin from within the 28 tribes. Our results disclose a granular map of population structure and have important implications for future genetic studies into Mendelian and common diseases in the region.

Loss-of-function genomic variants highlight potential therapeutic targets for cardiovascular disease

Pharmaceutical drugs targeting dyslipidemia and cardiovascular disease (CVD) may increase the risk of fatty liver disease and other metabolic disorders. To identify potential novel CVD drug targets without these adverse effects, we perform genome-wide analyses of participants in the HUNT Study in Norway (n = 69,479) to search for protein-altering variants with beneficial impact on quantitative blood traits related to cardiovascular disease, but without detrimental impact on liver function. We identify 76 (11 previously unreported) presumed causal protein-altering variants associated with one or more CVD- or liver-related blood traits. Nine of the variants are predicted to result in loss-of-function of the protein. This includes ZNF529:p.K405X, which is associated with decreased low-density-lipoprotein (LDL) cholesterol (P = 1.3 × 10-8) without being associated with liver enzymes or non-fasting blood glucose. Silencing of ZNF529 in human hepatoma cells results in upregulation of LDL receptor and increased LDL uptake in the cells. This suggests that inhibition of ZNF529 or its gene product should be prioritized as a novel candidate drug target for treating dyslipidemia and associated CVD.

Loss-of-function genomic variants highlight potential therapeutic targets for cardiovascular disease

Pharmaceutical drugs targeting dyslipidemia and cardiovascular disease (CVD) may increase the risk of fatty liver disease and other metabolic disorders. To identify potential novel CVD drug targets without these adverse effects, we perform genome-wide analyses of participants in the HUNT Study in Norway (n = 69,479) to search for protein-altering variants with beneficial impact on quantitative blood traits related to cardiovascular disease, but without detrimental impact on liver function. We identify 76 (11 previously unreported) presumed causal protein-altering variants associated with one or more CVD- or liver-related blood traits. Nine of the variants are predicted to result in loss-of-function of the protein. This includes ZNF529:p.K405X, which is associated with decreased low-density-lipoprotein (LDL) cholesterol (P = 1.3 × 10^-8) without being associated with liver enzymes or non-fasting blood glucose. Silencing of ZNF529 in human hepatoma cells results in upregulation of LDL receptor and increased LDL uptake in the cells. This suggests that inhibition of ZNF529 or its gene product should be prioritized as a novel candidate drug target for treating dyslipidemia and associated CVD.

The tip of the iceberg for diagnostic dilemmas: Performance of current diagnostics and future complementary screening approaches

Genetic testing is currently the leading edge of clinical care when it comes to diagnostics. However, many questions remain unanswered even when employing next-generation sequencing techniques due to our inability to decode genetic variations and our limited repertoire of available diagnoses. Accordingly, diagnostic yields for current genomic screenings are <50% and fail to provide the whole picture, leaving the remaining patients without a definitive diagnosis. Human phenotypic/disease expression is explained by alterations not only at the genome, but also at the transcriptome, proteome and metabolome levels. These "higher" complexity levels represent at wealth of information, and diagnostic screenings tests at these levels have been shown to significantly improve diagnostic yields in specific populations compared to conventional diagnostic workup or gold standards in use (7-30% increase in diagnostic yields, depending on the population, approach and gold standard being compared against). However, these are not yet routinely available to clinicians. Due to their dynamic and modifiable nature, tapping into data from different omics will improve our understanding of the pathophysiological bases underlying (many yet to characterize) human disorders. We herein review how alterations at these levels (e.g. post-transcriptional and post-translational) may be pathogenic, how such tests may be implemented and in which situations they are of significant utility.

Complex genetic signatures in immune cells underlie autoimmunity and inform therapy

We report on the influence of ~22 million variants on 731 immune cell traits in a cohort of 3,757 Sardinians. We detected 122 significant (P < 1.28 × 10^-11) independent association signals for 459 cell traits at 70 loci (53 of them novel) identifying several molecules and mechanisms involved in cell regulation. Furthermore, 53 signals at 36 loci overlapped with previously reported disease-associated signals, predominantly for autoimmune disorders, highlighting intermediate phenotypes in pathogenesis. Collectively, our findings illustrate complex genetic regulation of immune cells with highly selective effects on autoimmune disease risk at the cell-subtype level. These results identify drug-targetable pathways informing the design of more specific treatments for autoimmune diseases.

Whole genome sequencing of elite athletes

Whole genome sequencing (WGS) has great potential to explore all possible DNA variants associated with physical performance, psychological traits and health conditions of athletes. Here we present, for the first time, annotation of genomic variants of elite athletes, based on the WGS of 20 Tatar male wrestlers. The maximum number of high-quality variants per sample was over 3.8 M for single nucleotide polymorphisms (SNPs) and about 0.64 M for indels. The maximum number of nonsense mutations was 148 single nucleotide variants (SNVs) per individual. Athletes' genomes on average contained 18.9 nonsense SNPs in a homozygous state per sample, while non-athletes' exomes (Tatar controls, n = 19) contained 18 nonsense SNPs. Finally, we applied genomic data for the association analysis and used reaction time (RT) as an example. Out of 1884 known genome-wide significant SNPs related to RT, we identified four SNPs (KIF27 rs10125715, APC rs518013, TMEM229A rs7783359, LRRN3 rs80054135) associated with RT in wrestlers. The cumulative number of favourable alleles (KIF27 A, APC A, TMEM229A T, LRRN3 T) was significantly correlated with RT both in wrestlers (P = 0.0003) and an independent cohort (n = 43) of physically active subjects (P = 0.029). Furthermore, we found that the frequencies of the APC A (53.3 vs 44.0%, P = 0.033) and LRRN3 T (7.5 vs 2.8%, P = 0.009) alleles were significantly higher in elite athletes (n = 107) involved in sports with RT as an essential component of performance (combat sports, table tennis and volleyball) compared to less successful (n = 176) athletes. The LRRN3 T allele was also over-represented in elite athletes (7.5%) in comparison with 189 controls (2.9%, P = 0.009). In conclusion, we present the first WGS study of athletes showing that WGS can be applied in sport and exercise science.

Sequencing and imputation in GWAS: Cost-effective strategies to increase power and genomic coverage across diverse populations

A key aim for current genome-wide association studies (GWAS) is to interrogate the full spectrum of genetic variation underlying human traits, including rare variants, across populations. Deep whole-genome sequencing is the gold standard to fully capture genetic variation, but remains prohibitively expensive for large sample sizes. Array genotyping interrogates a sparser set of variants, which can be used as a scaffold for genotype imputation to capture a wider set of variants. However, imputation quality depends crucially on reference panel size and genetic distance from the target population. Here, we consider sequencing a subset of GWAS participants and imputing the rest using a reference panel that includes both sequenced GWAS participants and an external reference panel. We investigate how imputation quality and GWAS power are affected by the number of participants sequenced for admixed populations (African and Latino Americans) and European population isolates (Sardinians and Finns), and identify powerful, cost-effective GWAS designs given current sequencing and array costs. For populations that are well-represented in existing reference panels, we find that array genotyping alone is cost-effective and well-powered to detect common- and rare-variant associations. For poorly represented populations, sequencing a subset of participants is often most cost-effective, and can substantially increase imputation quality and GWAS power.

Characterization of Exome Variants and Their Metabolic Impact in 6,716 American Indians from the Southwest US

Applying exome sequencing to populations with unique genetic architecture has the potential to reveal novel genes and variants associated with traits and diseases. We sequenced and analyzed the exomes of 6,716 individuals from a Southwestern American Indian (SWAI) population with well-characterized metabolic traits. We found that the SWAI population has distinct allelic architecture compared to populations of European and East Asian ancestry, and there were many predicted loss-of-function (pLOF) and nonsynonymous variants that were highly enriched or private in the SWAI population. We used pLOF and nonsynonymous variants in the SWAI population to evaluate gene-burden associations of candidate genes from European genome-wide association studies (GWASs) for type 2 diabetes, body mass index, and four major plasma lipids. We found 19 significant gene-burden associations for 11 genes, providing additional evidence for prioritizing candidate effector genes of GWAS signals. Interestingly, these associations were mainly driven by pLOF and nonsynonymous variants that are unique or highly enriched in the SWAI population. Particularly, we found four pLOF or nonsynonymous variants in APOB, APOE, PCSK9, and TM6SF2 that are private or enriched in the SWAI population and associated with low-density lipoprotein (LDL) cholesterol levels. Their large estimated effects on LDL cholesterol levels suggest strong impacts on protein function and potential clinical implications of these variants in cardiovascular health. In summary, our study illustrates the utility and potential of exome sequencing in genetically unique populations, such as the SWAI population, to prioritize candidate effector genes within GWAS loci and to find additional variants in known disease genes with potential clinical impact.

Evidence of Polygenic Adaptation in Sardinia at Height-Associated Loci Ascertained from the Biobank Japan

Adult height is one of the earliest putative examples of polygenic adaptation in humans. However, this conclusion was recently challenged because residual uncorrected stratification from large-scale consortium studies was considered responsible for the previously noted genetic difference. It thus remains an open question whether height loci exhibit signals of polygenic adaptation in any human population. We re-examined this question, focusing on one of the shortest European populations, the Sardinians, in addition to mainland European populations. We utilized height-associated loci from the Biobank Japan (BBJ) dataset to further alleviate concerns of biased ascertainment of GWAS loci and showed that the Sardinians remain significantly shorter than expected under neutrality (∼0.22 standard deviation shorter than Utah residents with ancestry from northern and western Europe [CEU] on the basis of polygenic height scores, p = 3.89 × 10-4). We also found the trajectory of polygenic height scores between the Sardinian and the British populations diverged over at least the last 10,000 years (p = 0.0082), consistent with a signature of polygenic adaptation driven primarily by the Sardinian population. Although the polygenic score-based analysis showed a much subtler signature in mainland European populations, we found a clear and robust adaptive signature in the UK population by using a haplotype-based statistic, the trait singleton density score (tSDS), driven by the height-increasing alleles (p = 9.1 × 10-4). In summary, by ascertaining height loci in a distant East Asian population, we further supported the evidence of polygenic adaptation at height-associated loci among the Sardinians. In mainland Europeans, the adaptive signature was detected in haplotype-based analysis but not in polygenic score-based analysis.

Genetics of Atrial Fibrillation in 2020: GWAS, Genome Sequencing, Polygenic Risk, and Beyond

Atrial fibrillation is a common heart rhythm disorder that leads to an increased risk for stroke and heart failure. Atrial fibrillation is a complex disease with both environmental and genetic risk factors that contribute to the arrhythmia. Over the last decade, rapid progress has been made in identifying the genetic basis for this common condition. In this review, we provide an overview of the primary types of genetic analyses performed for atrial fibrillation, including linkage studies, genome-wide association studies, and studies of rare coding variation. With these results in mind, we aim to highlighting the existing knowledge gaps and future directions for atrial fibrillation genetics research.

A population-specific low-frequency variant of SLC22A12 (p.W258*) explains nearby genome-wide association signals for serum uric acid concentrations among Koreans

Prolonged hyperuricemia is a cause of gout and an independent risk factor for chronic health conditions including diabetes and chronic kidney diseases. Genome-wide association studies (GWASs) for serum uric acid (SUA) concentrations have repeatedly confirmed genetic loci including those encoding uric acid transporters such as ABCG2 and SLC9A2. However, many single nucleotide polymorphisms (SNPs) found in GWASs have been common variants with small effects and unknown functions. In addition, there is still much heritability to be explained. To identify the causative genetic variants for SUA concentrations in Korean subjects, we conducted a GWAS (1902 males) and validation study (2912 males and females) and found four genetic loci reaching genome-wide significance on chromosomes 4 (ABCG2) and 11 (FRMD8, EIF1AD and SLC22A12-NRXN2). Three loci on chromosome 11 were distributed within a distance of 1.3 megabases and they were in weak or moderate linkage disequilibrium (LD) states (r² = 0.02–0.68). In a subsequent association analysis on the GWAS loci of chromosome 11 using closely positioned markers derived from whole genome sequencing data, the most significant variant to be linked with the nearby GWAS signal was rs121907892 (c.774G>A, p.W258*) of the SLC22A12 gene. This variant, and each of the three GWAS SNPs, were in LD (r² = 0.06–0.32). The strength of association of SNPs with SUA concentration (negative logarithm of P-values) and the genetic distance (r² of LD) between rs121907892 and the surrounding SNPs showed a quantitative correlation. This variant has been found only in Korean and Japanese subjects and is known to lower the SUA concentration in the general population. Thus, this low-frequency variant, rs121907892, known to regulate SUA concentrations in previous studies, is responsible for the nearby GWAS signals.

MTHFR, XRCC1 and OGG1 genetic polymorphisms in breast cancer: a case-control study in a population from North Sardinia

BACKGROUND

Despite conflicting results, considerable evidence suggests the association between single nucleotide polymorphisms in MTHFR, XRCC1 and OGG1 genes and, risk of developing breast cancer. Here a case-control study is reported, including 135 breat cancer patients and 112 healthy women, all representative of Northern Sardinian population.

METHODS

Polymerase chain reaction/restriction fragment length polymorphism method was used to determine the genotypes of five polymorphisms: MTHFR C677T (rs1801133) and A1298C (rs1801131), XRCC1 Arg194Trp (rs1799782) and Arg399Gln (rs25487) and OGG1 Ser326Cys (rs1052133). Allelic, genotypic and haplotype association analyses with disease risk and clinicopathological parameters were performed.

RESULTS

A nominally significant association with breast cancer risk was observed for MTHFR C677T polymorphism heterozygous genotype in the codominant model (OR: 0.57, 95% CI: 0.32-1.00, p = 0.049) and for Cys/Cys genotype of the OGG1 Ser326Cys polymorphism in the recessive model (OR: 0.23, 95% CI: 0.05-1.11, p = 0.0465). No significant differences were found at genotype-level for A1298C polymorphism of the MTHFR gene and Arg194Trp and Arg399Gln of the XRCC1 gene. Furthermore, the OGG1 and XRCC1 rs25487 polymorphisms were nominally associated with PgR, Her2 status and with sporadic breast cancer, respectively.

CONCLUSIONS

Based on genetic characteristics of individuals included in this study, results suggest that MTHFR CT and OGG1 Cys/Cys genotypes have a protective effect that may have an influence on breast cancer risk in a representative Northern Sardinian population.

Genetic history from the Middle Neolithic to present on the Mediterranean island of Sardinia

The island of Sardinia has been of particular interest to geneticists for decades. The current model for Sardinia's genetic history describes the island as harboring a founder population that was established largely from the Neolithic peoples of southern Europe and remained isolated from later Bronze Age expansions on the mainland. To evaluate this model, we generate genome-wide ancient DNA data for 70 individuals from 21 Sardinian archaeological sites spanning the Middle Neolithic through the Medieval period. The earliest individuals show a strong affinity to western Mediterranean Neolithic populations, followed by an extended period of genetic continuity on the island through the Nuragic period (second millennium BCE). Beginning with individuals from Phoenician/Punic sites (first millennium BCE), we observe spatially-varying signals of admixture with sources principally from the eastern and northern Mediterranean. Overall, our analysis sheds light on the genetic history of Sardinia, revealing how relationships to mainland populations shifted over time.

A polygenic predictor of treatment-resistant depression using whole exome sequencing and genome-wide genotyping

Treatment-resistant depression (TRD) occurs in ~30% of patients with major depressive disorder (MDD) but the genetics of TRD was previously poorly investigated. Whole exome sequencing and genome-wide genotyping were available in 1209 MDD patients after quality control. Antidepressant response was compared to non-response to one treatment and non-response to two or more treatments (TRD). Differences in the risk of carrying damaging variants were tested. A score expressing the burden of variants in genes and pathways was calculated weighting each variant for its functional (Eigen) score and frequency. Gene-based and pathway-based scores were used to develop predictive models of TRD and non-response using gradient boosting in 70% of the sample (training) which were tested in the remaining 30% (testing), evaluating also the addition of clinical predictors. Independent replication was tested in STAR*D and GENDEP using exome array-based data. TRD and non-responders did not show higher risk to carry damaging variants compared to responders. Genes/pathways associated with TRD included those modulating cell survival and proliferation, neurodegeneration, and immune response. Genetic models showed significant prediction of TRD vs. response and they were improved by the addition of clinical predictors, but they were not significantly better than clinical predictors alone. Replication results were driven by clinical factors, except for a model developed in subjects treated with serotonergic antidepressants, which showed a clear improvement in prediction at the extremes of the genetic score distribution in STAR*D. These results suggested relevant biological mechanisms implicated in TRD and a new methodological approach to the prediction of TRD.

The stabilization of equilibria in evolutionary game dynamics through mutation: mutation limits in evolutionary games

The multi-population replicator dynamics is a dynamic approach to coevolving populations and multi-player games and is related to Cross learning. In general, not every equilibrium is a Nash equilibrium of the underlying game, and the convergence is not guaranteed. In particular, no interior equilibrium can be asymptotically stable in the multi-population replicator dynamics, e.g. resulting in cyclic orbits around a single interior Nash equilibrium. We introduce a new notion of equilibria of replicator dynamics, called mutation limits, based on a naturally arising, simple form of mutation, which is invariant under the specific choice of mutation parameters. We prove the existence of mutation limits for a large class of games, and consider a particularly interesting subclass called attracting mutation limits. Attracting mutation limits are approximated in every (mutation-)perturbed replicator dynamics, hence they offer an approximate dynamic solution to the underlying game even if the original dynamic is not convergent. Thus, mutation stabilizes the system in certain cases and makes attracting mutation limits near attainable. Hence, attracting mutation limits are relevant as a dynamic solution concept of games. We observe that they have some similarity to Q-learning in multi-agent reinforcement learning. Attracting mutation limits do not exist in all games, however, raising the question of their characterization.

Critical assessment of bioinformatics methods for the characterization of pathological repeat expansions with single-molecule sequencing data

A number of studies have reported the successful application of single-molecule sequencing technologies to the determination of the size and sequence of pathological expanded microsatellite repeats over the last 5 years. However, different custom bioinformatics pipelines were employed in each study, preventing meaningful comparisons and somewhat limiting the reproducibility of the results. In this review, we provide a brief summary of state-of-the-art methods for the characterization of expanded repeats alleles, along with a detailed comparison of bioinformatics tools for the determination of repeat length and sequence, using both real and simulated data. Our reanalysis of publicly available human genome sequencing data suggests a modest, but statistically significant, increase of the error rate of single-molecule sequencing technologies at genomic regions containing short tandem repeats. However, we observe that all the methods herein tested, irrespective of the strategy used for the analysis of the data (either based on the alignment or assembly of the reads), show high levels of sensitivity in both the detection of expanded tandem repeats and the estimation of the expansion size, suggesting that approaches based on single-molecule sequencing technologies are highly effective for the detection and quantification of tandem repeat expansions and contractions.

A bird's-eye view of Italian genomic variation through whole-genome sequencing

The genomic variation of the Italian peninsula populations is currently under characterised: the only Italian whole-genome reference is represented by the Tuscans from the 1000 Genome Project. To address this issue, we sequenced a total of 947 Italian samples from three different geographical areas. First, we defined a new Italian Genome Reference Panel (IGRP1.0) for imputation, which improved imputation accuracy, especially for rare variants, and we tested it by GWAS analysis on red blood traits. Furthermore, we extended the catalogue of genetic variation investigating the level of population structure, the pattern of natural selection, the distribution of deleterious variants and occurrence of human knockouts (HKOs). Overall the results demonstrate a high level of genomic differentiation between cohorts, different signatures of natural selection and a distinctive distribution of deleterious variants and HKOs, confirming the necessity of distinct genome references for the Italian population.

Improved power and precision with whole genome sequencing data in genome-wide association studies of inflammatory biomarkers

Genome-wide association studies (GWAS) have identified associations between thousands of common genetic variants and human traits. However, common variants usually explain a limited fraction of the heritability of a trait. A powerful resource for identifying trait-associated variants is whole genome sequencing (WGS) data in cohorts comprised of families or individuals from a limited geographical area. To evaluate the power of WGS compared to imputations, we performed GWAS on WGS data for 72 inflammatory biomarkers, in a kinship-structured cohort. When using WGS data, we identified 18 novel associations that were not detected when analyzing the same biomarkers with genotyped or imputed SNPs. Five of the novel top variants were low frequency variants with a minor allele frequency (MAF) of <5%. Our results suggest that, even when applying a GWAS approach, we gain power and precision using WGS data, presumably due to more accurate determination of genotypes. The lack of a comparable dataset for replication of our results is a limitation in our study. However, this further highlights that there is a need for more genetic epidemiological studies based on WGS data.

Kinpute: using identity by descent to improve genotype imputation

MOTIVATION

Genotype imputation, though generally accurate, often results in many genotypes being poorly imputed, particularly in studies where the individuals are not well represented by standard reference panels. When individuals in the study share regions of the genome identical by descent (IBD), it is possible to use this information in combination with a study-specific reference panel (SSRP) to improve the imputation results. Kinpute uses IBD information-due to recent, familial relatedness or distant, unknown ancestors-in conjunction with the output from linkage disequilibrium (LD) based imputation methods to compute more accurate genotype probabilities. Kinpute uses a novel method for IBD imputation, which works even in the absence of a pedigree, and results in substantially improved imputation quality.

RESULTS

Given initial estimates of average IBD between subjects in the study sample, Kinpute uses a novel algorithm to select an optimal set of individuals to sequence and use as an SSRP. Kinpute is designed to use as input both this SSRP and the genotype probabilities output from other LD-based imputation software, and uses a new method to combine the LD imputed genotype probabilities with IBD configurations to substantially improve imputation. We tested Kinpute on a human population isolate where 98 individuals have been sequenced. In half of this sample, whose sequence data was masked, we used Impute2 to perform LD-based imputation and Kinpute was used to obtain higher accuracy genotype probabilities. Measures of imputation accuracy improved significantly, particularly for those genotypes that Impute2 imputed with low certainty.

AVAILABILITY AND IMPLEMENTATION

Kinpute is an open-source and freely available C++ software package that can be downloaded from.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

NARD: whole-genome reference panel of 1779 Northeast Asians improves imputation accuracy of rare and low-frequency variants

Here, we present the Northeast Asian Reference Database (NARD), including whole-genome sequencing data of 1779 individuals from Korea, Mongolia, Japan, China, and Hong Kong. NARD provides the genetic diversity of Korean (n = 850) and Mongolian (n = 384) ancestries that were not present in the 1000 Genomes Project Phase 3 (1KGP3). We combined and re-phased the genotypes from NARD and 1KGP3 to construct a union set of haplotypes. This approach established a robust imputation reference panel for Northeast Asians, which yields the greatest imputation accuracy of rare and low-frequency variants compared with the existing panels. NARD imputation panel is available at https://nard.macrogen.com/ .

Meta-MultiSKAT: Multiple phenotype meta-analysis for region-based association test

The power of genetic association analyses can be increased by jointly meta-analyzing multiple correlated phenotypes. Here, we develop a meta-analysis framework, Meta-MultiSKAT, that uses summary statistics to test for association between multiple continuous phenotypes and variants in a region of interest. Our approach models the heterogeneity of effects between studies through a kernel matrix and performs a variance component test for association. Using a genotype kernel, our approach can test for rare-variants and the combined effects of both common and rare-variants. To achieve robust power, within Meta-MultiSKAT, we developed fast and accurate omnibus tests combining different models of genetic effects, functional genomic annotations, multiple correlated phenotypes, and heterogeneity across studies. In addition, Meta-MultiSKAT accommodates situations where studies do not share exactly the same set of phenotypes or have differing correlation patterns among the phenotypes. Simulation studies confirm that Meta-MultiSKAT can maintain the type-I error rate at the exome-wide level of 2.5 × 10^-6 . Further simulations under different models of association show that Meta-MultiSKAT can improve the power of detection from 23% to 38% on average over single phenotype-based meta-analysis approaches. We demonstrate the utility and improved power of Meta-MultiSKAT in the meta-analyses of four white blood cell subtype traits from the Michigan Genomics Initiative (MGI) and SardiNIA studies.

A nonparametric test for association with multiple loci in the retrospective case-control study

The genome-wide association studies aim at identifying common or rare variants associated with common diseases and explaining more heritability. It is well known that common diseases are influenced by multiple single nucleotide polymorphisms (SNPs) that are usually correlated in location or function. In order to powerfully detect association signals, it is highly desirable to take account of correlations or linkage disequilibrium (LD) information among multiple SNPs in testing for association. In this article, we propose a test SLIDE that depicts the difference of the average multi-locus genotypes between cases and controls and derive its variance-covariance matrix in the retrospective design. This matrix is composed of the pairwise LD between SNPs. Thus SLIDE can borrow the strength from an external database in the population of interest with a few thousands to hundreds of thousands individuals to improve the power for detecting association. Extensive simulations show that SLIDE has apparent superiority over the existing methods, especially in the situation involving both common and rare variants, both protective and deleterious variants. Furthermore, the efficiency of the proposed method is demonstrated in the application to the data from the Wellcome Trust Case Control Consortium.

Interrogation of human hematopoiesis at single-cell and single-variant resolution

Widespread linkage disequilibrium and incomplete annotation of cell-to-cell state variation represent substantial challenges to elucidating mechanisms of trait-associated genetic variation. Here, we perform genetic fine-mapping for blood cell traits in the UK Biobank to identify putative causal variants. These variants are enriched in genes encoding for proteins in trait-relevant biological pathways and in accessible chromatin of hematopoietic progenitors. For regulatory variants, we explore patterns of developmental enhancer activity, predict molecular mechanisms, and identify likely target genes. In several instances, we localize multiple independent variants to the same regulatory element or gene. We further observe that variants with pleiotropic effects preferentially act in common progenitor populations to direct the production of distinct lineages. Finally, we leverage fine-mapped variants in conjunction with continuous epigenomic annotations to identify trait-cell type enrichments within closely related populations and in single cells. Our study provides a comprehensive framework for single-variant and single-cell analyses of genetic associations.

Fine mapping of blood cell traits in UK Biobank identifies putative causal variants and enrichment of fine-mapped variants in accessible chromatin of hematopoietic progenitor cells. The study provides an analytical framework for single-variant and single-cell analyses of genetic associations.

Prevalence of multiple sclerosis in Sardinia: A systematic cross-sectional multi-source survey

Background:

Partial surveys in sub-regions of Sardinia have suggested a high prevalence of multiple sclerosis (MS) on the island, relative to other Mediterranean populations. We assessed the island-wide prevalence of MS and its detailed distribution in Sardinia.

Methods:

The study population consisted of 5677 MS patients, 1735 men and 3942 women, living in Sardinia. Neurologists retrospectively examined electronic and paper-based records of patients with a diagnosis of MS. The data were then linked to the administrative health information systems. Crude, age-, and sex-specific prevalence estimates of disease were calculated.

Results:

The overall age-adjusted MS prevalence was 330 per 100,000 (95% confidence interval (CI) 321–338) in individuals older than 15 years, 447 in women (95% CI 433–461), and 205 in men (95% CI 195–214). The prevalence was highest in the Ogliastra and Nuoro districts, respectively, 425 (95% CI 372–478) and 419 (95% CI 387–451), and lowest in the Olbia-Tempio district, 217 (95% CI 195–239). Most cases had relapsing–remitting MS (79.3%), 16.3% had secondary-progressive MS, and 4.4% had primary-progressive MS.

Conclusion:

These prevalence are among the highest reported so far worldwide. They provide estimates for comparative analyses in other populations and are essential for public health interventions.

Parent-of-origin effects on quantitative phenotypes in a large Hutterite pedigree

The impact of the parental origin of associated alleles in GWAS has been largely ignored. Yet sequence variants could affect traits differently depending on whether they are inherited from the mother or the father, as in imprinted regions, where identical inherited DNA sequences can have different effects based on the parental origin. To explore parent-of-origin effects (POEs), we studied 21 quantitative phenotypes in a large Hutterite pedigree to identify variants with single parent (maternal-only or paternal-only) effects, and then variants with opposite parental effects. Here we show that POEs, which can be opposite in direction, are relatively common in humans, have potentially important clinical effects, and will be missed in traditional GWAS. We identified POEs with 11 phenotypes, most of which are risk factors for cardiovascular disease. Many of the loci identified are characteristic of imprinted regions and are associated with the expression of nearby genes.

A scientometric review of genome-wide association studies

This scientometric review of genome-wide association studies (GWAS) from 2005 to 2018 (3639 studies; 3508 traits) reveals extraordinary increases in sample sizes, rates of discovery and traits studied. A longitudinal examination shows fluctuating ancestral diversity, still predominantly European Ancestry (88% in 2017) with 72% of discoveries from participants recruited from three countries (US, UK, Iceland). US agencies, primarily NIH, fund 85% and women are less often senior authors. We generate a unique GWAS H-Index and reveal a tight social network of prominent authors and frequently used data sets. We conclude with 10 evidence-based policy recommendations for scientists, research bodies, funders, and editors.

Relative impact of indels versus SNPs on complex disease

It is unclear whether insertions and deletions (indels) are more likely to influence complex traits than abundant single‐nucleotide polymorphisms (SNPs). We sought to understand which category of variation is more likely to impact health. Using the SardiNIA study as an exemplar, we characterized 478,876 common indels and 8,246,244 common SNPs in up to 5,949 well‐phenotyped individuals from an isolated valley in Sardinia. We assessed association between 120 traits, resulting in 89 nonoverlapping‐associated loci.We evaluated whether indels were enriched among credible sets of potential causal variants. These credible sets included 1,319 SNPs and 88 indels. We did not find indels to be significantly enriched. Indels were the most likely causal variant in seven loci, including one locus associated with monocyte count where an indel with causality and mechanism previously demonstrated (rs200748895:TGCTG/T) had a 0.999 posterior probability. Overall, our results show a very modest and nonsignificant enrichment for common indels in associated loci.

Population-scale genomics-Enabling precision public health

The current excitement for affordable genomics technologies and national precision medicine initiatives marks a turning point in worldwide healthcare practices. The last decade of global population sequencing efforts has defined the enormous extent of genetic variation in the human population resulting in insights into differential disease burden and response to therapy within and between populations. Population-scale pharmacogenomics helps to provide insights into the choice of optimal therapies and an opportunity to estimate, predict and minimize adverse events. Such an approach can potentially empower countries to formulate national selection and dosing policies for therapeutic agents thereby promoting public health with precision. We review the breadth and depth of worldwide population-scale sequencing efforts and its implications for the implementation of clinical pharmacogenetics toward making precision medicine a reality.

Genetics of blood lipids among ~300,000 multi-ethnic participants of the Million Veteran Program

The Million Veteran Program (MVP) was established in 2011 as a national research initiative to determine how genetic variation influences the health of US military veterans. Here we genotyped 312,571 MVP participants using a custom biobank array and linked the genetic data to laboratory and clinical phenotypes extracted from electronic health records covering a median of 10.0 years of follow-up. Among 297,626 veterans with at least one blood lipid measurement, including 57,332 black and 24,743 Hispanic participants, we tested up to around 32 million variants for association with lipid levels and identified 118 novel genome-wide significant loci after meta-analysis with data from the Global Lipids Genetics Consortium (total n > 600,000). Through a focus on mutations predicted to result in a loss of gene function and a phenome-wide association study, we propose novel indications for pharmaceutical inhibitors targeting PCSK9 (abdominal aortic aneurysm), ANGPTL4 (type 2 diabetes) and PDE3B (triglycerides and coronary disease).

A Comprehensive Map of Genetic Variation in the World's Largest Ethnic Group-Han Chinese

As are most non-European populations, the Han Chinese are relatively understudied in population and medical genetics studies. From low-coverage whole-genome sequencing of 11,670 Han Chinese women we present a catalog of 25,057,223 variants, including 548,401 novel variants that are seen at least 10 times in our data set. Individuals from this data set came from 24 out of 33 administrative divisions across China (including 19 provinces, 4 municipalities, and 1 autonomous region), thus allowing us to study population structure, genetic ancestry, and local adaptation in Han Chinese. We identified previously unrecognized population structure along the East-West axis of China, demonstrated a general pattern of isolation-by-distance among Han Chinese, and reported unique regional signals of admixture, such as European influences among the Northwestern provinces of China. Furthermore, we identified a number of highly differentiated, putatively adaptive, loci (e.g., MTHFR, ADH7, and FADS, among others) that may be driven by immune response, climate, and diet in the Han Chinese. Finally, we have made available allele frequency estimates stratified by administrative divisions across China in the Geography of Genetic Variant browser for the broader community. By leveraging the largest currently available genetic data set for Han Chinese, we have gained insights into the history and population structure of the world's largest ethnic group.

Is the third molar maturity index (I3M) useful for a genetic isolate population? Study of a Sardinian sample of children and young adults

This work aims to assess the validity of the cut-off value (0.08) of the third molar maturity index (I_3M) for discriminating minors from adults in Sardinian population. A sample of 336 digital panoramic radiographs of healthy Sardinian children and young minors (165 females and 171 males), aged between 15 and 23 years (mean age, 19.35 years in females and 18.80 years in males), was retrospectively evaluated. The left lower third molars were analysed by applying a specific cut-off value of 0.08 determined by Cameriere et al. in 2008. The reliability and reproducibility of the test was also studied: the intra-class correlation coefficients (ICC) were 0.91 (95% CI, 0.89-0.93) and 0.88 (95% CI, 0.86-0.90), for the intra- and inter-observer reliability, respectively. The I_3M gradually decreased as the real age gradually increased in both sexes. According to the pooled results of the diagnostic test, the accuracy (ACC) was 0.86 (95% CI, 0.82-0.89); the proportion of correctly classified subjects (Se = sensitivity) was 0.82 (95% CI: 0.76-0.86); and specificity (Sp = specificity) was 0.95 (95% CI, 0.89-0.97). The positive predictive values (PPV) and the negative predictive values (NPVs) were 0.97 (95% CI, 0.94-0.99) and 0.70 (95% CI, 0.62-0.77). The LR+ and the LR- were 17.12 (95% CI, 7.27 to 40.36) and 0.19 (95% CI, 0.14 to 0.25). In spite of this, significant differences in the early mineralisation of the third molar were found between sexes as well as in the results of the diagnostic test, showing a better sensitivity in males than in females. The results showed that, although the third molar teeth are highly variable in development, and with differences between females and males as compare to other teeth, the I_3M is a reliable method to distinguish between minors and adults even in such a genetic isolate population.

Genomic history of the Sardinian population

The population of the Mediterranean island of Sardinia has made important contributions to genome-wide association studies of complex disease traits and, based on ancient DNA (aDNA) studies of mainland Europe, Sardinia is hypothesized to be a unique refuge for early Neolithic ancestry. To provide new insights on the genetic history of this flagship population, we analyzed 3,514 whole-genome sequenced individuals from Sardinia. We find Sardinian samples show elevated levels of shared ancestry with Basque individuals, especially samples from the more historically isolated regions of Sardinia. Our analysis also uniquely illuminates how levels of genetic similarity with mainland aDNA samples varies subtly across the island. Together, our results indicate within-island sub-structure and sex-biased processes have substantially impacted the genetic history of Sardinia. These results give new insight to the demography of ancestral Sardinians and help further the understanding of sharing of disease risk alleles between Sardinia and mainland populations.

An efficient analytic approach in genome-wide identification of methylation quantitative trait loci response to fenofibrate treatment

Background

The study of DNA methylation quantitative trait loci (meQTLs) helps dissect regulatory mechanisms underlying genetic associations of human diseases. In this study, we conducted the first genome-wide examination of genetic drivers of methylation variation in response to a triglyceride-lowering treatment with fenofibrate (response-meQTL) by using an efficient analytic approach.

Methods

Subjects (n = 429) from the GAW20 real data set with genotype and both pre- (visit 2) and post- (visit 4) fenofibrate treatment methylation measurements were included. Following the quality control steps of removing certain cytosine-phosphate-guanine (CpG) probes, the post−/premethylation changes (post/pre) were log transformed and the association was performed on 208,449 CpG sites. An additive linear mixed-effects model was used to test the association between each CpG probe and single nucleotide polymorphisms (SNPs) around ±1 Mb region, with age, sex, smoke, batch effect, and principal components included as covariates. Bonferroni correction was applied to define the significance threshold (p < 5.6 × 10^− 10, given a total of 89,217,303 tests). Finally, we integrated our response-meQTL (re-meQTL) findings with the published genome-wide association study (GWAS) catalog of human diseases/traits.

Results

We identified 1087 SNPs as cis re-meQTLs associated with 610 CpG probes/sites located in 351 unique gene loci. Among these 1087 cis re-meQTL SNPs, 229 were unique and 6 were co-localized at 8 unique disease/trait loci reported in the GWAS catalog (enrichment p = 1.51 × 10^− 23). Specifically, a lipid SNP, rs10903129, located in intron regions of gene TMEM57, was a re-meQTL (p = 3.12 × 10^− 36) associated with the CpG probe cg09222892, which is in the upstream region of the gene RHCE, indicating a new target gene for rs10903129. In addition, we found that SNP rs12710728 has a suggestive association with cg17097782 (p = 1.77 × 10^− 4), and that this SNP is in high linkage disequilibrium (LD) (R² > 0.8) with rs7443270, which was previously reported to be associated with fenofibrate response (p = 5.00 × 10^− 6).

Conclusions

By using a novel analytic approach, we efficiently identified thousands of cis re-meQTLs that provide a unique resource for further characterizing functional roles and gene targets of the SNPs that are most responsive to fenofibrate treatment. Our efficient analytic approach can be extended to large response quantitative trait locus studies with large sample sizes and multiple time points data.

Exploring rare and low-frequency variants in the Saguenay-Lac-Saint-Jean population identified genes associated with asthma and allergy traits

The Saguenay–Lac-Saint-Jean (SLSJ) region is located in northeastern Quebec and is known for its unique demographic history and founder effect. As founder populations are enriched with population-specific variants, we characterized the variants distribution in SLSJ and compared it with four European populations (Finnish, Sweden, United Kingdom and France), of which the Finnish population is another founder population. Targeted sequencing of the coding and non-coding immune regulatory regions of the SLSJ asthma familial cohort and the four European populations were performed. Rare and low-frequency coding and non-coding regulatory variants identified in the SLSJ population were then investigated for variant- and gene-level associations with asthma and allergy-related traits (eosinophil percentage, immunoglobulin (Ig) E levels and lung function). Our data showed that (1) rare or deleterious variants were not enriched in the two founder populations as compared with the three non-founder European populations; (2) a larger proportion of founder population-specific variants occurred with higher frequencies; and (3) low-frequency variants appeared to be more deleterious. Furthermore, a rare variant, rs1386931, located in the 3ʹ-UTR of CXCR6 and intron of FYCO1 was found to be associated with eosinophil percentage. Gene-based analyses identified NRP2, MRPL44 and SERPINE2 to be associated with various asthma and allergy-related traits. Our study demonstrated the usefulness of using a founder population to identify new genes associated with asthma and allergy-related traits; thus better understand the genes and pathways implicated in pathophysiology.

Discovery, fine-mapping, and conditional analyses of genetic variants associated with C-reactive protein in multiethnic populations using the Metabochip in the Population Architecture using Genomics and Epidemiology (PAGE) study

C-reactive protein (CRP) is a circulating biomarker indicative of systemic inflammation. We aimed to evaluate genetic associations with CRP levels among non-European-ancestry populations through discovery, fine-mapping and conditional analyses. A total of 30 503 non-European-ancestry participants from 6 studies participating in the Population Architecture using Genomics and Epidemiology study had serum high-sensitivity CRP measurements and ∼200 000 single nucleotide polymorphisms (SNPs) genotyped on the Metabochip. We evaluated the association between each SNP and log-transformed CRP levels using multivariate linear regression, with additive genetic models adjusted for age, sex, the first four principal components of genetic ancestry, and study-specific factors. Differential linkage disequilibrium patterns between race/ethnicity groups were used to fine-map regions associated with CRP levels. Conditional analyses evaluated for multiple independent signals within genetic regions. One hundred and sixty-three unique variants in 12 loci in overall or race/ethnicity-stratified Metabochip-wide scans reached a Bonferroni-corrected P-value <2.5E-7. Three loci have no (HACL1, OLFML2B) or only limited (PLA2G6) previous associations with CRP levels. Six loci had different top hits in race/ethnicity-specific versus overall analyses. Fine-mapping refined the signal in six loci, particularly in HNF1A. Conditional analyses provided evidence for secondary signals in LEPR, IL1RN and HNF1A, and for multiple independent signals in CRP and APOE. We identified novel variants and loci associated with CRP levels, generalized known CRP associations to a multiethnic study population, refined association signals at several loci and found evidence for multiple independent signals at several well-known loci. This study demonstrates the benefit of conducting inclusive genetic association studies in large multiethnic populations.

Genetic-Driven Druggable Target Identification and Validation

Choosing the right biological target is the critical primary decision for the development of new drugs. Systematic genetic association testing of both human diseases and quantitative traits, along with resultant findings of coincident associations between them, is becoming a powerful approach to infer drug targetable candidates and generate in vitro tests to identify compounds that can modulate them therapeutically. Here, we discuss opportunities and challenges, and infer criteria for the optimal use of genetic findings in the drug discovery pipeline.

Genetic variants in mRNA untranslated regions

Genome Wide Association Studies (GWAS) have mapped thousands of genetic variants associated with complex disease risk and regulating quantitative traits, thus exploiting an unprecedented high-resolution genetic characterization of the human genome. A small fraction (3.7%) of the identified associations is located in untranslated regions (UTRs), and the molecular mechanism has been elucidated for few of them. Genetic variations at UTRs may modify regulatory elements affecting the interaction of the UTRs with proteins and microRNAs. The overall functional consequences include modulation of messenger RNA (mRNA) transcription, secondary structure, stability, localization, translation, and access to regulators like microRNAs (miRNAs) and RNA-binding proteins (RBPs). Alterations of these regulatory mechanisms are known to modify molecular pathways and cellular processes, potentially leading to disease processes. Here, we analyze some examples of genetic risk variants mapping in the UTR regulatory elements. We describe a recently identified genetic variant localized in the 3'UTR of the TNFSF13B gene, associated with autoimmunity risk and responsible of an increased stability and translation of TNFSF13B mRNA. We discuss how the correct use and interpretation of public GWAS repositories could lead to a better understanding of etiopathogenetic mechanisms and the generation of robust biological hypothesis as starting point for further functional studies. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Evolution and Genomics > Computational Analyses of RNA RNA in Disease and Development > RNA in Disease.

A large electronic-health-record-based genome-wide study of serum lipids

A genome-wide association study of 94,674 multi-ethnic Kaiser Permanente members utilizing 478,866 longitudinal untreated serum lipid electronic-health-record-derived measurements (EHRs) empowered multiple novel findings: 121 new SNP associations (46 primary, 15 conditional, 60 in meta-analysis with Global Lipids Genetic Consortium); increase of 33-42% in variance explained with multiple measurements; sex differences in genetic impact (greater in females for LDL, HDL, TC, the opposite for TG); differences in variance explained amongst non-Hispanic whites, Latinos, African Americans, and East Asians; genetic dominance and epistasis, with strong evidence for both at ABOxFUT2 for LDL; and eQTL tissue-enrichment implicating the liver, adipose, and pancreas. Utilizing EHR pharmacy data, both LDL and TG genetic risk scores (477 SNPs) were strongly predictive of age-at-initiation of lipid-lowering treatment. These findings highlight the value of longitudinal EHRs for identifying novel genetic features of cholesterol and lipoprotein metabolism with implications for lipid treatment and risk of coronary heart disease.