Despite the growing number of genome-wide association studies (GWASs), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWASs are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate, although this is not significant, most likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits.

Following muscle injury, fibro-adipogenic progenitors (FAPs) are rapidly activated and undergo apoptosis at the resolution stage, which is required for proper muscle regeneration. When excessive FAPs remain, it contributes to fibrotic and fatty infiltration, impairing muscle recovery. Mechanisms controlling FAP apoptosis remain poorly defined. We hypothesized that AMP-activated protein kinase (AMPK) in FAPs mediates their apoptosis during the muscle regeneration.

Hepatocellular carcinoma (HCC) is one of the most common primary hepatic malignancies. E2F transcription factors play an important role in the tumorigenesis and progression of HCC, mainly through the RB/E2F pathway. Prognostic models for HCC based on gene signatures have been developed rapidly in recent years; however, their discriminating ability at the single-cell level remains elusive, which could reflect the underlying mechanisms driving the sample bifurcation. In this study, we constructed and validated a predictive model based on E2F expression, successfully stratifying patients with HCC into two groups with different survival risks. Then we used a single-cell dataset to test the discriminating ability of the predictive model on infiltrating T cells, demonstrating remarkable cellular heterogeneity as well as altered cell fates. We identified distinct cell subpopulations with diverse molecular characteristics. We also found that the distribution of cell subpopulations varied considerably across onset stages among patients, providing a fundamental basis for patient-oriented precision evaluation. Moreover, single-sample gene set enrichment analysis revealed that subsets of CD8+ T cells with significantly different cell adhesion levels could be associated with different patterns of tumor cell dissemination. Therefore, our findings linked the conventional prognostic gene signature to the immune microenvironment and cellular heterogeneity at the single-cell level, thus providing deeper insights into the understanding of HCC tumorigenesis.

Hypertrophic cardiomyopathy (HCM) represents one of the most common inherited heart diseases. To identify key molecules involved in the development of HCM, gene expression patterns of the heart tissue samples in HCM patients from multiple microarray and RNA-seq platforms were investigated.

Gastric cancer (GC) remains an important cancer worldwide. Further understanding of the molecular mechanisms of gastric carcinogenesis will enhance the diagnosis and treatment of GC.

Impaired lung function is often caused by cigarette smoking, making it challenging to disentangle its role in lung cancer susceptibility. Investigation of the shared genetic basis of these phenotypes in the UK Biobank and International Lung Cancer Consortium (29,266 cases, 56,450 controls) shows that lung cancer is genetically correlated with reduced forced expiratory volume in one second (FEV1: rg = 0.098, p = 2.3 × 10-8) and the ratio of FEV1 to forced vital capacity (FEV1/FVC: rg = 0.137, p = 2.0 × 10-12). Mendelian randomization analyses demonstrate that reduced FEV1 increases squamous cell carcinoma risk (odds ratio (OR) = 1.51, 95% confidence intervals: 1.21-1.88), while reduced FEV1/FVC increases the risk of adenocarcinoma (OR = 1.17, 1.01-1.35) and lung cancer in never smokers (OR = 1.56, 1.05-2.30). These findings support a causal role of pulmonary impairment in lung cancer etiology. Integrative analyses reveal that pulmonary function instruments, including 73 novel variants, influence lung tissue gene expression and implicate immune-related pathways in mediating the observed effects on lung carcinogenesis.

Acquired uniparental disomy (aUPD) leads to homozygosity facilitating identification of monoallelically expressed genes. We analyzed single-nucleotide polymorphism array-based genotyping data of 448 head and neck squamous cell carcinoma (HNSCC) samples from The Cancer Genome Atlas to determine the frequency and distribution of aUPD regions and their association with survival, as well as to gain a better understanding of their influence on the tumor genome. We used expression data from the same dataset to identify differentially expressed genes between groups with and without aUPD. Univariate and multivariable Cox proportional hazards models were performed for survival analysis. We found that 82.14% of HNSCC samples carried aUPD; the most common regions were in chromosome 17p (31.25%), 9p (30.13%), and 9q (27.46%). In univariate analysis, five independent aUPD regions at chromosome 9p, two regions at chromosome 9q, and the CDKN2A region were associated with poor overall survival in all groups, including training and test sets and human papillomavirus (HPV)-negative samples. Forty-three genes in areas of aUPD including PD-L1 and CDKN2A were differentially expressed in samples with aUPD compared to samples without aUPD. In multivariable analysis, aUPD at the CDKN2A region was a significant predictor of overall survival in the whole cohort and in patients with HPV-negative HNSCC. aUPD at specific regions in the genome influences clinical outcomes of HNSCC and may be beneficial for selection of personalized therapy to prolong survival in patients with this disease.

Danon disease (DD) is an X-linked dominant multisystem disorder that is associated with cardiomyopathy, skeletal myopathy, and varying degrees of intellectual disability. It results from mutations in the lysosome-associated membrane protein 2 (LAMP2) gene.

Neo-tetraploid rice (NTR) is a new tetraploid rice germplasm that developed from the crossing and directional selection of different autotetraploid rice lines, which showed high fertility and promising yield potential. However, systematic yield assessment, genome composition and functional variations associated with fertility and yield remain elusive.

Neo-tetraploid rice (NTR) is a useful new germplasm that developed from the descendants of the autotetraploid rice (ATR) hybrids. NTR showed improved fertility and yield potential, and produced high yield heterosis when crossed with indica ATR for commercial utilization. However, their classification, population structure and genomic feature remain elusive.

Although genome-wide association studies have been conducted to investigate genetic variation of lung tumorigenesis, little is known about gene-gene (G × G) interactions that may influence the risk of non-small cell lung cancer (NSCLC).

Autotetraploid rice exhibited hybrid vigor and greater genetic variation compared to diploid rice, but low pollen fertility is a major hindrance for its utilization. Our previous analysis revealed that large number of pollen fertility genes were exhibited down-regulation in autotetraploid rice. Hence, it is of utmost importance to reveal the expression patterns of pollen fertility genes with high accuracy. To find stable reference genes for autotetraploid rice, we compared the pollen development stages between diploid and autotetraploid rice, and 14 candidate genes were selected based on transcriptome analysis to evaluate their expression levels. Autotetraploid rice (i.e. Taichung65-4x) displayed lower seed set (40.40%) and higher percentage of abnormalities during the pollen development process than its diploid counterpart. To detect the candidate reference genes for pollen development of autotetraploid and diploid rice, we used five different algorithms, including NormFinder, BestKeeper, ΔCt method, geNorm and Re-Finder to evaluate their expression patterns stability. Consequently, we identified two genes, Cytochrome b5 and CPI, as the best candidate reference genes for qRT-PCR normalization in autotetraploid and diploid rice during pre-meiosis, meiosis, single microspore and bicellular pollen development stages. However, Cytochrome b5 was found to be the most stably expressed gene during different pollen development stages in autotetraploid rice. The results of our study provide a platform for subsequent gene expression analyses in autotetraploid rice, which could also be used in other polyploid plants.

This study examined individual-level determinants of self-reported changes in healthy (diet and physical activity) and addictive (alcohol use, smoking, and vaping) lifestyle behaviors during the initial COVID-19 lockdown period in the USA. A national online survey was administered between May and June 2020 that targeted a representative U.S. sample and yielded data from 1276 respondents, including 58% male and 50% racial/ethnic minorities. We used univariate and multivariable linear regression models to examine the associations of sociodemographic, mental health, and behavioral determinants with self-reported changes in lifestyle behaviors. Some study participants reported increases in healthy lifestyle behaviors since the pandemic (i.e., 36% increased healthy eating behaviors, and 33% increased physical activity). However, they also reported increases in addictive lifestyle behaviors including alcohol use (40%), tobacco use (41%), and vaping (46%). With regard to individual-level determinants, individuals who reported adhering to social distancing guidelines were also more likely to report increases in healthy lifestyle behaviors (β = 0.12, 95% CI 0.04 to 0.21). Conversely, women (β = -0.37, 95% CI -0.62 to -0.12), and unemployed individuals (β = -0.33, 95% CI -0.64 to -0.02) were less likely to report increases in healthy lifestyle behaviors. In addition, individuals reporting anxiety were more likely to report increases in addictive behaviors (β = 0.26, 95% CI 0.09 to 0.43). Taken together, these findings suggest that women and unemployed individuals may benefit from interventions targeting diet and physical activity, and that individuals reporting anxiety may benefit from interventions targeting smoking and alcohol cessation to address lifestyle changes during the pandemic.

The genetic architecture of non-small cell lung cancer (NSCLC) is relevant to smoking status. However, the genetic contribution of long-term smoking cessation to the prognosis of NSCLC patients remains largely unknown. We conducted a genome-wide association study primarily on the prognosis of 1299 NSCLC patients of long-term former smokers from independent discovery (n = 566) and validation (n = 733) sets, and used in-silico function prediction and multi-omics analysis to identify single nucleotide polymorphisms (SNPs) on prognostics with NSCLC. We further detected SNPs with at least moderate association strength on survival within each group of never, short-term former, long-term former, and current smokers, and compared their genetic similarity at the SNP, gene, expression quantitative trait loci (eQTL), enhancer, and pathway levels. We identified two SNPs, rs34211819TNS3 at 7p12.3 (P = 3.90 × 10-9) and rs1143149SEPT7 at 7p14.2 (P = 9.75 × 10-9), were significantly associated with survival of NSCLC patients who were long-term former smokers. Both SNPs had significant interaction effects with years of smoking cessation (rs34211819TNS3: Pinteraction = 8.0 × 10-4; rs1143149SEPT7: Pinteraction = 0.003). In addition, in silico function prediction and multi-omics analysis provided evidence that these QTLs were associated with survival. Moreover, comparison analysis found higher genetic similarity between long-term former smokers and never-smokers, compared to short-term former smokers or current smokers. Pathway enrichment analysis indicated a unique pattern among long-term former smokers that was related to immune pathways. This study provides important insights into the genetic architecture associated with long-term former smoking NSCLC.

The X-chromosome is among the largest human chromosomes. It differs from autosomes by a number of important features including hemizygosity in males, an almost complete inactivation of one copy in females, and unique patterns of recombination. We used data from the Catalog of Published Genome Wide Association Studies to compare densities of the GWAS-detected SNPs on the X-chromosome and autosomes. The density of GWAS-detected SNPs on the X-chromosome is 6-fold lower compared to the density of the GWAS-detected SNPs on autosomes. Differences between the X-chromosome and autosomes cannot be explained by differences in the overall SNP density, lower X-chromosome coverage by genotyping platforms or low call rate of X-chromosomal SNPs. Similar differences in the density of GWAS-detected SNPs were found in female-only GWASs (e.g. ovarian cancer GWASs). We hypothesized that the lower density of GWAS-detected SNPs on the X-chromosome compared to autosomes is not a result of a methodological bias, e.g. differences in coverage or call rates, but has a real underlying biological reason-a lower density of functional SNPs on the X-chromosome versus autosomes. This hypothesis is supported by the observation that (i) the overall SNP density of X-chromosome is lower compared to the SNP density on autosomes and that (ii) the density of genic SNPs on the X-chromosome is lower compared to autosomes while densities of intergenic SNPs are similar.

Multi-omics datasets are becoming more common, necessitating better integration methods to realize their revolutionary potential. Here, we introduce multi-set correlation and factor analysis (MCFA), an unsupervised integration method tailored to the unique challenges of high-dimensional genomics data that enables fast inference of shared and private factors. We used MCFA to integrate methylation markers, protein expression, RNA expression, and metabolite levels in 614 diverse samples from the Trans-Omics for Precision Medicine/Multi-Ethnic Study of Atherosclerosis multi-omics pilot. Samples cluster strongly by ancestry in the shared space, even in the absence of genetic information, while private spaces frequently capture dataset-specific technical variation. Finally, we integrated genetic data by conducting a genome-wide association study (GWAS) of our inferred factors, observing that several factors are enriched for GWAS hits and trans-expression quantitative trait loci. Two of these factors appear to be related to metabolic disease. Our study provides a foundation and framework for further integrative analysis of ever larger multi-modal genomic datasets.

Blood lipid traits are treatable and heritable risk factors for heart disease, a leading cause of mortality worldwide. Although genome-wide association studies (GWAS) have discovered hundreds of variants associated with lipids in humans, most of the causal mechanisms of lipids remain unknown. To better understand the biological processes underlying lipid metabolism, we investigated the associations of plasma protein levels with total cholesterol (TC), triglycerides (TG), high-density lipoprotein cholesterol (HDL), and low-density lipoprotein cholesterol (LDL) in blood. We trained protein prediction models based on samples in the Multi-Ethnic Study of Atherosclerosis (MESA) and applied them to conduct proteome-wide association studies (PWAS) for lipids using the Global Lipids Genetics Consortium (GLGC) data. Of the 749 proteins tested, 42 were significantly associated with at least one lipid trait. Furthermore, we performed transcriptome-wide association studies (TWAS) for lipids using 9,714 gene expression prediction models trained on samples from peripheral blood mononuclear cells (PBMCs) in MESA and 49 tissues in the Genotype-Tissue Expression (GTEx) project. We found that although PWAS and TWAS can show different directions of associations in an individual gene, 40 out of 49 tissues showed a positive correlation between PWAS and TWAS signed p-values across all the genes, which suggests a high-level consistency between proteome-lipid associations and transcriptome-lipid associations.

Ocular choroidal melanoma (OCM) is the most common ocular primary malignant tumor in adults, and there is an increasing emphasis on its early detection and treatment worldwide. The main obstacle in early detection of OCM is its overlapping clinical features with benign choroidal nevus. Thus, we propose ultrasound localization microscopy (ULM) based on the image deconvolution algorithm to assist the diagnosis of small OCM in early stages. Furthermore, we develop ultrasound (US) plane wave imaging based on three-frame difference algorithm to guide the placement of the probe on the field of view. A high-frequency Verasonics Vantage system and an L22-14v linear array transducer were used to perform experiments on both custom-made modules in vitro and a SD rat with ocular choroidal melanoma in vivo. The results demonstrate that our proposed deconvolution method implement more robust microbubble (MB) localization, reconstruction of microvasculature network in a finer grid and more precise flow velocity estimation. The excellent performance of US plane wave imaging was successfully validated on the flow phantom and in an in vivo OCM model. In the future, the super-resolution ULM, a critical complementary imaging modality, can provide doctors with conclusive suggestions for early diagnosis of OCM, which is significant for the treatment and prognosis of patients.

During January-August 2021, the Community Prevalence of SARS-CoV-2 Study used time/location sampling to recruit a cross-sectional, population-based cohort to estimate SARS-CoV-2 seroprevalence and nasal swab sample PCR positivity across 15 US communities. Survey-weighted estimates of SARS-CoV-2 infection and vaccine willingness among participants at each site were compared within demographic groups by using linear regression models with inverse variance weighting. Among 22,284 persons >2 months of age and older, median prevalence of infection (prior, active, or both) was 12.9% across sites and similar across age groups. Within each site, average prevalence of infection was 3 percentage points higher for Black than White persons and average vaccine willingness was 10 percentage points lower for Black than White persons and 7 percentage points lower for Black persons than for persons in other racial groups. The higher prevalence of SARS-CoV-2 infection among groups with lower vaccine willingness highlights the disparate effect of COVID-19 and its complications.

Juvenile polyposis (JP) is an autosomal dominant syndrome in which affected patients develop upper- and/or lower-gastrointestinal (GI) polyps. A subset of families with JP have germline mutations in the SMAD4 (MADH4) gene and are at increased risk of GI cancers. To date, six families with JP have been described as having the same SMAD4 deletion (1244-1247delAGAC). The objective of the present study is to determine whether this deletion is a common ancestral mutation or a mutational hotspot. DNA from members of four families with JP, from Iowa, Mississippi, Texas, and Finland, that had this 4-bp deletion was used to genotype 15 simple tandem repeat polymorphism (STRP) markers flanking the SMAD4 gene, including 2 new STRPs within 6.3 and 70.9 kb of the deletion. Haplotypes cosegregating with JP in each family were constructed, and the distances of the closest markers were determined from the draft sequence of the human genome. No common haplotype was observed in these four families with JP. A 14-bp region containing the deletion had four direct repeats and one inverted repeat. Because no common ancestor was suggested by haplotype analysis and the sequence flanking the deletion contains repeats frequently associated with microdeletions, this common SMAD4 deletion in JP most likely represents a mutational hotspot.