Large-scale whole-genome sequencing of the Icelandic population. Nat Genet. 2015;47:435-444. [PMID: 25807286 DOI: 10.1038/ng.3247]

Myosin light-chain 4 gene-transfer attenuates atrial fibrosis while correcting autophagic flux dysregulation

OBJECTIVES

To determine the role of MYL4 regulation of lysosomal function and its disturbance in fibrotic atrial cardiomyopathy.

BACKGROUND

We have previously demonstrated that the atrial-specific essential light chain protein MYL4 is required for atrial contractile, electrical, and structural integrity. MYL4 mutation/dysfunction leads to atrial fibrosis, standstill, and dysrhythmia. However, the underlying pathogenic mechanisms remain unclear.

METHODS AND RESULTS

Rats subjected to knock-in of a pathogenic MYL4 mutant (p.E11K) developed fibrotic atrial cardiomyopathy. Proteome analysis and single-cell RNA sequencing indicate enrichment of autophagy pathways in mutant-MYL4 atrial dysfunction. Immunofluorescence and electron microscopy revealed undegraded autophagic vesicles accumulated in MYL4p.E11K rat atrium. Next, we identified that dysfunctional MYL4 protein impairs autophagy flux in vitro and in vivo. Cardiac lysosome positioning and mobility were regulated by MYL4 in cardiomyocytes, which affected lysosomal acidification and maturation of lysosomal cathepsins. We then examined the effects of MYL4 overexpression via adenoviral gene-transfer on atrial cardiomyopathy induced by MYL4 mutation: MYL4 protein overexpression attenuated atrial structural remodeling and autophagy dysfunction.

CONCLUSIONS

MYL4 regulates autophagic flux in atrial cardiomyocytes via lysosomal mobility. MYL4 overexpression attenuates MYL4 p.E11K induced fibrotic atrial cardiomyopathy, while correcting autophagy and lysosomal function. These results provide a molecular basis for MYL4-mutant induced fibrotic atrial cardiomyopathy and identify a potential biological-therapy approach for the treatment of atrial fibrosis.

Molecular and Clinical Links between Drug-Induced Cholestasis and Familial Intrahepatic Cholestasis

Idiosyncratic Drug-Induced Liver Injury (iDILI) represents an actual health challenge, accounting for more than 40% of hepatitis cases in adults over 50 years and more than 50% of acute fulminant hepatic failure cases. In addition, approximately 30% of iDILI are cholestatic (drug-induced cholestasis (DIC)). The liver's metabolism and clearance of lipophilic drugs depend on their emission into the bile. Therefore, many medications cause cholestasis through their interaction with hepatic transporters. The main canalicular efflux transport proteins include: 1. the bile salt export pump (BSEP) protein (ABCB11); 2. the multidrug resistance protein-2 (MRP2, ABCC2) regulating the bile salts' independent flow by excretion of glutathione; 3. the multidrug resistance-1 protein (MDR1, ABCB1) that transports organic cations; 4. the multidrug resistance-3 protein (MDR3, ABCB4). Two of the most known proteins involved in bile acids' (BAs) metabolism and transport are BSEP and MDR3. BSEP inhibition by drugs leads to reduced BAs' secretion and their retention within hepatocytes, exiting in cholestasis, while mutations in the ABCB4 gene expose the biliary epithelium to the injurious detergent actions of BAs, thus increasing susceptibility to DIC. Herein, we review the leading molecular pathways behind the DIC, the links with the other clinical forms of familial intrahepatic cholestasis, and, finally, the main cholestasis-inducing drugs.

Developmental language disorder - a comprehensive study of more than 46,000 individuals

Developmental language disorder (DLD) is characterized by enduring low language abilities with a significant functional impact, in the absence of biomedical conditions in which language impairment is part of a complex of impairments. There is a lack of awareness of DLD even among healthcare professionals. Here we estimated the prevalence of DLD and its links to reading and learning difficulties and physical and mental health in the Danish Blood Donor Study (N = 46,547), where DLD-related information is based on questionnaires (self-report). We compared the questionnaire-derived DLD status with the relevant language-related diagnoses from hospital registers. We also investigated the genetic architecture of DLD in a subset of the cohort (N = 18,380). DLD was significantly associated with reading and learning difficulties and poorer mental and physical health. DLD prevalence was 3.36%-3.70% based on questionnaires, compared with 0.04% in hospital registers. Our genetic analyses identified one genome-wide significant locus, but not a significant heritability estimate. Our study shows that DLD has health-related implications that may last into adulthood, and that DLD may be undiagnosed in general healthcare. Furthermore, DLD is likely more genetically heterogeneous than narrower developmental language phenotypes. Our results emphasize the need to raise awareness of DLD and consider criteria for molecular studies of DLD to reduce case heterogeneity.

Connexin hemichannels as candidate targets for cardioprotective and anti-arrhythmic treatments

Connexins are crucial cardiac proteins that form hemichannels and gap junctions. Gap junctions are responsible for the propagation of electrical and chemical signals between myocardial cells and cells of the specialized conduction system in order to synchronize the cardiac cycle and steer cardiac pump function. Gap junctions are normally open, while hemichannels are closed, but pathological circumstances may close gap junctions and open hemichannels, thereby perturbing cardiac function and homeostasis. Current evidence demonstrates an emerging role of hemichannels in myocardial ischemia and arrhythmia, and tools are now available to selectively inhibit hemichannels without inhibiting gap junctions as well as to stimulate hemichannel incorporation into gap junctions. We review available experimental evidence for hemichannel contributions to cellular pro-arrhythmic events in ventricular and atrial cardiomyocytes, and link these to insights at the level of molecular control of connexin-43-based hemichannel opening. We conclude that a double-edged approach of both preventing hemichannel opening and preserving gap junctional function will be key for further research and development of new connexin-based experimental approaches for treating heart disease.

A crowdsourcing database for the copy-number variation of the Spanish population

BACKGROUND

Despite being a very common type of genetic variation, the distribution of copy-number variations (CNVs) in the population is still poorly understood. The knowledge of the genetic variability, especially at the level of the local population, is a critical factor for distinguishing pathogenic from non-pathogenic variation in the discovery of new disease variants.

RESULTS

Here, we present the SPAnish Copy Number Alterations Collaborative Server (SPACNACS), which currently contains copy number variation profiles obtained from more than 400 genomes and exomes of unrelated Spanish individuals. By means of a collaborative crowdsourcing effort whole genome and whole exome sequencing data, produced by local genomic projects and for other purposes, is continuously collected. Once checked both, the Spanish ancestry and the lack of kinship with other individuals in the SPACNACS, the CNVs are inferred for these sequences and they are used to populate the database. A web interface allows querying the database with different filters that include ICD10 upper categories. This allows discarding samples from the disease under study and obtaining pseudo-control CNV profiles from the local population. We also show here additional studies on the local impact of CNVs in some phenotypes and on pharmacogenomic variants. SPACNACS can be accessed at: http://csvs.clinbioinfosspa.es/spacnacs/ .

CONCLUSION

SPACNACS facilitates disease gene discovery by providing detailed information of the local variability of the population and exemplifies how to reuse genomic data produced for other purposes to build a local reference database.

Heart Disease and Stroke Statistics-2023 Update: A Report From the American Heart Association

BACKGROUND

The American Heart Association, in conjunction with the National Institutes of Health, annually reports the most up-to-date statistics related to heart disease, stroke, and cardiovascular risk factors, including core health behaviors (smoking, physical activity, diet, and weight) and health factors (cholesterol, blood pressure, and glucose control) that contribute to cardiovascular health. The Statistical Update presents the latest data on a range of major clinical heart and circulatory disease conditions (including stroke, congenital heart disease, rhythm disorders, subclinical atherosclerosis, coronary heart disease, heart failure, valvular disease, venous disease, and peripheral artery disease) and the associated outcomes (including quality of care, procedures, and economic costs).

METHODS

The American Heart Association, through its Epidemiology and Prevention Statistics Committee, continuously monitors and evaluates sources of data on heart disease and stroke in the United States to provide the most current information available in the annual Statistical Update with review of published literature through the year before writing. The 2023 Statistical Update is the product of a full year's worth of effort in 2022 by dedicated volunteer clinicians and scientists, committed government professionals, and American Heart Association staff members. The American Heart Association strives to further understand and help heal health problems inflicted by structural racism, a public health crisis that can significantly damage physical and mental health and perpetuate disparities in access to health care, education, income, housing, and several other factors vital to healthy lives. This year's edition includes additional COVID-19 (coronavirus disease 2019) publications, as well as data on the monitoring and benefits of cardiovascular health in the population, with an enhanced focus on health equity across several key domains.

RESULTS

Each of the chapters in the Statistical Update focuses on a different topic related to heart disease and stroke statistics.

CONCLUSIONS

The Statistical Update represents a critical resource for the lay public, policymakers, media professionals, clinicians, health care administrators, researchers, health advocates, and others seeking the best available data on these factors and conditions.

The cancer-risk variant frequency among Polish population reported by the first national whole-genome sequencing study

Introduction

Population-based cancer screening has raised many controversies in recent years, not only regarding the costs but also regarding the ethical nature and issues related to variant interpretation. Nowadays, genetic cancer screening standards are different in every country and usually encompass only individuals with a personal or family history of relevant cancer.

Methods

Here we performed a broad genetic screening for cancer-related rare germline variants on population data from the Thousand Polish Genomes database based on 1076 Polish unrelated individuals that underwent whole genome sequencing (WGS).

Results

We identified 19 551 rare variants in 806 genes related to oncological diseases, among them 89% have been located in non-coding regions. The combined BRCA1/BRCA2 pathogenic/likely pathogenic according to ClinVar allele frequency in the unselected population of 1076 Poles was 0.42%, corresponding to nine carriers.

Discussion

Altogether, on the population level, we found especially problematic the assessment of the pathogenicity of variants and the relation of ACMG guidelines to the population frequency. Some of the variants may be overinterpreted as disease-causing due to their rarity or lack of annotation in the databases. On the other hand, some relevant variants may have been overseen given that there is little pooled population whole genome data on oncology. Before population WGS screening will become a standard, further studies are needed to assess the frequency of the variants suspected to be pathogenic on the population level and with reporting of likely benign variants.

WGS Data Collections: How Do Genomic Databases Transform Medicine?

As a scientific community we assumed that exome sequencing will elucidate the basis of most heritable diseases. However, it turned out it was not the case; therefore, attention has been increasingly focused on the non-coding sequences that encompass 98% of the genome and may play an important regulatory function. The first WGS-based datasets have already been released including underrepresented populations. Although many databases contain pooled data from several cohorts, recently the importance of local databases has been highlighted. Genomic databases are not only collecting data but may also contribute to better diagnostics and therapies. They may find applications in population studies, rare diseases, oncology, pharmacogenetics, and infectious and inflammatory diseases. Further data may be analysed with Al technologies and in the context of other omics data. To exemplify their utility, we put a highlight on the Polish genome database and its practical application.

Genome-wide association study of varicose veins identifies a protective missense variant in GJD3 enriched in the Finnish population

Varicose veins is the most common manifestation of chronic venous disease that displays female-biased incidence. To identify protein-inactivating variants that could guide identification of drug target genes for varicose veins and genetic evidence for the disease prevalence difference between the sexes, we conducted a genome-wide association study of varicose veins in Finns using the FinnGen dataset with 17,027 cases and 190,028 controls. We identified 50 associated genetic loci (P < 5.0 × 10-8) of which 29 were novel including one near ERG with female-specificity (rs2836405-G, OR[95% CI] = 1.09[1.05-1.13], P = 3.1 × 10-8). These also include two X-chromosomal (ARHGAP6 and SRPX) and two autosomal novel loci (TGFB2 and GJD3) with protein-coding lead variants enriched above 56-fold in Finns over non-Finnish non-Estonian Europeans. A low-frequency missense variant in GJD3 (p.Pro59Thr) is exclusively associated with a lower risk for varicose veins (OR = 0.62 [0.55-0.70], P = 1.0 × 10-14) in a phenome-wide scan of the FinnGen data. The absence of observed pleiotropy and its membership of the connexin gene family underlines GJD3 as a potential connexin-modulating therapeutic strategy for varicose veins. Our results provide insights into varicose veins etiopathology and highlight the power of isolated populations, including Finns, to discover genetic variants that inform therapeutic development.

Population Genome Programs across the Middle East and North Africa: Successes, Challenges, and Future Directions

In this review, we discuss the current state of population genome programs (PGPs) conducted in the Middle East and North African (MENA) region. This region has high prevalence of genetic diseases and significant health challenges as well as being a significantly underrepresented population in public genetic databases. The majority of ongoing PGPs represent regions in Europe, North and South America, South Asia, Australia, and Africa, with little to no descriptive information highlighted only on the MENA Region when it comes to genome programs databases, outcomes, or the challenges that MENA region countries may face establishing their own national programs. This review has identified 6 PGPs currently underway in the MENA region, namely in the Kingdom of Saudi Arabia, Qatar, Egypt, the United Arab Emirates, Bahrain, and Iran. Due to the rapidly growing involvement of the MENA region in national-scale genomic data collection, an increase in representation in public genetic databases is to be expected to occur in the near future. Whilst significant progress is being made in some MENA countries, future initiatives as well as ongoing programs will be facing several challenges related to collaboration, finance, infrastructure and institutional data access, data analysis, sustainability, health records, and biobanks. The review also reiterates the need for ensuring ethical and regulated genomic initiatives which can drive developments in personalized medicine treatments to improve patient prognosis and quality of life.

No evidence of a genetic causal relationship between ankylosing spondylitis and iron homeostasis: A two-sample Mendelian randomization study

Background

Ankylosing spondylitis (AS) is an immune-mediated chronic inflammatory disease that leads to bone hyperplasia and spinal ankylosis. Iron homeostasis plays a very important role in the inflammatory response and is closely related to the pathogenesis of AS. This study aimed to use large-scale genome-wide association study (GWAS) summary data to study the genetic causal relationship between AS and iron homeostasis using Mendelian randomization (MR).

Methods

Genome-wide association study summary data of AS and iron homeostasis-related indicators were obtained from the FinnGen consortium and the DeCODE genetics database, respectively. We used four iron homeostasis-related indicators: ferritin, serum iron, total iron binding capacity (TIBC), and transferrin saturation (TSAT) for two-sample MR analyses to test for genetic causal association with AS using the "TwoSampleMR" package of the R software (version 4.1.2). The random-effects inverse variance weighted (IVW) method was the main analysis method used for MR. We examined the MR analysis results for heterogeneity, horizontal pleiotropy, and possible outliers. In addition, we confirmed the robustness of the MR analysis by testing whether the results were affected by a single SNP and whether they followed a normal distribution.

Results

The random-effects IVW results showed that ferritin [p = 0.225, OR 95% confidence interval (CI) = 0.836 (0.627-1.116)], serum iron [p = 0.714, OR 95% CI = 0.948 (0.714-1.260)], TIBC [p = 0.380, OR 95% CI = 0.917 (0.755-1.113)], and TSAT [p = 0.674, OR 95% CI = 0.942 (0.713-1.244)] have no genetic causal relationship with AS. We detected no heterogeneity，horizontal pleiotropy and possible outliers in our MR analysis (p > 0.05). In addition, our MR analysis results were not affected by a single SNP, and were normally distributed.

Conclusion

Our study did not detect a genetic causal relationship between AS and iron homeostasis. Nonetheless, this does not rule out a relationship between the two at other mechanistic levels.

Challenges Related to the Use of Next-Generation Sequencing for the Optimization of Drug Therapy

Over the last decade, next-generation sequencing (NGS) methods have become increasingly used in various areas of human genomics. In routine clinical care, their use is already implemented in oncology to profile the mutational landscape of a tumor, as well as in rare disease diagnostics. However, its utilization in pharmacogenomics is largely lacking behind. Recent population-scale genome data has revealed that human pharmacogenes carry a plethora of rare genetic variations that are not interrogated by conventional array-based profiling methods and it is estimated that these variants could explain around 30% of the genetically encoded functional pharmacogenetic variability.To interpret the impact of such variants on drug response a multitude of computational tools have been developed, but, while there have been major advancements, it remains to be shown whether their accuracy is sufficient to improve personalized pharmacogenetic recommendations in robust trials. In addition, conventional short-read sequencing methods face difficulties in the interrogation of complex pharmacogenes and high NGS test costs require stringent evaluations of cost-effectiveness to decide about reimbursement by national healthcare programs. Here, we illustrate current challenges and discuss future directions toward the clinical implementation of NGS to inform genotype-guided decision-making.

Correlation between Serum Myosin Light Chain 4 Levels and Recurrence after Radiofrequency Ablation in Patients with Atrial Fibrillation

Atrial fibrillation (AF) is the most common arrhythmia that is harmful to human health. This study aims to explore the relationship between myosin light chain 4 (MYL4) and AF recurrence after radiofrequency ablation (RFA). Patients with AF (n = 85) were enrolled, and healthy subjects (n = 90) with normal sinus rhythm and no previous history of AF were selected as controls. The serum levels of MYL4, transforming growth factor (TGF) -β1, and procollagen type-I C-terminal propeptide (PICP) were determined. The correlation between MYL4 and atrial fibrosis remodeling indicators (TGF-β1/PICP) and left atrial diameter (LAD) was analyzed. The influence of MYL4 on AF recurrence after RFA was evaluated, and the independent correlation between them was assessed. Patients with AF and the controls showed no significant differences in age, gender, body mass index, systolic blood pressure, diastolic blood pressure, left ventricular ejection fraction, triglycerides, total cholesterol, high-density lipoprotein, low-density lipoprotein, white blood cell count, neutrophil/lymphocyte ratio, brain natriuretic peptide, and history of smoking, drinking, hypertension, and diabetes (P > 0.05), but with increased LAD in patients with AF (P < 0.01). Serum MYL4 level was reduced in patients with AF (0.6 ± 0.2) compared with that of controls (0.1 ± 0.6) (P < 0.01), and it was negatively correlated with TGF-β1, PICP, and LAD (r = -0.2389, P < 0.05; r = -0.5174, P < 0.01; r = -0.3191; P < 0.01). Low levels of MYL4 increased the risk of AF recurrence after RFA (χ² = 16.64; P < 0.0001). A low MYL4 level in patients with AF showed a poorer prognosis. Serum MYL4 level and AF type were independent risk factors affecting AF recurrence after RFA.

Whole genome sequencing identifies structural variants contributing to hematologic traits in the NHLBI TOPMed program

Genome-wide association studies have identified thousands of single nucleotide variants and small indels that contribute to variation in hematologic traits. While structural variants are known to cause rare blood or hematopoietic disorders, the genome-wide contribution of structural variants to quantitative blood cell trait variation is unknown. Here we utilized whole genome sequencing data in ancestrally diverse participants of the NHLBI Trans Omics for Precision Medicine program (N = 50,675) to detect structural variants associated with hematologic traits. Using single variant tests, we assessed the association of common and rare structural variants with red cell-, white cell-, and platelet-related quantitative traits and observed 21 independent signals (12 common and 9 rare) reaching genome-wide significance. The majority of these associations (N = 18) replicated in independent datasets. In genome-editing experiments, we provide evidence that a deletion associated with lower monocyte counts leads to disruption of an S1PR3 monocyte enhancer and decreased S1PR3 expression.

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.

Polygenic risk score for ACE-inhibitor-associated cough based on the discovery of new genetic loci

AIMS

To search for sequence variants associated with ACEi discontinuation and to test their association with ACEi-associated adverse drug reactions (ADRs).

METHODS AND RESULTS

A genome-wide association study (GWAS) on ACEi discontinuation was conducted, including 33 959 ACEi-discontinuers and 44 041 controls. Cases were defined as persons who switched from an ACEi treatment to an angiotensin receptor blocker. Controls were defined as persons who continued ACEi treatment for at least 1 year. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were computed for ACEi discontinuation risk by mixed model regression analysis. Summary statistics from the individual cohorts were meta-analyzed with a fixed-effects model. To test for association with specific ACEi-associated ADRs, any genome-wide significant (P < 5 × 10-8) ACEi discontinuation variants was tested for association with ACEi-associated cough and angioedema. A polygenetic risk score (PRS) based on ACEi discontinuation GWAS data was constructed and tested for association with ACEi-associated cough and angioedema in two population-based samples. In total, seven genetic genome-wide loci were identified, of which six were previously unreported. The strongest association with ACEi discontinuation was at 20q13.3 (NTSR1; OR: 1.21; 95% CI: 1.17-1.24; P = 2.1 × 10-34). Five of seven lead variants were associated with ACEi-associated cough, whereas none were associated with ACEi-associated angioedema. The ACEi discontinuation PRS was associated with ACEi-associated cough in a dose-response manner but not with ACEi-associated angioedema. ACEi discontinuation was genetically correlated with important causes for cough, including gastro-esophageal reflux disease, allergic rhinitis, hay fever, and asthma, which indicates partly shared genetic underpinning between these traits.

CONCLUSION

This study showed the advantage of using prescription patterns to discover genetic links with ADRs. In total, seven genetic loci that associated with ACEi discontinuation were identified. There was evidence of a strong association between our ADR phenotype and ACEi-associated cough. Taken together, these findings increase insight into the pathophysiological processes that underlie ACEi-associated ADRs.

A common variant in the hepatobiliary phospholipid transporter ABCB4 modulates liver injury in PBC but not in PSC: prospective analysis in 867 patients

BACKGROUND

The ATP-binding cassette subfamily B member 4 (ABCB4) gene encodes the hepatic phospholipid transporter. Variants in the ABCB4 gene are associated with various cholestatic phenotypes, some of which progress to liver fibrosis and cirrhosis. The aim of our study was to investigate the role of the cholestasis-associated variant ABCB4 c.711A > T (p.I237I, rs2109505) in patients with primary biliary cholangitis (PBC) and primary sclerosing cholangitis (PSC).

RESULTS

Two cohorts of Polish patients took part in this study. The Szczecin cohort comprised 196 patients with PBC (174 females, 38% with cirrhosis) and 135 patients with PSC (39 females, 39% with cirrhosis). The Warsaw cohort consisted of 260 patients with PBC (241 females, 44% with cirrhosis) and 276 patients with PSC (97 females, 33% with cirrhosis). Two control cohorts-150 healthy blood donors and 318 patients without liver disease, were recruited in Szczecin and in Warsaw, respectively. The ABCB4 c.711A > T polymorphism was genotyped using TaqMan assay. In both PBC cohorts, carriers of the risk variant presented more frequently with cirrhosis (Szczecin: OR = 1.841, P = 0.025; Warsaw: OR = 1.528, P = 0.039). The risk allele was associated with increased serum AST, GGT and ALP (all P < 0.05) at inclusion. During the follow-up, patients in both cohorts significantly improved their laboratory results, independently of their ABCB4 c.711A > T genotype (P > 0.05). During 8 ± 4 years follow-up, a total of 22 patients in the Szczecin PBC group developed cirrhosis, and this risk was higher among carriers of the risk variant (OR = 5.65, P = 0.04). In contrast to PBC, we did not detect any association of ABCB4 c.711A > T with a liver phenotype in PSC cohorts.

CONCLUSIONS

The frequent pro-cholestatic variant ABCB4 c.711A > T modulates liver injury in PBC, but not in PSC. In particular, carriers of the major allele are at increased risk of progressive liver scarring.

Multiomics study of nonalcoholic fatty liver disease

Nonalcoholic fatty liver (NAFL) and its sequelae are growing health problems. We performed a genome-wide association study of NAFL, cirrhosis and hepatocellular carcinoma, and integrated the findings with expression and proteomic data. For NAFL, we utilized 9,491 clinical cases and proton density fat fraction extracted from 36,116 liver magnetic resonance images. We identified 18 sequence variants associated with NAFL and 4 with cirrhosis, and found rare, protective, predicted loss-of-function variants in MTARC1 and GPAM, underscoring them as potential drug targets. We leveraged messenger RNA expression, splicing and predicted coding effects to identify 16 putative causal genes, of which many are implicated in lipid metabolism. We analyzed levels of 4,907 plasma proteins in 35,559 Icelanders and 1,459 proteins in 47,151 UK Biobank participants, identifying multiple proteins involved in disease pathogenesis. We show that proteomics can discriminate between NAFL and cirrhosis. The present study provides insights into the development of noninvasive evaluation of NAFL and new therapeutic options.

Vasor: Accurate prediction of variant effects for amino acid substitutions in multidrug resistance protein 3

The phosphatidylcholine floppase multidrug resistance protein 3 (MDR3) is an essential hepatobiliary transport protein. MDR3 dysfunction is associated with various liver diseases, ranging from severe progressive familial intrahepatic cholestasis to transient forms of intrahepatic cholestasis of pregnancy and familial gallstone disease. Single amino acid substitutions are often found as causative of dysfunction, but identifying the substitution effect in in vitro studies is time and cost intensive. We developed variant assessor of MDR3 (Vasor), a machine learning-based model to classify novel MDR3 missense variants into the categories benign or pathogenic. Vasor was trained on the largest data set to date that is specific for benign and pathogenic variants of MDR3 and uses general predictors, namely Evolutionary Models of Variant Effects (EVE), EVmutation, PolyPhen-2, I-Mutant2.0, MUpro, MAESTRO, and PON-P2 along with other variant properties, such as half-sphere exposure and posttranslational modification site, as input. Vasor consistently outperformed the integrated general predictors and the external prediction tool MutPred2, leading to the current best prediction performance for MDR3 single-site missense variants (on an external test set: F1-score, 0.90; Matthew's correlation coefficient, 0.80). Furthermore, Vasor predictions cover the entire sequence space of MDR3. Vasor is accessible as a webserver at https://cpclab.uni-duesseldorf.de/mdr3_predictor/ for users to rapidly obtain prediction results and a visualization of the substitution site within the MDR3 structure. The MDR3-specific prediction tool Vasor can provide reliable predictions of single-site amino acid substitutions, giving users a fast way to initially assess whether a variant is benign or pathogenic.

Gallbladder cancer

Gallbladder cancer (GBC) is the most common cancer of the biliary tract, characterized by a very poor prognosis when diagnosed at advanced stages owing to its aggressive behaviour and limited therapeutic options. Early detection at a curable stage remains challenging because patients rarely exhibit symptoms; indeed, most GBCs are discovered incidentally following cholecystectomy for symptomatic gallbladder stones. Long-standing chronic inflammation is an important driver of GBC, regardless of the lithiasic or non-lithiasic origin. Advances in omics technologies have provided a deeper understanding of GBC pathogenesis, uncovering mechanisms associated with inflammation-driven tumour initiation and progression. Surgical resection is the only treatment with curative intent for GBC but very few cases are suitable for resection and most adjuvant therapy has a very low response rate. Several unmet clinical needs require to be addressed to improve GBC management, including discovery and validation of reliable biomarkers for screening, therapy selection and prognosis. Standardization of preneoplastic and neoplastic lesion nomenclature, as well as surgical specimen processing and sampling, now provides reproducible and comparable research data that provide a basis for identifying and implementing early detection strategies and improving drug discovery. Advances in the understanding of next-generation sequencing, multidisciplinary care for GBC, neoadjuvant and adjuvant strategies, and novel systemic therapies including chemotherapy and immunotherapies are gradually changing the treatment paradigm and prognosis of this recalcitrant cancer.

Multi-ancestry genome-wide association study of gestational diabetes mellitus highlights genetic links with type 2 diabetes

Gestational diabetes mellitus (GDM) is associated with increased risk of pregnancy complications and adverse perinatal outcomes. GDM often reoccurs and is associated with increased risk of subsequent diagnosis of type 2 diabetes (T2D). To improve our understanding of the aetiological factors and molecular processes driving the occurrence of GDM, including the extent to which these overlap with T2D pathophysiology, the GENetics of Diabetes In Pregnancy Consortium assembled genome-wide association studies of diverse ancestry in a total of 5485 women with GDM and 347 856 without GDM. Through multi-ancestry meta-analysis, we identified five loci with genome-wide significant association (P < 5 × 10-8) with GDM, mapping to/near MTNR1B (P = 4.3 × 10-54), TCF7L2 (P = 4.0 × 10-16), CDKAL1 (P = 1.6 × 10-14), CDKN2A-CDKN2B (P = 4.1 × 10-9) and HKDC1 (P = 2.9 × 10-8). Multiple lines of evidence pointed to the shared pathophysiology of GDM and T2D: (i) four of the five GDM loci (not HKDC1) have been previously reported at genome-wide significance for T2D; (ii) significant enrichment for associations with GDM at previously reported T2D loci; (iii) strong genetic correlation between GDM and T2D and (iv) enrichment of GDM associations mapping to genomic annotations in diabetes-relevant tissues and transcription factor binding sites. Mendelian randomization analyses demonstrated significant causal association (5% false discovery rate) of higher body mass index on increased GDM risk. Our results provide support for the hypothesis that GDM and T2D are part of the same underlying pathology but that, as exemplified by the HKDC1 locus, there are genetic determinants of GDM that are specific to glucose regulation in pregnancy.

De novo mutations, genetic mosaicism and human disease

Mosaicism—the existence of genetically distinct populations of cells in a particular organism—is an important cause of genetic disease. Mosaicism can appear as de novo DNA mutations, epigenetic alterations of DNA, and chromosomal abnormalities. Neurodevelopmental or neuropsychiatric diseases, including autism—often arise by de novo mutations that usually not present in either of the parents. De novo mutations might occur as early as in the parental germline, during embryonic, fetal development, and/or post-natally, through ageing and life. Mutation timing could lead to mutation burden of less than heterozygosity to approaching homozygosity. Developmental timing of somatic mutation attainment will affect the mutation load and distribution throughout the body. In this review, we discuss the timing of de novo mutations, spanning from mutations in the germ lineage (all ages), to post-zygotic, embryonic, fetal, and post-natal events, through aging to death. These factors can determine the tissue specific distribution and load of de novo mutations, which can affect disease. The disease threshold burden of somatic de novo mutations of a particular gene in any tissue will be important to define.

Genomic prediction with whole-genome sequence data in intensely selected pig lines

Background

Early simulations indicated that whole-genome sequence data (WGS) could improve the accuracy of genomic predictions within and across breeds. However, empirical results have been ambiguous so far. Large datasets that capture most of the genomic diversity in a population must be assembled so that allele substitution effects are estimated with high accuracy. The objectives of this study were to use a large pig dataset from seven intensely selected lines to assess the benefits of using WGS for genomic prediction compared to using commercial marker arrays and to identify scenarios in which WGS provides the largest advantage.

Methods

We sequenced 6931 individuals from seven commercial pig lines with different numerical sizes. Genotypes of 32.8 million variants were imputed for 396,100 individuals (17,224 to 104,661 per line). We used BayesR to perform genomic prediction for eight complex traits. Genomic predictions were performed using either data from a standard marker array or variants preselected from WGS based on association tests.

Results

The accuracies of genomic predictions based on preselected WGS variants were not robust across traits and lines and the improvements in prediction accuracy that we achieved so far with WGS compared to standard marker arrays were generally small. The most favourable results for WGS were obtained when the largest training sets were available and standard marker arrays were augmented with preselected variants with statistically significant associations to the trait. With this method and training sets of around 80k individuals, the accuracy of within-line genomic predictions was on average improved by 0.025. With multi-line training sets, improvements of 0.04 compared to marker arrays could be expected.

Conclusions

Our results showed that WGS has limited potential to improve the accuracy of genomic predictions compared to marker arrays in intensely selected pig lines. Thus, although we expect that larger improvements in accuracy from the use of WGS are possible with a combination of larger training sets and optimised pipelines for generating and analysing such datasets, the use of WGS in the current implementations of genomic prediction should be carefully evaluated against the cost of large-scale WGS data on a case-by-case basis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00756-0.

The Pioneer Advantage: Filling the blank spots on the map of genome diversity in Europe

Documenting genome diversity is important for the local biomedical communities and instrumental in developing precision and personalized medicine. Currently, tens of thousands of whole-genome sequences from Europe are publicly available, but most of these represent populations of developed countries of Europe. The uneven distribution of the available data is further impaired by the lack of data sharing. Recent whole-genome studies in Eastern Europe, one in Ukraine and one in Russia, demonstrated that local genome diversity and population structure from Eastern Europe historically had not been fully represented. An unexpected wealth of genomic variation uncovered in these studies was not so much a consequence of high variation within their population, but rather due to the "pioneer advantage." We discovered more variants because we were the first to prospect in the Eastern European genome pool. This simple comparison underscores the importance of removing the remaining geographic genome deserts from the rest of the world map of the human genome diversity.

HLA alleles, disease severity, and age associate with T-cell responses following infection with SARS-CoV-2

Memory T-cell responses following SARS-CoV-2 infection have been extensively investigated but many studies have been small with a limited range of disease severity. Here we analyze SARS-CoV-2 reactive T-cell responses in 768 convalescent SARS-CoV-2-infected (cases) and 500 uninfected (controls) Icelanders. The T-cell responses are stable three to eight months after SARS-CoV-2 infection, irrespective of disease severity and even those with the mildest symptoms induce broad and persistent T-cell responses. Robust CD4⁺ T-cell responses are detected against all measured proteins (M, N, S and S1) while the N protein induces strongest CD8⁺ T-cell responses. CD4⁺ T-cell responses correlate with disease severity, humoral responses and age, whereas CD8⁺ T-cell responses correlate with age and functional antibodies. Further, CD8⁺ T-cell responses associate with several class I HLA alleles. Our results, provide new insight into HLA restriction of CD8⁺ T-cell immunity and other factors contributing to heterogeneity of T-cell responses following SARS-CoV-2 infection.

A study of 768 convalescent SARS CoV-2-infected and 500 uninfected Icelanders reveals broad and stable T-cell responses 3-8 months from infection. HLA alleles, disease severity, and age contribute to the heterogeneity of cellular immunity.

Association Study between Polymorphisms in DNA Methylation-Related Genes and Testicular Germ Cell Tumor Risk

BACKGROUND

Testicular germ cell tumors (TGCT), histologically classified as seminomas and nonseminomas, are believed to arise from primordial gonocytes, with the maturation process blocked when they are subjected to DNA methylation reprogramming. SNPs in DNA methylation machinery and folate-dependent one-carbon metabolism genes have been postulated to influence the proper establishment of DNA methylation.

METHODS

In this pathway-focused investigation, we evaluated the association between 273 selected tag SNPs from 28 DNA methylation-related genes and TGCT risk. We carried out association analysis at individual SNP and gene-based level using summary statistics from the Genome Wide Association Study meta-analysis recently conducted by the international Testicular Cancer Consortium on 10,156 TGCT cases and 179,683 controls.

RESULTS

In individual SNP analyses, seven SNPs, four mapping within MTHFR, were associated with TGCT risk after correction for multiple testing (q ≤ 0.05). Queries of public databases showed that three of these SNPs were associated with MTHFR changes in enzymatic activity (rs1801133) or expression level in testis tissue (rs12121543, rs1476413). Gene-based analyses revealed MTHFR (q = 8.4 × 10-4), methyl-CpG-binding protein 2 (MECP2; q = 2 × 10-3), and ZBTB4 (q = 0.03) as the top TGCT-associated genes. Stratifying by tumor histology, four MTHFR SNPs were associated with seminoma. In gene-based analysis MTHFR was associated with risk of seminoma (q = 2.8 × 10-4), but not with nonseminomatous tumors (q = 0.22).

CONCLUSIONS

Genetic variants within MTHFR, potentially having an impact on the DNA methylation pattern, are associated with TGCT risk.

IMPACT

This finding suggests that TGCT pathogenesis could be associated with the folate cycle status, and this relation could be partly due to hereditary factors.

Extent to which array genotyping and imputation with large reference panels approximate deep whole-genome sequencing

Understanding the genetic basis of human diseases and traits is dependent on the identification and accurate genotyping of genetic variants. Deep whole-genome sequencing (WGS), the gold standard technology for SNP and indel identification and genotyping, remains very expensive for most large studies. Here, we quantify the extent to which array genotyping followed by genotype imputation can approximate WGS in studies of individuals of African, Hispanic/Latino, and European ancestry in the US and of Finnish ancestry in Finland (a population isolate). For each study, we performed genotype imputation by using the genetic variants present on the Illumina Core, OmniExpress, MEGA, and Omni 2.5M arrays with the 1000G, HRC, and TOPMed imputation reference panels. Using the Omni 2.5M array and the TOPMed panel, ≥90% of bi-allelic single-nucleotide variants (SNVs) are well imputed (r2 > 0.8) down to minor-allele frequencies (MAFs) of 0.14% in African, 0.11% in Hispanic/Latino, 0.35% in European, and 0.85% in Finnish ancestries. There was little difference in TOPMed-based imputation quality among the arrays with >700k variants. Individual-level imputation quality varied widely between and within the three US studies. Imputation quality also varied across genomic regions, producing regions where even common (MAF > 5%) variants were consistently not well imputed across ancestries. The extent to which array genotyping and imputation can approximate WGS therefore depends on reference panel, genotype array, sample ancestry, and genomic location. Imputation quality by variant or genomic region can be queried with our new tool, RsqBrowser, now deployed on the Michigan Imputation Server.

Comparative Analyses of Turkish Variome and Widely Used Genomic Variation Databases for the Evaluation of Rare Sequence Variants in Turkish Individuals: Idiopathic Hypogonadotropic Hypogonadism as a Disease Model

Objective

With the increasing use of whole-exome sequencing, one of the challenges in identifying the causal allele for a Mendelian disease is the lack of availability of population-specific human genetic variation reference databases. The people of Turkey were not represented in GnomAD or other publicly available large databases until recently, when the first comprehensive genomic variation database, Turkish Variome (TRV), was published. The aim of this study was to evaluate whether TRV or other publicly available large genomic variation databases can reliably be used for rare disease variant evaluation in Turkish individuals.

Methods

Sixty non-disease-causing, non-synonymous variants (minor allele frequencies >1%) were identified in 58 genes that are known to be associated with idiopathic hypogonadotropic hypogonadism from a large Turkish patient cohort. The allelic frequencies of these variants were then compared with those in various public genomic variation databases, including TRV.

Results

Our cohort variants showed the highest correlations with those in the TRV, Iranome, and The Greater Middle East Variome, in decreasing order.

Conclusion

These results suggest that the TRV is the appropriate database to use for rare genomic variant evaluations in the Turkish population. Our data also suggest that variomes from geographic neighborhoods may serve as substitute references for populations devoid of their own genomic variation databases.

GAWMerge expands GWAS sample size and diversity by combining array-based genotyping and whole-genome sequencing

Genome-wide association studies (GWAS) have made impactful discoveries for complex diseases, often by amassing very large sample sizes. Yet, GWAS of many diseases remain underpowered, especially for non-European ancestries. One cost-effective approach to increase sample size is to combine existing cohorts, which may have limited sample size or be case-only, with public controls, but this approach is limited by the need for a large overlap in variants across genotyping arrays and the scarcity of non-European controls. We developed and validated a protocol, Genotyping Array-WGS Merge (GAWMerge), for combining genotypes from arrays and whole-genome sequencing, ensuring complete variant overlap, and allowing for diverse samples like Trans-Omics for Precision Medicine to be used. Our protocol involves phasing, imputation, and filtering. We illustrated its ability to control technology driven artifacts and type-I error, as well as recover known disease-associated signals across technologies, independent datasets, and ancestries in smoking-related cohorts. GAWMerge enables genetic studies to leverage existing cohorts to validly increase sample size and enhance discovery for understudied traits and ancestries.

Multiomics analysis of rheumatoid arthritis yields sequence variants that have large effects on risk of the seropositive subset

OBJECTIVES

To find causal genes for rheumatoid arthritis (RA) and its seropositive (RF and/or ACPA positive) and seronegative subsets.

METHODS

We performed a genome-wide association study (GWAS) of 31 313 RA cases (68% seropositive) and ~1 million controls from Northwestern Europe. We searched for causal genes outside the HLA-locus through effect on coding, mRNA expression in several tissues and/or levels of plasma proteins (SomaScan) and did network analysis (Qiagen).

RESULTS

We found 25 sequence variants for RA overall, 33 for seropositive and 2 for seronegative RA, altogether 37 sequence variants at 34 non-HLA loci, of which 15 are novel. Genomic, transcriptomic and proteomic analysis of these yielded 25 causal genes in seropositive RA and additional two overall. Most encode proteins in the network of interferon-alpha/beta and IL-12/23 that signal through the JAK/STAT-pathway. Highlighting those with largest effect on seropositive RA, a rare missense variant in STAT4 (rs140675301-A) that is independent of reported non-coding STAT4-variants, increases the risk of seropositive RA 2.27-fold (p=2.1×10-9), more than the rs2476601-A missense variant in PTPN22 (OR=1.59, p=1.3×10-160). STAT4 rs140675301-A replaces hydrophilic glutamic acid with hydrophobic valine (Glu128Val) in a conserved, surface-exposed loop. A stop-mutation (rs76428106-C) in FLT3 increases seropositive RA risk (OR=1.35, p=6.6×10-11). Independent missense variants in TYK2 (rs34536443-C, rs12720356-C, rs35018800-A, latter two novel) associate with decreased risk of seropositive RA (ORs=0.63-0.87, p=10-9-10-27) and decreased plasma levels of interferon-alpha/beta receptor 1 that signals through TYK2/JAK1/STAT4.

CONCLUSION

Sequence variants pointing to causal genes in the JAK/STAT pathway have largest effect on seropositive RA, while associations with seronegative RA remain scarce.

CMDB: the comprehensive population genome variation database of China

A high-quality genome variation database derived from a large-scale population is one of the most important infrastructures for genomics, clinical and translational medicine research. Here, we developed the Chinese Millionome Database (CMDB), a database that contains 9.04 million single nucleotide variants (SNV) with allele frequency information derived from low-coverage (0.06×-0.1×) whole-genome sequencing (WGS) data of 141 431 unrelated healthy Chinese individuals. These individuals were recruited from 31 out of the 34 administrative divisions in China, covering Han and 36 other ethnic minorities. CMDB, housing the WGS data of a multi-ethnic Chinese population featuring wide geographical distribution, has become the most representative and comprehensive Chinese population genome database to date. Researchers can quickly search for variant, gene or genomic regions to obtain the variant information, including mutation basic information, allele frequency, genic annotation and overview of frequencies in global populations. Furthermore, the CMDB also provides information on the association of the variants with a range of phenotypes, including height, BMI, maternal age and twin pregnancy. Based on these data, researchers can conduct meta-analysis of related phenotypes. CMDB is freely available at https://db.cngb.org/cmdb/.

New Drug Development and Clinical Trial Design by Applying Genomic Information Management

Depending on the patients' genotype, the same drug may have different efficacies or side effects. With the cost of genomic analysis decreasing and reliability of analysis methods improving, vast amount of genomic information has been made available. Several studies in pharmacology have been based on genomic information to select the optimal drug, determine the dose, predict efficacy, and prevent side effects. This paper reviews the tissue specificity and genomic information of cancer. If the tissue specificity of cancer is low, cancer is induced in various organs based on a single gene mutation. Basket trials can be performed for carcinomas with low tissue specificity, confirming the efficacy of one drug for a single gene mutation in various carcinomas. Conversely, if the tissue specificity of cancer is high, cancer is induced in only one organ based on a single gene mutation. An umbrella trial can be performed for carcinomas with a high tissue specificity. Some drugs are effective for patients with a specific genotype. A companion diagnostic strategy that prescribes a specific drug for patients selected with a specific genotype is also reviewed. Genomic information is used in pharmacometrics to identify the relationship among pharmacokinetics, pharmacodynamics, and biomarkers of disease treatment effects. Utilizing genomic information, sophisticated clinical trials can be designed that will be better suited to the patients of specific genotypes. Genomic information also provides prospects for innovative drug development. Through proper genomic information management, factors relating to drug response and effects can be determined by selecting the appropriate data for analysis and by understanding the structure of the data. Selecting pre-processing and appropriate machine-learning libraries for use as machine-learning input features is also necessary. Professional curation of the output result is also required. Personalized medicine can be realized using a genome-based customized clinical trial design.

Genetics in Familial Intrahepatic Cholestasis: Clinical Patterns and Development of Liver and Biliary Cancers: A Review of the Literature

The family of inherited intrahepatic cholestasis includes autosomal recessive cholestatic rare diseases of childhood involved in bile acids secretion or bile transport defects. Specific genetic pathways potentially cause many otherwise unexplained cholestasis or hepatobiliary tumours in a healthy liver. Lately, next-generation sequencing and whole-exome sequencing have improved the diagnostic procedures of familial intrahepatic cholestasis (FIC), as well as the discovery of several genes responsible for FIC. Moreover, mutations in these genes, even in the heterozygous status, may be responsible for cryptogenic cholestasis in both young and adults. Mutations in FIC genes can influence serum and hepatic levels of bile acids. Experimental studies on the NR1H4 gene have shown that high bile acids concentrations cause excessive production of inflammatory cytokines, resistance to apoptosis, and increased cell regeneration, all risk conditions for developing hepatocellular carcinoma (HCC) and cholangiocarcinoma (CCA). NR1H4 gene encodes farnesoid X-activated receptor having a pivotal role in bile salts synthesis. Moreover, HCC and CCA can emerge in patients with several FIC genes such as ABCB11, ABCB4 and TJP2. Herein, we reviewed the available data on FIC-related hepatobiliary cancers, reporting on genetics to the pathophysiology, the risk factors and the clinical presentation.

Whole-Genome Sequencing and Genomic Variant Analysis of Kazakh Individuals

Kazakhstan, the ninth-largest country in the world, is located along the Great Silk Road and connects Europe with Asia. Historically, its territory has been inhabited by nomadic tribes, and modern-day Kazakhstan is a multiethnic country with a dominant Kazakh population. We sequenced and analyzed the genomes of five ethnic Kazakhs at high coverage using the Illumina HiSeq2000 next-generation sequencing platform. The five Kazakhs yielded a total number of base pairs ranging from 87,308,581,400 to 107,526,741,301. On average, 99.06% were properly mapped. Based on the Het/Hom and Ti/Tv ratios, the quality of the genomic data ranged from 1.35 to 1.49 and from 2.07 to 2.08, respectively. Genetic variants were identified and annotated. Functional analysis of the genetic variants identified several variants that were associated with higher risks of metabolic and neurogenerative diseases. The present study showed high levels of genetic admixture of Kazakhs that were comparable to those of other Central Asians. These whole-genome sequence data of healthy Kazakhs could contribute significantly to biomedical studies of common diseases as their findings could allow better insight into the genotype–phenotype relations at the population level.

The sequences of 150,119 genomes in the UK Biobank

Detailed knowledge of how diversity in the sequence of the human genome affects phenotypic diversity depends on a comprehensive and reliable characterization of both sequences and phenotypic variation. Over the past decade, insights into this relationship have been obtained from whole-exome sequencing or whole-genome sequencing of large cohorts with rich phenotypic data^1,2. Here we describe the analysis of whole-genome sequencing of 150,119 individuals from the UK Biobank³. This constitutes a set of high-quality variants, including 585,040,410 single-nucleotide polymorphisms, representing 7.0% of all possible human single-nucleotide polymorphisms, and 58,707,036 indels. This large set of variants allows us to characterize selection based on sequence variation within a population through a depletion rank score of windows along the genome. Depletion rank analysis shows that coding exons represent a small fraction of regions in the genome subject to strong sequence conservation. We define three cohorts within the UK Biobank: a large British Irish cohort, a smaller African cohort and a South Asian cohort. A haplotype reference panel is provided that allows reliable imputation of most variants carried by three or more sequenced individuals. We identified 895,055 structural variants and 2,536,688 microsatellites, groups of variants typically excluded from large-scale whole-genome sequencing studies. Using this formidable new resource, we provide several examples of trait associations for rare variants with large effects not found previously through studies based on whole-exome sequencing and/or imputation.

Germline De Novo Mutations as a Cause of Childhood Cancer

Germline de novo mutations (DNMs) represent one of the important topics that need extensive attention from epidemiologists, geneticists, and other relevant stakeholders. Advances in next-generation sequencing technologies allowed examination of parent-offspring trios to ascertain the frequency of germline DNMs. Many epidemiological risk factors for childhood cancer are indicative of DNMs as a mechanism. The aim of this review was to give an overview of germline DNMs, their causes in general, and to discuss their relation to childhood cancer risk. In addition, we highlighted existing gaps in knowledge in many topics of germline DNMs in childhood cancer that need exploration and collaborative efforts.

Advances and Trends in Omics Technology Development

The human history has witnessed the rapid development of technologies such as high-throughput sequencing and mass spectrometry that led to the concept of “omics” and methodological advancement in systematically interrogating a cellular system. Yet, the ever-growing types of molecules and regulatory mechanisms being discovered have been persistently transforming our understandings on the cellular machinery. This renders cell omics seemingly, like the universe, expand with no limit and our goal toward the complete harness of the cellular system merely impossible. Therefore, it is imperative to review what has been done and is being done to predict what can be done toward the translation of omics information to disease control with minimal cell perturbation. With a focus on the “four big omics,” i.e., genomics, transcriptomics, proteomics, metabolomics, we delineate hierarchies of these omics together with their epiomics and interactomics, and review technologies developed for interrogation. We predict, among others, redoxomics as an emerging omics layer that views cell decision toward the physiological or pathological state as a fine-tuned redox balance.

Germline variants of ATG7 in familial cholangiocarcinoma alter autophagy and p62

Autophagy is a housekeeping mechanism tasked with eliminating misfolded proteins and damaged organelles to maintain cellular homeostasis. Autophagy deficiency results in increased oxidative stress, DNA damage and chronic cellular injury. Among the core genes in the autophagy machinery, ATG7 is required for autophagy initiation and autophagosome formation. Based on the analysis of an extended pedigree of familial cholangiocarcinoma, we determined that all affected family members had a novel germline mutation (c.2000C>T p.Arg659* (p.R659*)) in ATG7. Somatic deletions of ATG7 were identified in the tumors of affected individuals. We applied linked-read sequencing to one tumor sample and demonstrated that the ATG7 somatic deletion and germline mutation were located on distinct alleles, resulting in two hits to ATG7. From a parallel population genetic study, we identified a germline polymorphism of ATG7 (c.1591C>G p.Asp522Glu (p.D522E)) associated with increased risk of cholangiocarcinoma. To characterize the impact of these germline ATG7 variants on autophagy activity, we developed an ATG7-null cell line derived from the human bile duct. The mutant p.R659* ATG7 protein lacked the ability to lipidate its LC3 substrate, leading to complete loss of autophagy and increased p62 levels. Our findings indicate that germline ATG7 variants have the potential to impact autophagy function with implications for cholangiocarcinoma development.

Genetic Analysis of ABCB4 Mutations and Variants Related to the Pathogenesis and Pathophysiology of Low Phospholipid-Associated Cholelithiasis

Clinical studies have revealed that the ABCB4 gene encodes the phospholipid transporter on the canalicular membrane of hepatocytes, and its mutations and variants are the genetic basis of low phospholipid-associated cholelithiasis (LPAC), a rare type of gallstone disease caused by a single-gene mutation or variation. The main features of LPAC include a reduction or deficiency of phospholipids in bile, symptomatic cholelithiasis at <40 years of age, intrahepatic sludge and microlithiasis, mild chronic cholestasis, a high cholesterol/phospholipid ratio in bile, and recurrence of biliary symptoms after cholecystectomy. Needle-like cholesterol crystals, putatively “anhydrous” cholesterol crystallization at low phospholipid concentrations in model and native bile, are characterized in ABCB4 knockout mice, a unique animal model for LPAC. Gallbladder bile with only trace amounts of phospholipids in these mice is supersaturated with cholesterol, with lipid composition plotting in the left two-phase zone of the ternary phase diagram, consistent with “anhydrous” cholesterol crystallization. In this review, we summarize the molecular biology and physiological functions of ABCB4 and comprehensively discuss the latest advances in the genetic analysis of ABCB4 mutations and variations and their roles in the pathogenesis and pathophysiology of LPAC in humans, based on the results from clinical studies and mouse experiments. To date, approximately 158 distinct LPAC-causing ABCB4 mutations and variants in humans have been reported in the literature, indicating that it is a monogenic risk factor for LPAC. The elucidation of the ABCB4 function in the liver, the identification of ABCB4 mutations and variants in LPAC patients, and the exploration of gene therapy for ABCB4 deficiency in animal models can help us to better understand the cellular, molecular, and genetic mechanisms underlying the onset of the disease, and will pave the way for early diagnosis and prevention of susceptible subjects and effective intervention for LPAC in patients.

Rare and population-specific functional variation across pig lines

BACKGROUND

It is expected that functional, mainly missense and loss-of-function (LOF), and regulatory variants are responsible for most phenotypic differences between breeds and genetic lines of livestock species that have undergone diverse selection histories. However, there is still limited knowledge about the existing missense and LOF variation in commercial livestock populations, in particular regarding population-specific variation and how it can affect applications such as across-breed genomic prediction.

METHODS

We re-sequenced the whole genome of 7848 individuals from nine commercial pig lines (average sequencing coverage: 4.1×) and imputed whole-genome genotypes for 440,610 pedigree-related individuals. The called variants were categorized according to predicted functional annotation (from LOF to intergenic) and prevalence level (number of lines in which the variant segregated; from private to widespread). Variants in each category were examined in terms of their distribution along the genome, alternative allele frequency, per-site Wright's fixation index (F_ST), individual load, and association to production traits.

RESULTS

Of the 46 million called variants, 28% were private (called in only one line) and 21% were widespread (called in all nine lines). Genomic regions with a low recombination rate were enriched with private variants. Low-prevalence variants (called in one or a few lines only) were enriched for lower allele frequencies, lower F_ST, and putatively functional and regulatory roles (including LOF and deleterious missense variants). On average, individuals carried fewer private deleterious missense alleles than expected compared to alleles with other predicted consequences. Only a small subset of the low-prevalence variants had intermediate allele frequencies and explained small fractions of phenotypic variance (up to 3.2%) of production traits. The significant low-prevalence variants had higher per-site F_ST than the non-significant ones. These associated low-prevalence variants were tagged by other more widespread variants in high linkage disequilibrium, including intergenic variants.

CONCLUSIONS

Most low-prevalence variants have low minor allele frequencies and only a small subset of low-prevalence variants contributed detectable fractions of phenotypic variance of production traits. Accounting for low-prevalence variants is therefore unlikely to noticeably benefit across-breed analyses, such as the prediction of genomic breeding values in a population using reference populations of a different genetic background.

Genetic architecture of band neutrophil fraction in Iceland

The characteristic lobulated nuclear morphology of granulocytes is partially determined by composition of nuclear envelope proteins. Abnormal nuclear morphology is primarily observed as an increased number of hypolobulated immature neutrophils, called band cells, during infection or in rare envelopathies like Pelger-Huët anomaly. To search for sequence variants affecting nuclear morphology of granulocytes, we performed a genome-wide association study using band neutrophil fraction from 88,101 Icelanders. We describe 13 sequence variants affecting band neutrophil fraction at nine loci. Five of the variants are at the Lamin B receptor (LBR) locus, encoding an inner nuclear membrane protein. Mutations in LBR are linked to Pelger-Huët anomaly. In addition, we identify cosegregation of a rare stop-gain sequence variant in LBR and Pelger Huët anomaly in an Icelandic eight generation pedigree, initially reported in 1963. Two of the other loci include genes which, like LBR, play a role in the nuclear membrane function and integrity. These GWAS results highlight the role proteins of the inner nuclear membrane have as important for neutrophil nuclear morphology.

Genomic analyses of 10,376 individuals in the Westlake BioBank for Chinese (WBBC) pilot project

We initiate the Westlake BioBank for Chinese (WBBC) pilot project with 4,535 whole-genome sequencing (WGS) individuals and 5,841 high-density genotyping individuals, and identify 81.5 million SNPs and INDELs, of which 38.5% are absent in dbSNP Build 151. We provide a population-specific reference panel and an online imputation server ( https://wbbc.westlake.edu.cn/ ) which could yield substantial improvement of imputation performance in Chinese population, especially for low-frequency and rare variants. By analyzing the singleton density of the WGS data, we find selection signatures in SNX29, DNAH1 and WDR1 genes, and the derived alleles of the alcohol metabolism genes (ADH1A and ADH1B) emerge around 7,000 years ago and tend to be more common from 4,000 years ago in East Asia. Genetic evidence supports the corresponding geographical boundaries of the Qinling-Huaihe Line and Nanling Mountains, which separate the Han Chinese into subgroups, and we reveal that North Han was more homogeneous than South Han.

Genetics of atrial fibrillation-an update of recent findings

Atrial fibrillation (AF) is a common cardiac arrhythmia and a major risk factor for stroke, heart failure, and premature death. AF has a strong genetic predisposition. This review highlights the recent findings on the genetics of AF from genome-wide association studies (GWAS) and high-throughput sequencing studies. The consensus from GWAS implies that AF is both polygenic and pleiotropic in nature. With the advent of whole-genome sequencing and whole-exome sequencing, rare variants associated with AF pathogenesis have been identified. The recent studies have contributed towards better understanding of AF pathogenesis.

Genome-wide analysis identifies gallstone-susceptibility loci including genes regulating gastrointestinal motility

BACKGROUND AND AIMS

Genome-wide association studies (GWAS) have identified several risk loci for gallstone disease. As with most polygenic traits, it is likely that many genetic determinants are undiscovered. The aim of this study was to identify genetic variants that represent new targets for gallstone research and treatment.

APPROACH AND RESULTS

We performed a GWAS of 28,627 gallstone cases and 348,373 controls in the UK Biobank, replicated findings in a Scottish cohort (1089 cases, 5228 controls), and conducted a GWA meta-analysis (43,639 cases, 506,798 controls) with the FinnGen cohort. We assessed pathway enrichment using gene-based then gene-set analysis and tissue expression of identified genes in Genotype-Tissue Expression project data. We constructed a polygenic risk score (PRS) and evaluated phenotypic traits associated with the score. Seventy-five risk loci were identified (p < 5 × 10-8 ), of which 46 were new. Pathway enrichment revealed associations with lipid homeostasis, glucuronidation, phospholipid metabolism, and gastrointestinal motility. Anoctamin 1 (ANO1) and transmembrane Protein 147 (TMEM147), both in novel, replicated loci, are expressed in the gallbladder and gastrointestinal tract. Both regulate gastrointestinal motility. The gallstone risk allele rs7599-A leads to suppression of hepatic TMEM147 expression, suggesting that the protein protects against gallstone formation. The highest decile of the PRS demonstrated a 6-fold increased odds of gallstones compared with the lowest decile. The PRS was strongly associated with increased body mass index, serum liver enzymes, and C-reactive protein concentrations, and decreased lipoprotein cholesterol concentrations.

CONCLUSIONS

This GWAS demonstrates the polygenic nature of gallstone risk and identifies 46 novel susceptibility loci. We implicate genes influencing gastrointestinal motility in the pathogenesis of gallstones.

Universal tumor screening in a population with MSH6- and PMS2-associated Lynch syndrome

PURPOSE

Universal screening for Lynch syndrome (LS) on resected colorectal carcinomas (CRCs) and endometrial carcinomas (ECs) was implemented in Iceland in 2017 using immunohistochemistry (IHC) for mismatch repair (MMR) proteins. We examined the efficacy of the universal screening algorithm to detect LS and the diagnostic accuracy of MMR IHC by comparing results with a population-based genotype database.

METHODS

All patients diagnosed with CRC or EC per the Icelandic Cancer Registry from 2017 to 2019 who had tumor MMR IHC performed were included. Pathology reports and patient charts were reviewed. MMR IHC stains were crossmatched with genotyping results obtained from the deCODE database.

RESULTS

IHC staining was done on 404 patients with CRC and 74 patients with EC. A total of 61 (15.1%) patients with CRC and 15 (20.3%) patients with EC were MMR-deficient. MMR IHC had 88.9% sensitivity in identifying patients with LS and a positive predictive value of 10.7%. Only 50% of individuals were appropriately referred for genetic testing, leading to underdiagnosis of LS.

CONCLUSION

Universal screening for LS using MMR protein IHC in CRC and EC accurately identified patients appropriate for genetic testing in a population with MSH6 and PMS2 LS predominance. Because of lack of referral to genetic counseling, only 50% of patients with LS were identified through the screening algorithm.