Rare and common variants: twenty arguments

Discovery of Resistance Genes in Rye by Targeted Long-Read Sequencing and Association Genetics

The majority of released rye cultivars are susceptible to leaf rust because of a low level of resistance in the predominant hybrid rye-breeding gene pools Petkus and Carsten. To discover new sources of leaf rust resistance, we phenotyped a diverse panel of inbred lines from the less prevalent Gülzow germplasm using six distinct isolates of Puccinia recondita f. sp. secalis and found that 55 out of 92 lines were resistant to all isolates. By performing a genome-wide association study using 261,406 informative SNP markers, we identified five resistance-associated QTLs on chromosome arms 1RS, 1RL, 2RL, 5RL and 7RS. To identify candidate Puccinia recondita (Pr) resistance genes in these QTLs, we sequenced the rye nucleotide-binding leucine-rich repeat (NLR) intracellular immune receptor complement using a Triticeae NLR bait-library and PacBio^® long-read single-molecule high-fidelity (HiFi) sequencing. Trait-genotype correlations across 10 resistant and 10 susceptible lines identified four candidate NLR-encoding Pr genes. One of these physically co-localized with molecular markers delimiting Pr3 on chromosome arm 1RS and the top-most resistance-associated QTL in the panel.

Gene Region Association Analysis of Longitudinal Quantitative Traits Based on a Function-On-Function Regression Model

In the process of growth and development in life, gene expressions that control quantitative traits will turn on or off with time. Studies of longitudinal traits are of great significance in revealing the genetic mechanism of biological development. With the development of ultra-high-density sequencing technology, the associated analysis has tremendous challenges to statistical methods. In this paper, a longitudinal functional data association test (LFDAT) method is proposed based on the function-on-function regression model. LFDAT can simultaneously treat phenotypic traits and marker information as continuum variables and analyze the association of longitudinal quantitative traits and gene regions. Simulation studies showed that: 1) LFDAT performs well for both linkage equilibrium simulation and linkage disequilibrium simulation, 2) LFDAT has better performance for gene regions (include common variants, low-frequency variants, rare variants and mixture), and 3) LFDAT can accurately identify gene switching in the growth and development stage. The longitudinal data of the Oryza sativa projected shoot area is analyzed by LFDAT. It showed that there is the advantage of quick calculations. Further, an association analysis was conducted between longitudinal traits and gene regions by integrating the micro effects of multiple related variants and using the information of the entire gene region. LFDAT provides a feasible method for studying the formation and expression of longitudinal traits.

Simulation Research on the Methods of Multi-Gene Region Association Analysis Based on a Functional Linear Model

Genome-wide association analysis is an important approach to identify genetic variants associated with complex traits. Complex traits are not only affected by single gene loci, but also by the interaction of multiple gene loci. Studies of association between gene regions and quantitative traits are of great significance in revealing the genetic mechanism of biological development. There have been a lot of studies on single-gene region association analysis, but the application of functional linear models in multi-gene region association analysis is still less. In this paper, a functional multi-gene region association analysis test method is proposed based on the functional linear model. From the three directions of common multi-gene region method, multi-gene region weighted method and multi-gene region loci weighted method, that test method is studied combined with computer simulation. The following conclusions are obtained through computer simulation: (a) The functional multi-gene region association analysis test method has higher power than the functional single gene region association analysis test method; (b) The functional multi-gene region weighted method performs better than the common functional multi-gene region method; (c) the functional multi-gene region loci weighted method is the best method for association analysis on three directions of the common multi-gene region method; (d) the performance of the Step method and Multi-gene region loci weighted Step for multi-gene regions is the best in general. Functional multi-gene region association analysis test method can theoretically provide a feasible method for the study of complex traits affected by multiple genes.

Whole-Exome Sequencing Implicates the USP34 rs777591A > G Intron Variant in Chronic Obstructive Pulmonary Disease in a Kashi Cohort

Genetic factors are important factors in chronic obstructive pulmonary disease (COPD) onset. Plenty of risk and new causative genes for COPD have been identified in patients of the Chinese Han population. In contrast, we know considerably little concerning the genetics in the Kashi COPD population (Uyghur). This study aims at clarifying the genetic maps regarding COPD susceptibility in Kashi (China). Whole-exome sequencing (WES) was used to analyze three Uyghur families with COPD in Kashi (eight patients and one healthy control). Sanger sequencing was also used to verify the WES results in 541 unrelated Uyghur COPD patients and 534 Uyghur healthy controls. WES showed 72 single nucleotide variants (SNVs), two deletions, and small insertions (InDels), 26 copy number variants (CNVs), and 34 structural variants (SVs), including g.71230620T > A (rs12449210T > A, NC_000,016.10) in the HYDIN axonemal central pair apparatus protein (HYDIN) gene and g.61190482A > G (rs777591A > G, NC_000002.12) in the ubiquitin-specific protease 34 (USP34) gene. After Sanger sequencing, we found that rs777591“AA” under different genetic models except for the dominant model (adjusted OR = 0.8559, 95%CI 0.6568–1.115, p > .05), could significantly reduce COPD risk, but rs12449210T > A was not related to COPD. In stratified analysis of smoking status, rs777591“AA” reduced COPD risk significantly among the nonsmoker group. Protein and mRNA expression of USP34 in cigarette smoke extract-treated BEAS-2b cells increased significantly compared with those in the control group. Our findings associate the USP34 rs777591“AA” genotype as a protector factor in COPD.

Definitive roles of TOMM40-APOE-APOC1 variants in the Alzheimer's risk

Despite advances, the roles of genetic variants from the APOE-harboring 19q13.32 region in Alzheimer's disease (AD) remain controversial. We leverage a comprehensive approach to gain insights into a more homogeneous genetic architecture of AD in this region. We use a sample of 2,673 AD-affected and 16,246 unaffected subjects from 4 studies and validate our main findings in the landmark Alzheimer's Disease Genetics Consortium cohort (3,662 AD-cases and 1,541 controls). We report the remarkably high excesses of the AD risk for carriers of the ε4 allele who also carry minor alleles of rs2075650 (TOMM40) and rs12721046 (APOC1) polymorphisms compared to carriers of their major alleles. The exceptionally high 4.37-fold (p=1.34 × 10^-3) excess was particularly identified for the minor allele homozygotes. The beneficial and adverse variants were significantly depleted and enriched, respectively, in the AD-affected families. This study provides compelling evidence for the definitive roles of the APOE-TOMM40-APOC1 variants in the AD risk.

Rare variant association tests for ancestry-matched case-control data based on conditional logistic regression

With the increasing volume of human sequencing data available, analysis incorporating external controls becomes a popular and cost-effective approach to boost statistical power in disease association studies. To prevent spurious association due to population stratification, it is important to match the ancestry backgrounds of cases and controls. However, rare variant association tests based on a standard logistic regression model are conservative when all ancestry-matched strata have the same case-control ratio and might become anti-conservative when case-control ratio varies across strata. Under the conditional logistic regression (CLR) model, we propose a weighted burden test (CLR-Burden), a variance component test (CLR-SKAT) and a hybrid test (CLR-MiST). We show that the CLR model coupled with ancestry matching is a general approach to control for population stratification, regardless of the spatial distribution of disease risks. Through extensive simulation studies, we demonstrate that the CLR-based tests robustly control type 1 errors under different matching schemes and are more powerful than the standard Burden, SKAT and MiST tests. Furthermore, because CLR-based tests allow for different case-control ratios across strata, a full-matching scheme can be employed to efficiently utilize all available cases and controls to accelerate the discovery of disease associated genes.

Constitutively active SARM1 variants that induce neuropathy are enriched in ALS patients

Background

In response to injury, neurons activate a program of organized axon self-destruction initiated by the NAD⁺ hydrolase, SARM1. In healthy neurons SARM1 is autoinhibited, but single amino acid changes can abolish autoinhibition leading to constitutively active SARM1 enzymes that promote degeneration when expressed in cultured neurons.

Methods

To investigate whether naturally occurring human variants might disrupt SARM1 autoinhibition and potentially contribute to risk for neurodegenerative disease, we assayed the enzymatic activity of all 42 rare SARM1 alleles identified among 8507 amyotrophic lateral sclerosis (ALS) patients and 9671 controls. We then intrathecally injected mice with virus expressing SARM1 constructs to test the capacity of an ALS-associated constitutively active SARM1 variant to promote neurodegeneration in vivo.

Results

Twelve out of 42 SARM1 missense variants or small in-frame deletions assayed exhibit constitutive NADase activity, including more than half of those that are unique to the ALS patients or that occur in multiple patients. There is a > 5-fold enrichment of constitutively active variants among patients compared to controls. Expression of constitutively active ALS-associated SARM1 alleles in cultured dorsal root ganglion (DRG) neurons is pro-degenerative and cytotoxic. Intrathecal injection of an AAV expressing the common SARM1 reference allele is innocuous to mice, but a construct harboring SARM1^V184G, the constitutively active variant found most frequently among the ALS patients, causes axon loss, motor dysfunction, and sustained neuroinflammation.

Conclusions

These results implicate rare hypermorphic SARM1 alleles as candidate genetic risk factors for ALS and other neurodegenerative conditions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13024-021-00511-x.

Quantitative Trait Loci Mapping of Adult Plant and Seedling Resistance to Stripe Rust (Puccinia striiformis Westend.) in a Multiparent Advanced Generation Intercross Wheat Population

Stripe rust caused by the biotrophic fungus Puccinia striiformis Westend. is one of the most important diseases of wheat worldwide, causing high yield and quality losses. Growing resistant cultivars is the most efficient way to control stripe rust, both economically and ecologically. Known resistance genes are already present in numerous cultivars worldwide. However, their effectiveness is limited to certain races within a rust population and the emergence of stripe rust races being virulent against common resistance genes forces the demand for new sources of resistance. Multiparent advanced generation intercross (MAGIC) populations have proven to be a powerful tool to carry out genetic studies on economically important traits. In this study, interval mapping was performed to map quantitative trait loci (QTL) for stripe rust resistance in the Bavarian MAGIC wheat population, comprising 394 F6 : 8 recombinant inbred lines (RILs). Phenotypic evaluation of the RILs was carried out for adult plant resistance in field trials at three locations across three years and for seedling resistance in a growth chamber. In total, 21 QTL for stripe rust resistance corresponding to 13 distinct chromosomal regions were detected, of which two may represent putatively new QTL located on wheat chromosomes 3D and 7D.

Domestication reshaped the genetic basis of inbreeding depression in a maize landrace compared to its wild relative, teosinte

Inbreeding depression is the reduction in fitness and vigor resulting from mating of close relatives observed in many plant and animal species. The extent to which the genetic load of mutations contributing to inbreeding depression is due to large-effect mutations versus variants with very small individual effects is unknown and may be affected by population history. We compared the effects of outcrossing and self-fertilization on 18 traits in a landrace population of maize, which underwent a population bottleneck during domestication, and a neighboring population of its wild relative teosinte. Inbreeding depression was greater in maize than teosinte for 15 of 18 traits, congruent with the greater segregating genetic load in the maize population that we predicted from sequence data. Parental breeding values were highly consistent between outcross and selfed offspring, indicating that additive effects determine most of the genetic value even in the presence of strong inbreeding depression. We developed a novel linkage scan to identify quantitative trait loci (QTL) representing large-effect rare variants carried by only a single parent, which were more important in teosinte than maize. Teosinte also carried more putative juvenile-acting lethal variants identified by segregation distortion. These results suggest a mixture of mostly polygenic, small-effect partially recessive effects in linkage disequilibrium underlying inbreeding depression, with an additional contribution from rare larger-effect variants that was more important in teosinte but depleted in maize following the domestication bottleneck. Purging associated with the maize domestication bottleneck may have selected against some large effect variants, but polygenic load is harder to purge and overall segregating mutational burden increased in maize compared to teosinte.

Inbreeding depression is the reduction in fitness and vigor resulting from mating of close relatives observed in many plant and animal species. Mating of close relatives increases the probability that an individual inherits two non-functioning mutations at the same gene, resulting in lower fitness of such matings. We do not know the extent to which inbreeding depression is due to mutations with large-effects versus small-effect polygenic variants. We compared the effects of outcrossing and self-fertilization on 18 traits in a landrace population of maize, which underwent a population bottleneck during domestication, and a neighboring population of its wild relative teosinte. Inbreeding depression was greater in maize than teosinte for 15 of 18 traits and we found that this was consistent with higher predicted ‘genetic load’ in maize based solely on the evolutionary conservation of the sequence variants observed in the population. We also mapped genome positions associated with inbreeding depression, identifying more and larger-effect genetic variants in teosinte than maize. These results suggest that during domestication, some of the rare large-effect variants in teosinte were bred out, but many genetic variants of small effects on inbreeding depression increased in frequency maize.

A comprehensive gene-centric pleiotropic association analysis for 14 psychiatric disorders with GWAS summary statistics

BACKGROUND

Recent genome-wide association studies (GWASs) have revealed the polygenic nature of psychiatric disorders and discovered a few of single-nucleotide polymorphisms (SNPs) associated with multiple psychiatric disorders. However, the extent and pattern of pleiotropy among distinct psychiatric disorders remain not completely clear.

METHODS

We analyzed 14 psychiatric disorders using summary statistics available from the largest GWASs by far. We first applied the cross-trait linkage disequilibrium score regression (LDSC) to estimate genetic correlation between disorders. Then, we performed a gene-based pleiotropy analysis by first aggregating a set of SNP-level associations into a single gene-level association signal using MAGMA. From a methodological perspective, we viewed the identification of pleiotropic associations across the entire genome as a high-dimensional problem of composite null hypothesis testing and utilized a novel method called PLACO for pleiotropy mapping. We ultimately implemented functional analysis for identified pleiotropic genes and used Mendelian randomization for detecting causal association between these disorders.

RESULTS

We confirmed extensive genetic correlation among psychiatric disorders, based on which these disorders can be grouped into three diverse categories. We detected a large number of pleiotropic genes including 5884 associations and 2424 unique genes and found that differentially expressed pleiotropic genes were significantly enriched in pancreas, liver, heart, and brain, and that the biological process of these genes was remarkably enriched in regulating neurodevelopment, neurogenesis, and neuron differentiation, offering substantial evidence supporting the validity of identified pleiotropic loci. We further demonstrated that among all the identified pleiotropic genes there were 342 unique ones linked with 6353 drugs with drug-gene interaction which can be classified into distinct types including inhibitor, agonist, blocker, antagonist, and modulator. We also revealed causal associations among psychiatric disorders, indicating that genetic overlap and causality commonly drove the observed co-existence of these disorders.

CONCLUSIONS

Our study is among the first large-scale effort to characterize gene-level pleiotropy among a greatly expanded set of psychiatric disorders and provides important insight into shared genetic etiology underlying these disorders. The findings would inform psychiatric nosology, identify potential neurobiological mechanisms predisposing to specific clinical presentations, and pave the way to effective drug targets for clinical treatment.

The correlation of substitution effects across populations and generations in the presence of nonadditive functional gene action

Allele substitution effects at quantitative trait loci (QTL) are part of the basis of quantitative genetics theory and applications such as association analysis and genomic prediction. In the presence of nonadditive functional gene action, substitution effects are not constant across populations. We develop an original approach to model the difference in substitution effects across populations as a first order Taylor series expansion from a "focal" population. This expansion involves the difference in allele frequencies and second-order statistical effects (additive by additive and dominance). The change in allele frequencies is a function of relationships (or genetic distances) across populations. As a result, it is possible to estimate the correlation of substitution effects across two populations using three elements: magnitudes of additive, dominance, and additive by additive variances; relationships (Nei's minimum distances or Fst indexes); and assumed heterozygosities. Similarly, the theory applies as well to distinct generations in a population, in which case the distance across generations is a function of increase of inbreeding. Simulation results confirmed our derivations. Slight biases were observed, depending on the nonadditive mechanism and the reference allele. Our derivations are useful to understand and forecast the possibility of prediction across populations and the similarity of GWAS effects.

Myeloid Toxicity of Radionuclide Cancer Therapy

Emergent genomic analytic techniques in patients with cancer offer the potential to define the risk of myelo dysplastic syndrome (MDS) and acute leukemia (AL) manifesting following targeted radionuclide therapy of metastatic lymphoma, neuroendocrine tumors (NETs), and prostate cancer. Characterization of the genetic profile will allow risk stratification of patients before theranostic radionuclide management of advanced cancers and offers the opportunity to minimize toxicity while preserving optimal individualized efficacy in the practice of personalized precision nuclear oncology. Our review of a single-center experience of prospective radionuclide theranostic management of metastatic non-Hodgkin lymphoma (NHL), NETs, and castration-resistant prostate cancer (metastatic castrate-resistant prostate cancer [mCRPC]) over the past decade, and comparison with published studies, shows that while the risk of significant myelotoxicity is generally low, at <3%, the consequences in the small minority of patients who develop MDS or AL are substantial, and survival is poor. Timely identification of patients at heightened risk of hematologic toxic complication, using novel genomic technology before institution of radionuclide therapy, will facilitate amelioration of myelotoxicity. In current clinical practice, the minimal hematological toxicity of chemo-free theranostic management of advanced cancer is significantly less compared with newly adopted chemotherapy -immunotherapy regimens, and the financial toxicity associated with these novel agents is avoided.

Germline breast cancer susceptibility genes, tumor characteristics, and survival

BACKGROUND

Mutations in certain genes are known to increase breast cancer risk. We study the relevance of rare protein-truncating variants (PTVs) that may result in loss-of-function in breast cancer susceptibility genes on tumor characteristics and survival in 8852 breast cancer patients of Asian descent.

METHODS

Gene panel sequencing was performed for 34 known or suspected breast cancer predisposition genes, of which nine genes (ATM, BRCA1, BRCA2, CHEK2, PALB2, BARD1, RAD51C, RAD51D, and TP53) were associated with breast cancer risk. Associations between PTV carriership in one or more genes and tumor characteristics were examined using multinomial logistic regression. Ten-year overall survival was estimated using Cox regression models in 6477 breast cancer patients after excluding older patients (≥75years) and stage 0 and IV disease.

RESULTS

PTV9genes carriership (n = 690) was significantly associated (p < 0.001) with more aggressive tumor characteristics including high grade (poorly vs well-differentiated, odds ratio [95% confidence interval] 3.48 [2.35-5.17], moderately vs well-differentiated 2.33 [1.56-3.49]), as well as luminal B [HER-] and triple-negative subtypes (vs luminal A 2.15 [1.58-2.92] and 2.85 [2.17-3.73], respectively), adjusted for age at diagnosis, study, and ethnicity. Associations with grade and luminal B [HER2-] subtype remained significant after excluding BRCA1/2 carriers. PTV25genes carriership (n = 289, excluding carriers of the nine genes associated with breast cancer) was not associated with tumor characteristics. However, PTV25genes carriership, but not PTV9genes carriership, was suggested to be associated with worse 10-year overall survival (hazard ratio [CI] 1.63 [1.16-2.28]).

CONCLUSIONS

PTV9genes carriership is associated with more aggressive tumors. Variants in other genes might be associated with the survival of breast cancer patients. The finding that PTV carriership is not just associated with higher breast cancer risk, but also more severe and fatal forms of the disease, suggests that genetic testing has the potential to provide additional health information and help healthy individuals make screening decisions.

Genetic analysis of four consanguineous multiplex families with inflammatory bowel disease

BACKGROUND

Family studies support a genetic predisposition to inflammatory bowel diseases (IBD), but known genetic variants only partially explain the disease heritability. Families with multiple affected individuals potentially harbour rare and high-impact causal variants. Long regions of homozygosity due to recent inbreeding may increase the risk of individuals bearing homozygous loss-of-function variants. This study aimed to identify rare and homozygous genetic variants contributing to IBD.

METHODS

Four families with known consanguinity and multiple cases of IBD were recruited. In a family-specific analysis, we utilised homozygosity mapping complemented by whole-exome sequencing.

RESULTS

We detected a single region of homozygosity shared by Crohn's disease cases from a family of Druze ancestry, spanning 2.6 Mb containing the NOD2 gene. Whole-exome sequencing did not identify any potentially damaging variants within the region, suggesting that non-coding variation may be involved. In addition, affected individuals in the families harboured several rare and potentially damaging homozygous variants in genes with a role in autophagy and innate immunity including LRRK1, WHAMM, DENND3, and C5.

CONCLUSION

This study examined the potential contribution of rare, high-impact homozygous variants in consanguineous families with IBD. While the analysis was not designed to achieve statistical significance, our findings highlight genes or loci that warrant further research. Non-coding variants affecting NOD2 may be of importance in Druze patients with Crohn's disease.

Clinical utility of polygenic risk scores for coronary artery disease

Over the past decade, substantial progress has been made in the discovery of alleles contributing to the risk of coronary artery disease. In addition to providing causal insights into disease, these endeavours have yielded and enabled the refinement of polygenic risk scores. These scores can be used to predict incident coronary artery disease in multiple cohorts and indicate the clinical response to some preventive therapies in post hoc analyses of clinical trials. These observations and the widespread ability to calculate polygenic risk scores from direct-to-consumer and health-care-associated biobanks have raised many questions about responsible clinical adoption. In this Review, we describe technical and downstream considerations for the derivation and validation of polygenic risk scores and current evidence for their efficacy and safety. We discuss the implementation of these scores in clinical medicine for uses including risk prediction and screening algorithms for coronary artery disease, prioritization of patient subgroups that are likely to derive benefit from treatment, and efficient prospective clinical trial designs.

From genotype to phenotype in Arabidopsis thaliana: in-silico genome interpretation predicts 288 phenotypes from sequencing data

In many cases, the unprecedented availability of data provided by high-throughput sequencing has shifted the bottleneck from a data availability issue to a data interpretation issue, thus delaying the promised breakthroughs in genetics and precision medicine, for what concerns Human genetics, and phenotype prediction to improve plant adaptation to climate change and resistance to bioagressors, for what concerns plant sciences. In this paper, we propose a novel Genome Interpretation paradigm, which aims at directly modeling the genotype-to-phenotype relationship, and we focus on A. thaliana since it is the best studied model organism in plant genetics. Our model, called Galiana, is the first end-to-end Neural Network (NN) approach following the genomes in/phenotypes out paradigm and it is trained to predict 288 real-valued Arabidopsis thaliana phenotypes from Whole Genome sequencing data. We show that 75 of these phenotypes are predicted with a Pearson correlation ≥0.4, and are mostly related to flowering traits. We show that our end-to-end NN approach achieves better performances and larger phenotype coverage than models predicting single phenotypes from the GWAS-derived known associated genes. Galiana is also fully interpretable, thanks to the Saliency Maps gradient-based approaches. We followed this interpretation approach to identify 36 novel genes that are likely to be associated with flowering traits, finding evidence for 6 of them in the existing literature.

Ultrafast homomorphic encryption models enable secure outsourcing of genotype imputation

Genotype imputation is a fundamental step in genomic data analysis, where missing variant genotypes are predicted using the existing genotypes of nearby "tag" variants. Although researchers can outsource genotype imputation, privacy concerns may prohibit genetic data sharing with an untrusted imputation service. Here, we developed secure genotype imputation using efficient homomorphic encryption (HE) techniques. In HE-based methods, the genotype data are secure while it is in transit, at rest, and in analysis. It can only be decrypted by the owner. We compared secure imputation with three state-of-the-art non-secure methods and found that HE-based methods provide genetic data security with comparable accuracy for common variants. HE-based methods have time and memory requirements that are comparable or lower than those for the non-secure methods. Our results provide evidence that HE-based methods can practically perform resource-intensive computations for high-throughput genetic data analysis. The source code is freely available for download at https://github.com/K-miran/secure-imputation.

Targeted re-sequencing and genome-wide association analysis for wood property traits in breeding population of Eucalyptus tereticornis × E. grandis

Globally, Eucalyptus plantations occupy 22 million ha area and is one of the preferred hardwood species due to their short rotation, rapid growth, adaptability and wood properties. In this study, we present results of GWAS in parents and 100 hybrids of Eucalyptus tereticornis × E. grandis using 762 genes presumably involved in wood formation. Comparative analysis between parents predicted 32,202 polymorphic SNPs with high average read depth of 269-562× per individual per nucleotide. Seventeen wood related traits were phenotyped across three diverse environments and GWAS was conducted using 13,610 SNPs. A total of 45 SNP-trait associations were predicted across two locations. Seven large effect markers were identified which explained more than 80% of phenotypic variation for fibre area. This study has provided an array of candidate genes which may govern fibre morphology in this genus and has predicted potential SNPs which can guide future breeding programs in tropical Eucalyptus.

Prediction of primary venous thromboembolism based on clinical and genetic factors within the U.K. Biobank

Both clinical and genetic factors drive the risk of venous thromboembolism. However, whether clinically recorded risk factors and genetic variants can be combined into a clinically applicable predictive score remains unknown. Using Cox proportional-hazard models, we analyzed the association of risk factors with the likelihood of venous thromboembolism in U.K. Biobank, a large prospective cohort. We then created a polygenic risk score of 36 single nucleotide polymorphisms and a clinical score determined by age, sex, body mass index, previous cancer diagnosis, smoking status, and fracture in the last 5 years. Participants were at significantly increased risk of venous thromboembolism if they were at high clinical risk (subhazard ratio, 4.37 [95% CI, 3.85-4.97]) or high genetic risk (subhazard ratio, 3.02 [95% CI, 2.63-3.47]) relative to participants at low clinical or genetic risk, respectively. The combined model, consisting of clinical and genetic components, was significantly better than either the clinical or the genetic model alone (P < 0.001). Participants at high risk in the combined score had nearly an eightfold increased risk of venous thromboembolism relative to participants at low risk (subhazard ratio, 7.51 [95% CI, 6.28-8.98]). This risk score can be used to guide decisions regarding venous thromboembolism prophylaxis, although external validation is needed.

A proximal LAVA method for genome-wide association and prediction of traits with mixed inheritance patterns

BACKGROUND

The genetic basis of phenotypic traits is highly variable and usually divided into mono-, oligo- and polygenic inheritance classes. Relatively few traits are known to be monogenic or oligogeneic. The majority of traits are considered to have a polygenic background. To what extent there are mixtures between these classes is unknown. The rapid advancement of genomic techniques makes it possible to directly map large amounts of genomic markers (GWAS) and predict unknown phenotypes (GWP). Most of the multi-marker methods for GWAS and GWP falls into one of two regularization frameworks. The first framework is based on [Formula: see text]-norm regularization (e.g. the LASSO) and is suitable for mono- and oligogenic traits, whereas the second framework regularize with the [Formula: see text]-norm (e.g. ridge regression; RR) and thereby is favourable for polygenic traits. A general framework for mixed inheritance is lacking.

RESULTS

We have developed a proximal operator algorithm based on the recent LAVA regularization method that jointly performs [Formula: see text]- and [Formula: see text]-norm regularization. The algorithm is built on the alternating direction method of multipliers and proximal translation mapping (LAVA ADMM). When evaluated on the simulated QTLMAS2010 data, it is shown that the LAVA ADMM together with Bayesian optimization of the regularization parameters provides an efficient approach with lower test prediction mean-squared-error (65.89) than the LASSO (66.11), Ridge regression (83.41) and Elastic net (66.11). For the real pig data the test MSE of the LAVA ADMM is 0.850 compared to the LASSO, RR and EN with 0.875, 0.853 and 0.853, respectively.

CONCLUSIONS

This study presents the LAVA ADMM that is capable of joint modelling of monogenic major genetic effects and polygenic minor genetic effects which can be used for both genome-wide assoiciation and prediction purposes. The statistical evaluations based on both simulated and real pig data set shows that the LAVA ADMM has better prediction properies than the LASSO, RR and EN. Julia code for the LAVA ADMM is available at: https://github.com/patwa67/LAVAADMM .

The population frequency of human mitochondrial DNA variants is highly dependent upon mutational bias

Next-generation sequencing can quickly reveal genetic variation potentially linked to heritable disease. As databases encompassing human variation continue to expand, rare variants have been of high interest, since the frequency of a variant is expected to be low if the genetic change leads to a loss of fitness or fecundity. However, the use of variant frequency when seeking genomic changes linked to disease remains very challenging. Here, I explored the role of selection in controlling human variant frequency using the HelixMT database, which encompasses hundreds of thousands of mitochondrial DNA (mtDNA) samples. I found that a substantial number of synonymous substitutions, which have no effect on protein sequence, were never encountered in this large study, while many other synonymous changes are found at very low frequencies. Further analyses of human and mammalian mtDNA datasets indicate that the population frequency of synonymous variants is predominantly determined by mutational biases rather than by strong selection acting upon nucleotide choice. My work has important implications that extend to the interpretation of variant frequency for non-synonymous substitutions.

No Casual Relationship Between T2DM and the Risk of Infectious Diseases: A Two-Sample Mendelian Randomization Study

Background

In epidemiological studies, it has been proven that the occurrence of type 2 diabetes mellitus (T2DM) is related to an increased risk of infectious diseases. However, it is still unclear whether the relationship is casual.

Methods

We employed a two-sample Mendelian randomization (MR) to clarify the causal effect of T2DM on high-frequency infectious diseases: sepsis, skin and soft tissue infections (SSTIs), urinary tract infections (UTIs), pneumonia, and genito-urinary infection (GUI) in pregnancy. And then, we analyzed the genome-wide association study (GWAS) meta-analysis of European-descent individuals and conducted T2DM-related single-nucleotide polymorphisms (SNPs) as instrumental variables (IVs) that were associated with genome-wide significance (p < 5 × 10^–8). MR estimates were obtained using the inverse variance-weighted (IVW), the MR-Egger regression, the simple mode (SM), weighted median, and weighted mode.

Results

The UK Biobank (UKB) cohort (n > 500,000) provided data for GWASs on infectious diseases. MR analysis showed little evidence of a causal relationship of T2DM with five mentioned infections’ (sepsis, SSTI, UTI, pneumonia, and GUI in pregnancy) susceptibility [odds ratio (OR) = 0.99999, p = 0.916; OR = 0.99986, p = 0.233; OR = 0.99973, p = 0.224; OR = 0.99997, p = 0.686; OR, 1.00002, p = 0.766]. Sensitivity analysis showed similar results, indicating the robustness of causality. There were no heterogeneity and pleiotropic bias.

Conclusion

T2DM would not be causally associated with high-frequency infectious diseases (including sepsis, SSTI, UTI, pneumonia, and GUI in pregnancy).

utr.annotation: a tool for annotating genomic variants that could influence post-transcriptional regulation

SUMMARY

Whole genome sequencing of patient populations is identifying thousands of new variants in UnTranslated Regions(UTRs). While the consequences of UTR mutations are not as easily predicted from primary sequence as coding mutations are, there are some known features of UTRs that modulate their function. utr.annotation is an R package that can be used to annotate potential deleterious variants in the UTR regions for both human and mouse species. Given a CSV or VCF format variant file, utr.annotation provides information of each variant on whether and how it alters known translational regulators including: upstream Open Reading Frames (uORFs), upstream Kozak sequences, polyA signals, Kozak sequences at the annotated translation start site, start codons, and stop codons, conservation scores in the variant position, and whether and how it changes ribosome loading based on a model derived from empirical data.

AVAILABILITY

utr.annotation is freely available on Bitbucket (https://bitbucket.org/jdlabteam/utr.annotation/src/master/) and CRAN (https://cran.r-project.org/web/packages/utr.annotation/index.html).

SUPPLEMENTARY INFORMATION

Supplementary data are available at https://wustl.box.com/s/yye99bryfin89nav45gv91l5k35fxo7z.

Germline CDH1 G212E Missense Variant: Combining Clinical, In Vitro and In Vivo Strategies to Unravel Disease Burden

Simple Summary

Hereditary diffuse gastric cancer (HDGC) is an inherited cancer syndrome associated with CDH1 germline mutations. The increasing detection of CDH1 genetic variants due to multigene panel testing poses a serious clinical challenge and urges the development of effective classification strategies. In this study, we describe the identification of the novel CDH1 G212E variant in a large family strongly affected by diffuse gastric cancer. Through a comprehensive characterization pipeline, we provide evidence of the damaging nature of this genetic alteration, thus impacting patient management and family screening.

Abstract

E-cadherin, encoded by CDH1, is an essential molecule for epithelial homeostasis, whose loss or aberrant expression results in disturbed cell–cell adhesion, increased cell invasion and metastasis. Carriers of CDH1 germline mutations have a high risk of developing diffuse gastric cancer and lobular breast cancer, associated with the cancer syndrome Hereditary Diffuse Gastric Cancer (HDGC). The ubiquitous availability of cancer panels has led to the identification of an increasing amount of “incidental” CDH1 genetic variants that pose a serious clinical challenge. This has sparked intensive research aiming at an accurate classification of the variants and consequent validation of their clinical relevance. The present study addressed the significance of a novel CDH1 variant, G212E, identified in an unusually large pedigree displaying strong aggregation of diffuse gastric cancer. We undertook a comprehensive pipeline encompassing family data, in silico predictions, in vitro assays and in vivo strategies, which validated the deleterious phenotype induced by this genetic alteration. In particular, we demonstrated that the G212E variant affects the stability and localization, as well as the adhesive and anti-invasive functions of E-cadherin, triggering epithelial disruption and disorganization. Our findings illustrate the clinical implication of a complementary approach for effective variant categorization and patient management.

Genome-wide association studies of cardiac electrical phenotypes

The genetic basis of cardiac electrical phenotypes has in the last 25 years been the subject of intense investigation. While in the first years, such efforts were dominated by the study of familial arrhythmia syndromes, in recent years, large consortia of investigators have successfully pursued genome-wide association studies (GWAS) for the identification of single-nucleotide polymorphisms that govern inter-individual variability in electrocardiographic parameters in the general population. We here provide a review of GWAS conducted on cardiac electrical phenotypes in the last 14 years and discuss the implications of these discoveries for our understanding of the genetic basis of disease susceptibility and variability in disease severity. Furthermore, we review functional follow-up studies that have been conducted on GWAS loci associated with cardiac electrical phenotypes and highlight the challenges and opportunities offered by such studies.

Natural human knockouts and Mendelian disorders: deep phenotyping in Italian isolates

Whole genome sequencing (WGS) allows the identification of human knockouts (HKOs), individuals in whom loss of function (LoF) variants disrupt both alleles of a given gene. HKOs are a valuable model for understanding the consequences of genes function loss. Naturally occurring biallelic LoF variants tend to be significantly enriched in "genetic isolates," making these populations specifically suited for HKO studies. In this work, a meticulous WGS data analysis combined with an in-depth phenotypic assessment of 947 individuals from three Italian genetic isolates led to the identification of ten biallelic LoF variants in ten OMIM genes associated with known autosomal recessive diseases. Notably, only a minority of the identified HKOs (C7, F12, and GPR68 genes) displayed the expected phenotype. For most of the genes, instead, (ACADSB, FANCL, GRK1, LGI4, MPO, PGAM2, and RP1L1), the carriers showed none or few of the signs and symptoms typically associated with the related diseases. Of particular interest is a case presenting with a FANCL biallelic LoF variant and a positive diepoxybutane test but lacking a full Fanconi anemia phenotypic spectrum. Identifying KO subjects displaying expected phenotypes suggests that the lack of correct genetic diagnoses may lead to inappropriate and delayed treatment. In contrast, the presence of HKOs with phenotypes deviating from the expected patterns underlines how LoF variants may be responsible for broader phenotypic spectra. Overall, these results highlight the importance of in-depth phenotypical characterization to understand the role of LoF variants and the advantage of studying these variants in genetic isolates.

Environmentally induced ribosomal DNA (rDNA) instability in human cells and populations exposed to hexavalent chromium [Cr (VI)]

Hexavalent Chromium [Cr (VI)] is an established toxicant, carcinogen, and a significant source of public health concern. The multicopy ribosomal DNA (rDNA) array is mechanistically linked to aging and cancer, is the most evolutionarily conserved segment of the human genome, and gives origin to nucleolus, a nuclear organelle where ribosomes are assembled. Here we show that exposure to Cr (VI) induces instability in the rDNA, triggering cycles of rapid, specific, and transient amplification and contraction of the array in human cells. The dynamic of environmentally responsive rDNA copy number (CN) amplification and contraction occurs at doses to which millions of individuals are regularly exposed. Finally, analyses of human populations occupationally exposed to Cr (VI) indicate that environmental exposure history and drinking habits but not age shape extensive naturally occurring rDNA copy number variation. Our observations identify a novel pathway of response to hexavalent chromium exposure and raise the prospect that a suite of environmental determinants of rDNA copy number remain to be discovered.

Polygenic Risk Score of Adolescent Idiopathic Scoliosis for Potential Clinical Use

Adolescent idiopathic scoliosis (AIS) is a common disease causing three-dimensional spinal deformity in as many as 3% of adolescents. Development of a method that can accurately predict the onset and progression of AIS is an immediate need for clinical practice. Because the heritability of AIS is estimated as high as 87.5% in twin studies, prediction of its onset and progression based on genetic data is a promising option. We show the usefulness of polygenic risk score (PRS) for the prediction of onset and progression of AIS. We used AIS genomewide association study (GWAS) data comprising 79,211 subjects in three cohorts and constructed a PRS based on association statistics in a discovery set including 31,999 female subjects. After calibration using a validation data set, we applied the PRS to a test data set. By integrating functional annotations showing heritability enrichment in the selection of variants, the PRS demonstrated an association with AIS susceptibility (p = 3.5 × 10^-40 with area under the receiver-operating characteristic [AUROC] = 0.674, sensitivity = 0.644, and specificity = 0.622). The decile with the highest PRS showed an odds ratio of as high as 3.36 (p = 1.4 × 10^-10 ) to develop AIS compared with the fifth in decile. The addition of a predictive model with only a single clinical parameter (body mass index) improved predictive ability for development of AIS (AUROC = 0.722, net reclassification improvement [NRI] 0.505 ± 0.054, p = 1.6 × 10^-8 ), potentiating clinical use of the prediction model. Furthermore, we found the Cobb angle (CA), the severity measurement of AIS, to be a polygenic trait that showed a significant genetic correlation with AIS susceptibility (rg = 0.6, p = 3.0 × 10^-4 ). The AIS PRS demonstrated a significant association with CA. These results indicate a shared polygenic architecture between onset and progression of AIS and the potential usefulness of PRS in clinical settings as a predictor to promote early intervention of AIS and avoid invasive surgery. © 2021 American Society for Bone and Mineral Research (ASBMR).

Evidence of gene-gene interaction in Hidradenitis suppurativa - A nationwide register study of Danish twins

BACKGROUND

Hidradenitis suppurativa (HS) is a recurrent inflammatory skin disease that, apart from rare causative loss-of-function mutations, has a widely unknown genetic aetiology. Our objective was to estimate the relative importance of genetic and environmental factors underlying HS susceptibility.

METHODS

Through the Danish Twin Registry and the Danish National Patient Registry we joined information on zygosity with that of HS status. HS cases were identified by International Code of Diseases 8 (705.91) and 10 (L73.2). Heritability was assessed by the classic biometric model and the possibility of gene-gene interaction through the multi-locus modeling approach.

RESULTS

Amongst 100,044 registered twins, we found 170 twins (from 163 pairs) diagnosed with HS. The seven concordant pairs were all monozygotic, and monozygotic twins had a casewise concordance rate of 28% (95% CI: 7%; 49%), corresponding to a familial risk of 73 (95% CI 13; 133) times that of the background population. The biometrical modelling suggested a heritability of 0.80 (95% CI 0.67; 0.93), and the multilocus index estimate was 230 (95% CI: 60; 400). This is highly indicative of gene-gene interactions, with the possibility of up to six interacting loci.

CONCLUSION

This twin study is substantially larger, and employs a more valid phenotype than prior studies. Genetics account for the majority of the HS susceptibility, and HS is most likely caused by gene-gene interactions rather than monogenetic mutations or solely additive genetic factors. New approaches aimed at assessing potential interactions at a SNP-SNP level should be implemented in future HS genome-wide association studies.

Contextualizing genetic risk score for disease screening and rare variant discovery

Studies of the genetic basis of complex traits have demonstrated a substantial role for common, small-effect variant polygenic burden (PB) as well as large-effect variants (LEV, primarily rare). We identify sufficient conditions in which GWAS-derived PB may be used for well-powered rare pathogenic variant discovery or as a sample prioritization tool for whole-genome or exome sequencing. Through extensive simulations of genetic architectures and generative models of disease liability with parameters informed by empirical data, we quantify the power to detect, among cases, a lower PB in LEV carriers than in non-carriers. Furthermore, we uncover clinically useful conditions wherein the risk derived from the PB is comparable to the LEV-derived risk. The resulting summary-statistics-based methodology (with publicly available software, PB-LEV-SCAN) makes predictions on PB-based LEV screening for 36 complex traits, which we confirm in several disease datasets with available LEV information in the UK Biobank, with important implications on clinical decision-making.

From QTL to gene: C. elegans facilitates discoveries of the genetic mechanisms underlying natural variation

Although many studies have examined quantitative trait variation across many species, only a small number of genes and thereby molecular mechanisms have been discovered. Without these data, we can only speculate about evolutionary processes that underlie trait variation. Here, we review how quantitative and molecular genetics in the nematode Caenorhabditis elegans led to the discovery and validation of 37 quantitative trait genes over the past 15 years. Using these data, we can start to make inferences about evolution from these quantitative trait genes, including the roles that coding versus noncoding variation, gene family expansion, common versus rare variants, pleiotropy, and epistasis play in trait variation across this species.

Recalling the pathology of Parkinson's disease; lacking exact figure of prevalence and genetic evidence in Asia with an alarming outcome: A time to step-up

Parkinson's disease (PD) is the second most common and progressive neurodegenerative disease globally, with major symptoms like bradykinesia, impaired posture, and tremor. Several genetic and environmental factors have been identified but elucidating the main factors have been challenging due to the disease's complex nature. Diagnosis, prognosis, and management of such diseases are challenging and require effective targeted attention in developing countries. Recently, PD is growing rapidly in many crowded Asian countries as an alarming threat with inadequate knowledge of its prevalence, genetic architecture, and geographic distribution. This study gave an in-depth overview of the prevalence, incidence and genomic/genetics studies published so far in the Asian population. To the best of our knowledge, PD has increased significantly in several Asian countries, including China, South Korea, Japan, Thailand, and Israel over the past few years, requiring a greater level of care and attention. Genetic screening of families with PD at national levels and establishing an official database of PD cases are essential to get a comprehensive and conclusive view of the exact prevalence and genetic diversity of PD in the Asian population to properly manage and treat the disease.

Candidate Rlm6 resistance genes against Leptosphaeria. maculans identified through a genome-wide association study in Brassica juncea (L.) Czern

One hundred and sixty-seven B. juncea varieties were genotyped on the 90K Brassica assay (42,914 SNPs), which led to the identification of sixteen candidate genes for Rlm6. Brassica species are at high risk of severe crop loss due to pathogens, especially Leptosphaeria maculans (the causal agent of blackleg). Brassica juncea (L.) Czern is an important germplasm resource for canola improvement, due to its good agronomic traits, such as heat and drought tolerance and high blackleg resistance. The present study is the first using genome-wide association studies to identify candidate genes for blackleg resistance in B. juncea based on genome-wide SNPs obtained from the Illumina Infinium 90 K Brassica SNP array. The verification of Rlm6 in B. juncea was performed through a cotyledon infection test. Genotyping 42,914 single nucleotide polymorphisms (SNPs) in a panel of 167 B. juncea lines revealed a total of seven SNPs significantly associated with Rlm6 on chromosomes A07 and B04 in B. juncea. Furthermore, 16 candidate Rlm6 genes were found in these regions, defined as nucleotide binding site leucine-rich-repeat (NLR), leucine-rich repeat RLK (LRR-RLK) and LRR-RLP genes. This study will give insights into the blackleg resistance in B. juncea and facilitate identification of functional blackleg resistance genes which can be used in Brassica breeding.

Genome-wide association studies: assessing trait characteristics in model and crop plants

GWAS involves testing genetic variants across the genomes of many individuals of a population to identify genotype–phenotype association. It was initially developed and has proven highly successful in human disease genetics. In plants genome-wide association studies (GWAS) initially focused on single feature polymorphism and recombination and linkage disequilibrium but has now been embraced by a plethora of different disciplines with several thousand studies being published in model and crop species within the last decade or so. Here we will provide a comprehensive review of these studies providing cases studies on biotic resistance, abiotic tolerance, yield associated traits, and metabolic composition. We also detail current strategies of candidate gene validation as well as the functional study of haplotypes. Furthermore, we provide a critical evaluation of the GWAS strategy and its alternatives as well as future perspectives that are emerging with the emergence of pan-genomic datasets.

SINE jumping contributes to large-scale polymorphisms in the pig genomes

BACKGROUND

Molecular markers based on retrotransposon insertion polymorphisms (RIPs) have been developed and are widely used in plants and animals. Short interspersed nuclear elements (SINEs) exert wide impacts on gene activity and even on phenotypes. However, SINE RIP profiles in livestock remain largely unknown, and not be revealed in pigs.

RESULTS

Our data revealed that SINEA1 displayed the most polymorphic insertions (22.5 % intragenic and 26.5 % intergenic), followed by SINEA2 (10.5 % intragenic and 9 % intergenic) and SINEA3 (12.5 % intragenic and 5.0 % intergenic). We developed a genome-wide SINE RIP mining protocol and obtained a large number of SINE RIPs (36,284), with over 80 % accuracy and an even distribution in chromosomes (14.5/Mb), and 74.34 % of SINE RIPs generated by SINEA1 element. Over 65 % of pig SINE RIPs overlap with genes, most of them (> 95 %) are in introns. Overall, about one forth (23.09 %) of the total genes contain SINE RIPs. Significant biases of SINE RIPs in the transcripts of protein coding genes were observed. Nearly half of the RIPs are common in these pig breeds. Sixteen SINE RIPs were applied for population genetic analysis in 23 pig breeds, the phylogeny tree and cluster analysis were generally consistent with the geographical distributions of native pig breeds in China.

CONCLUSIONS

Our analysis revealed that SINEA1-3 elements, particularly SINEA1, are high polymorphic across different pig breeds, and generate large-scale structural variations in the pig genomes. And over 35,000 SINE RIP markers were obtained. These data indicate that young SINE elements play important roles in creating new genetic variations and shaping the evolution of pig genome, and also provide strong evidences to support the great potential of SINE RIPs as genetic markers, which can be used for population genetic analysis and quantitative trait locus (QTL) mapping in pig.

Advances in colorectal cancer genomics and transcriptomics drive early detection and prevention

Colorectal carcinoma (CRC) is a high incidence cancer and leading cause of cancer mortality worldwide. The advances in genomics and transcriptomics in the past decades have improved the detection and prevention of CRC in familial CRC syndromes. Nevertheless, the ultimate goal of personalized medicine for sporadic CRC is still not within reach due no less to the difficulty in integrating population disparity and clinical data to combat what essentially is a very heterogenous disease. This minireview highlights the achievement of the past decades and present possible direction in the hope of early detection and metastasis prevention for reducing CRC-associated morbidity and mortality.

Sequencing at lymphoid neoplasm susceptibility loci maps six myeloma risk genes

Inherited genetic risk factors play a role in multiple myeloma (MM), yet considerable missing heritability exists. Rare risk variants at genome-wide association study (GWAS) loci are a new avenue to explore. Pleiotropy between lymphoid neoplasms (LNs) has been suggested in family history and genetic studies, but no studies have interrogated sequencing for pleiotropic genes or rare risk variants. Sequencing genetically enriched cases can help discover rarer variants. We analyzed exome sequencing in familial or early-onset MM cases to identify rare, functionally relevant variants near GWAS loci for a range of LNs. A total of 149 distinct and significant LN GWAS loci have been published. We identified six recurrent, rare, potentially deleterious variants within 5 kb of significant GWAS single nucleotide polymorphisms in 75 MM cases. Mutations were observed in BTNL2, EOMES, TNFRSF13B, IRF8, ACOXL and TSPAN32. All six genes replicated in an independent set of 255 early-onset MM or familial MM or precursor cases. Expansion of our analyses to the full length of these six genes resulted in a list of 39 rare and deleterious variants, seven of which segregated in MM families. Three genes also had significant rare variant burden in 733 sporadic MM cases compared with 935 control individuals: IRF8 (P = 1.0 × 10-6), EOMES (P = 6.0 × 10-6) and BTNL2 (P = 2.1 × 10-3). Together, our results implicate six genes in MM risk, provide support for genetic pleiotropy between LN subtypes and demonstrate the utility of sequencing genetically enriched cases to identify functionally relevant variants near GWAS loci.

Rare variants regulate expression of nearby individual genes in multiple tissues

The rapid decrease in sequencing cost has enabled genetic studies to discover rare variants associated with complex diseases and traits. Once this association is identified, the next step is to understand the genetic mechanism of rare variants on how the variants influence diseases. Similar to the hypothesis of common variants, rare variants may affect diseases by regulating gene expression, and recently, several studies have identified the effects of rare variants on gene expression using heritability and expression outlier analyses. However, identifying individual genes whose expression is regulated by rare variants has been challenging due to the relatively small sample size of expression quantitative trait loci studies and statistical approaches not optimized to detect the effects of rare variants. In this study, we analyze whole-genome sequencing and RNA-seq data of 681 European individuals collected for the Genotype-Tissue Expression (GTEx) project (v8) to identify individual genes in 49 human tissues whose expression is regulated by rare variants. To improve statistical power, we develop an approach based on a likelihood ratio test that combines effects of multiple rare variants in a nonlinear manner and has higher power than previous approaches. Using GTEx data, we identify many genes regulated by rare variants, and some of them are only regulated by rare variants and not by common variants. We also find that genes regulated by rare variants are enriched for expression outliers and disease-causing genes. These results suggest the regulatory effects of rare variants, which would be important in interpreting associations of rare variants with complex traits.

It has been shown that rare variants may affect many diseases including both rare and common diseases with the advent of next-generation sequencing technology. An important question is how rare variants affect diseases or traits, especially whether or how they regulate gene expression as they may affect diseases through gene regulation. However, it is challenging to identify the regulatory effects of rare variants because it often requires large sample sizes and the existing statistical approaches are not optimized for it. Here, we develop a novel method, LRT-q, based on a likelihood ratio test that aggregates the effects of multiple rare variants nonlinearly to achieve higher statistical power than previous rare variant association methods. We apply LRT-q to the latest GTEx v8 dataset and identify regulatory effect of rare variants on individual genes. We also observe that genes regulated by rare variants are likely to be disease-causing genes. These results demonstrate the functional effects of rare variants, especially on gene expression, which provides important biological insights in understanding the genetic mechanism of rare variants in complex traits and diseases.

Childhood risk factors for adulthood chronic kidney disease

Chronic kidney disease (CKD) is a major public health challenge, affecting as much as 8 to 18% of the world population. Identifying childhood risk factors for future CKD may help clinicians make early diagnoses and initiation of preventive interventions for CKD and its attendant comorbidities as well as monitoring for complications. The purpose of this review is to describe childhood risk factors that may predict development of overt kidney disease later in life. Currently, there are multiple childhood risk factors associated with future onset and progression of CKD. These risk factors can be grouped into five categories: genetic factors (e.g., monogenic or risk alleles), perinatal factors (e.g., low birth weight and prematurity), childhood kidney diseases (e.g., congenital anomalies, glomerular diseases, and renal cystic ciliopathies), childhood onset of chronic conditions (e.g., cancer, diabetes, hypertension, dyslipidemia, and obesity), and different lifestyle factors (e.g., physical activity, diet, and factors related to socioeconomic status). The available published information suggests that the lifelong risk for CKD can be attributed to multiple factors that appear already during childhood. However, results are conflicting on the effects of childhood physical activity, diet, and dyslipidemia on future renal function. On the other hand, there is consistent evidence to support follow-up of high-risk groups.

Advances and challenges in quantitative delineation of the genetic architecture of complex traits

Background

Genome-wide association studies (GWAS) have been widely adopted in studies of human complex traits and diseases.

Results

This review surveys areas of active research: quantifying and partitioning trait heritability, fine mapping functional variants and integrative analysis, genetic risk prediction of phenotypes, and the analysis of sequencing studies that have identified millions of rare variants. Current challenges and opportunities are highlighted.

Conclusion

GWAS have fundamentally transformed the field of human complex trait genetics. Novel statistical and computational methods have expanded the scope of GWAS and have provided valuable insights on the genetic architecture underlying complex phenotypes.

Statistical methods for mediation analysis in the era of high-throughput genomics: Current successes and future challenges

Mediation analysis investigates the intermediate mechanism through which an exposure exerts its influence on the outcome of interest. Mediation analysis is becoming increasingly popular in high-throughput genomics studies where a common goal is to identify molecular-level traits, such as gene expression or methylation, which actively mediate the genetic or environmental effects on the outcome. Mediation analysis in genomics studies is particularly challenging, however, thanks to the large number of potential mediators measured in these studies as well as the composite null nature of the mediation effect hypothesis. Indeed, while the standard univariate and multivariate mediation methods have been well-established for analyzing one or multiple mediators, they are not well-suited for genomics studies with a large number of mediators and often yield conservative p-values and limited power. Consequently, over the past few years many new high-dimensional mediation methods have been developed for analyzing the large number of potential mediators collected in high-throughput genomics studies. In this work, we present a thorough review of these important recent methodological advances in high-dimensional mediation analysis. Specifically, we describe in detail more than ten high-dimensional mediation methods, focusing on their motivations, basic modeling ideas, specific modeling assumptions, practical successes, methodological limitations, as well as future directions. We hope our review will serve as a useful guidance for statisticians and computational biologists who develop methods of high-dimensional mediation analysis as well as for analysts who apply mediation methods to high-throughput genomics studies.

A genetic polymorphism that is associated with mitochondrial energy metabolism increases risk of fibromyalgia

Alterations in cellular energy metabolism have been implicated in chronic pain, suggesting a role for mitochondrial DNA. Previous studies reported associations of a limited number of mitochondrial DNA polymorphisms with specific pain conditions. In this study, we examined the full mitochondrial genomes of people with a variety of chronic pain conditions. A discovery cohort consisting of 609 participants either with or without a complex persistent pain conditions (CPPCs) was examined. Mitochondrial DNA was subjected to deep sequencing for identification of rare mutations, common variants, haplogroups, and heteroplasmy associated with 5 CPPCs: episodic migraine, irritable bowel syndrome, fibromyalgia, vulvar vestibulitis, or temporomandibular disorders. The strongest association found was the presence of the C allele at the single nucleotide polymorphism m.2352T>C (rs28358579) that significantly increased the risk for fibromyalgia (odds ratio [OR] = 4.6, P = 4.3 × 10). This relationship was even stronger in women (OR = 5.1, P = 2.8 × 10), and m.2352T>C was associated with all other CPPCs in a consistent risk-increasing fashion. This finding was replicated in another cohort (OR = 4.3, P = 2.6 × 10) of the Orofacial Pain: Prospective Evaluation and Risk Assessment study consisting of 1754 female participants. To gain insight into the cellular consequences of the associated genetic variability, we conducted an assay testing metabolic reprogramming in human cell lines with defined genotypes. The minor allele C was associated with decreased mitochondrial membrane potential under conditions where oxidative phosphorylation is required, indicating a role of oxidative phosphorylation in pathophysiology of chronic pain. Our results suggest that cellular energy metabolism, modulated by m.2352T>C, contributes to fibromyalgia and possibly other chronic pain conditions.

Germline and somatic mutations in the pathology of pineal cyst: A whole-exome sequencing study of 93 individuals

Background

Pineal cyst is a benign lesion commonly occurring in people of any age. Until now, the underlying molecular alterations have not been explored.

Methods

We performed whole exome sequencing of 93 germline samples and 21 pineal cyst tissue samples to illustrate its genetic architecture and somatic mutations. The dominant and recessive inheritance modes were considered, and a probability was calculated to evaluate the significance of variant overrepresentation.

Results

By analyzing pineal cyst as a Mendelian disease with a dominant inheritance pattern, we identified 42,325 rare germline variants, and NM_001004711.1:c.476A>G was highly enriched (FDR<0.2). By analyzing it as a recessive disorder, we identified 753 homozygous rare variants detected in at least one pineal cyst sample each. One STIM2 rare variant, NM_001169117.1:c.1652C>T, was overrepresented (FDR<0.05). Analyzing at a gene‐based level, we identified a list of the most commonlymutated germline genes, including POP4, GNGT2 and TMEM254. A somatic mutation analysis of 21 samples identified 16 variants in 15 genes, which mainly participated in the biological processes of gene expression and epigenetic regulation, immune response modulation, and transferase activity.

Conclusion

These molecular profiles are novel for this condition and provide data for investigators interested in pineal cysts.

Genetic variability in COVID-19-related genes in the Brazilian population

SARS-CoV-2 utilizes the angiotensin-converting enzyme 2 (ACE2) receptor and transmembrane serine protease (TMPRSS2) to infect human lung cells. Previous studies have suggested that different host ACE2 and TMPRSS2 genetic backgrounds might contribute to differences in the rate of SARS-CoV-2 infection or COVID-19 severity. Recent studies have also shown that variants in 15 genes related to type I interferon immunity to influenza virus might predispose patients toward life-threatening COVID-19 pneumonia. Other genes (SLC6A20, LZTFL1, CCR9, FYCO1, CXCR6, XCR1, IL6, CTSL, ABO, and FURIN) and HLA alleles have also been implicated in the response to infection with SARS-CoV-2. Currently, Brazil has recorded the third-highest number of COVID-19 cases worldwide. We aimed to investigate the genetic variation present in COVID-19-related genes in the Brazilian population. We analyzed 27 candidate genes and HLA alleles in 954 admixed Brazilian exomes. We used the information available in two public databases (http://www.bipmed.org and http://abraom.ib.usp.br/) and additional exomes from individuals born in southeast Brazil, the region of the country with the highest number of COVID-19 patients. Variant allele frequencies were compared with the 1000 Genomes Project phase 3 (1KGP) and gnomAD databases. We detected 395 nonsynonymous variants; of these, 325 were also found in the 1KGP and/or gnomAD. Six of these variants were previously reported to influence the rate of infection or clinical prognosis of COVID-19. The remaining 70 variants were identified exclusively in the Brazilian sample, with a mean allele frequency of 0.0025. In silico analysis revealed that seven of these variants are predicted to affect protein function. Furthermore, we identified HLA alleles previously associated with the COVID-19 response at loci DQB1 and DRB1. Our results showed genetic variability common to other populations and rare and ultrarare variants exclusively found in the Brazilian population. These findings might lead to differences in the rate of infection or response to infection by SARS-CoV-2 and should be further investigated in patients with this disease.

Genetic variants in the Brazilian population do not appear to explain the relatively high overall rates of COVID-19 cases compared to other countries, but could be relevant for risk assessment at the individual level. Iscia Lopes-Cendes of the University of Campinas and colleagues in Brazil examined gene sequences for variations in 27 genes reported to influence COVID-19 infection in Brazilians of mixed heritage and of people born in the southeast where case numbers are high. These studies included the ACE2 gene that codes for a receptor allowing SARS-CoV-2 to enter host cells. Genetic variants were also identified in two international genome databases. Seventy variants were unique to Brazilians, but computer modeling predicted only seven variants that could affect protein function. The team also looked for and found variants that could affect individual immunity to COVID-19. Similar genetic studies could help identify at-risk individuals and therapeutic targets.

Genetic analysis of 39 erythrocytosis and hereditary hemochromatosis-associated genes in the Slovenian family with idiopathic erythrocytosis

BACKGROUND

Erythrocytosis is a condition with an excessive number of erythrocytes, accompanied by an elevated haemoglobin and/or haematocrit value. Congenital erythrocytosis has a diverse genetic background with several genes involved in erythropoiesis. In clinical practice, nine genes are usually examined, but in approximately 70% of patients, no causative mutation can be identified. In this study, we screened 39 genes, aiming to identify potential disease-driving variants in the family with erythrocytosis of unknown cause.

PATIENTS AND METHODS

Two affected family members with elevated haemoglobin and/or haematocrit and negative for acquired causes and one healthy relative from the same family were selected for molecular-genetic analysis of 24 erythrocytosis and 15 hereditary haemochromatosis-associated genes with targeted NGS. The identified variants were further analysed for pathogenicity using various bioinformatic tools and review of the literature.

RESULTS

Of the 12 identified variants, two heterozygous variants, the missense variant c.471G>C (NM_022051.2) (p.(Gln157His)) in the EGLN1 gene and the intron variant c.2572-13A>G (NM_004972.3) in the JAK2 gene, were classified as low-frequency variants in European population. None of the two variants were present in a healthy family member. Variant c.2572-13A>G has potential impact on splicing by one prediction tool.

CONCLUSION

For the first time, we included 39 genes in the erythrocytosis clinical panel and identified two potential disease-driving variants in the Slovene family studied. Based on the reported functional in vitro studies combined with our bioinformatics analysis, we suggest further functional analysis of variant in the JAK2 gene and evaluation of a cumulative effect of both variants.