Negative Genetic correlation in Common Factor GWAS

45 views
Skip to first unread message

Yu - Yvonne Huang

unread,
Feb 8, 2024, 12:38:11 AMFeb 8
to Genomic SEM Users
Dear All,

I'm running a CommonFactorGWAS for 4 related traits. The LDSC result shows some traits have negative genetic correlations with others.

Genetic Correlation between CRAE_UKB and AVR_UKB: 0.3179 (0.1664)     
Genetic Correlation between CRAE_UKB and SBP_EVAN: -0.2648 (0.0462)     
Genetic Correlation between CRAE_UKB and DBP_EVAN: -0.3031 (0.0468)     
Genetic Correlation between AVR_UKB and SBP_EVAN: -0.3144 (0.0603)     
Genetic Correlation between AVR_UKB and DBP_EVAN: -0.346 (0.0607)     
Genetic Correlation between SBP_EVAN and DBP_EVAN: 0.8099 (0.0392)


 I run the CFA model with one factor, and the model fits well:

CommonFactor_ML
$modelfit
      chisq df    p_chisq      AIC       CFI       SRMR
df 5.576171  2 0.06153892 21.57617 0.9995867 0.06217045

$results
       lhs op      rhs Unstandardized_Estimate Unstandardized_SE
1       F1 =~ CRAE_UKB              0.11715047       0.017421367
2       F1 =~  AVR_UKB              0.09745155       0.016105531
3       F1 =~ SBP_EVAN             -0.32194645       0.015622345
4       F1 =~ DBP_EVAN             -0.34598825       0.018050336
6 CRAE_UKB ~~ CRAE_UKB              0.11487112       0.018135033
7  AVR_UKB ~~  AVR_UKB              0.05827442       0.016975698
8 SBP_EVAN ~~ SBP_EVAN              0.03626659       0.009312294
9 DBP_EVAN ~~ DBP_EVAN              0.01589492       0.010970629

  Standardized_Est Standardized_SE      p_value
1        0.3266867      0.04858135 1.761628e-11
2        0.3743399      0.06186605 1.441171e-09
3       -0.8606960      0.04176498 2.322930e-94
4       -0.9395655      0.04901747 6.853501e-82
6        0.8932758      0.14102402 2.385594e-10
7        0.8598696      0.25048534 5.973511e-04
8        0.2592024      0.06655626 9.840754e-05
9        0.1172168      0.08090267 1.473763e-01

Then I ran:
 VESSELWIDTH_factor <- commonfactorGWAS(covstruc = LDSCoutput, SNPs = p_sumstats1, estimation = "ML", cores = 20, parallel = TRUE,smooth_check=TRUE)

In independent GWAS, each of the traits has several significant hits. However, the CommonFactorGWAS (F1 ~ SNP) did not give any significant hits....

Can this happen?  I assume extracting the common factor from different traits will be similar to conducting MTAG, which will boost the power of GWAS...But in my case, it didn't...

Best wishes,
Yvonne 
Reply all
Reply to author
Forward
0 new messages