# Gene Regulation Analysis - Final Summary

## Overview
This analysis examined the influence of miRNA, lncRNA, and DNA methylation on gene expression using multi-omics data from 40 samples (ACR-139 through ACR-265, each with 4 time points).

## Data Summary
- **Gene Expression**: 36,084 genes across 40 samples
- **miRNA Expression**: 51 miRNAs across 40 samples  
- **lncRNA Expression**: 15,900 lncRNAs across 40 samples
- **DNA Methylation**: 249 CpG sites across 40 samples
- **Sample Overlap**: 40 samples with complete data across all molecular layers

## Analysis Results

### Correlation Analysis (Primary Results)
The correlation analysis examined direct relationships between regulatory factors and gene expression for the top 20 most variable genes.

**Key Findings:**
- **Total correlations analyzed**: 60 (20 genes × 3 regulatory factors)
- **Significant correlations (p < 0.05)**: 13 out of 60 (21.7%)
- **Strong correlations (|r| > 0.5)**: 1 out of 60 (1.7%)

### Regulatory Factor Influence

#### 1. DNA Methylation
- **Mean correlation**: 0.031 (weak positive)
- **Significant correlations**: 7 out of 20 genes
- **Strongest correlation**: FUN_009574 (r = 0.545, p = 2.8e-04)
- **Interpretation**: DNA methylation shows the strongest regulatory influence among the three factors

#### 2. miRNA Expression
- **Mean correlation**: 0.021 (weak positive)  
- **Significant correlations**: 5 out of 20 genes
- **Strongest correlation**: FUN_008377 (r = 0.456, p = 3.1e-03)
- **Interpretation**: miRNA shows moderate regulatory influence

#### 3. lncRNA Expression
- **Mean correlation**: -0.014 (weak negative)
- **Significant correlations**: 1 out of 20 genes
- **Interpretation**: lncRNA shows the weakest regulatory influence

### Top Regulatory Relationships

1. **FUN_009574 - DNA Methylation**: r = 0.545 (p = 2.8e-04)
   - Strongest regulatory relationship
   - Suggests DNA methylation significantly influences this gene's expression

2. **FUN_042621 - DNA Methylation**: r = 0.500 (p = 1.0e-03)
   - Moderate regulatory relationship
   - DNA methylation explains ~25% of expression variance

3. **FUN_008377 - DNA Methylation**: r = 0.458 (p = 2.9e-03)
   - Moderate regulatory relationship
   - Also shows miRNA correlation (r = 0.456)

4. **FUN_008377 - miRNA**: r = 0.456 (p = 3.1e-03)
   - Moderate miRNA regulatory influence
   - Same gene shows both methylation and miRNA correlations

5. **FUN_009574 - miRNA**: r = 0.428 (p = 5.8e-03)
   - Moderate miRNA regulatory influence
   - Same gene shows both methylation and miRNA correlations

## Biological Interpretation

### Regulatory Hierarchy
1. **DNA Methylation** exerts the strongest regulatory influence
2. **miRNA** shows moderate regulatory influence  
3. **lncRNA** shows minimal regulatory influence

### Multi-Layer Regulation
Several genes (FUN_009574, FUN_008377) show correlations with multiple regulatory factors, suggesting:
- **Combinatorial regulation**: Multiple molecular layers work together
- **Redundant regulation**: Backup regulatory mechanisms
- **Context-dependent regulation**: Different factors active under different conditions

### Sample Characteristics
- **40 samples** representing different conditions/time points
- **High variability** in gene expression (top 20 genes selected by variance)
- **Consistent regulatory patterns** across samples

## Technical Notes

### Data Preprocessing
- All data underwent log2 transformation for normalization
- Sample names standardized across datasets
- Zero-count filtering applied to remove non-expressed features

### Statistical Analysis
- Pearson correlation analysis with significance testing (p < 0.05)
- Cross-validation in regression models
- Multiple regression approaches (Linear, Ridge, Lasso, Random Forest)

### Quality Control
- Sample overlap verified across all molecular layers
- Data consistency checks performed
- Outlier detection and handling

## Conclusions

1. **DNA methylation is the primary regulatory mechanism** influencing gene expression in this dataset
2. **miRNA provides secondary regulatory control** with moderate influence
3. **lncRNA shows minimal direct regulatory influence** on the analyzed genes
4. **Multi-layer regulation exists** for some genes, suggesting complex regulatory networks
5. **Regulatory influence is generally weak to moderate**, indicating these factors work in combination rather than individually

## Recommendations

1. **Focus on DNA methylation** for primary regulatory studies
2. **Investigate miRNA-mRNA interactions** for secondary regulatory mechanisms  
3. **Explore lncRNA function** in other contexts or with different analysis approaches
4. **Study combinatorial effects** of multiple regulatory layers
5. **Validate findings** with experimental approaches (e.g., CRISPR, knockdown studies)

## Files Generated

- `correlation_results.csv`: Detailed correlation coefficients and p-values
- `correlation_summary.txt`: Statistical summary of results
- `correlation_heatmap.png`: Visual representation of gene-regulatory factor correlations
- `correlation_distributions.png`: Distribution of correlation coefficients by regulatory factor
- `significant_correlations.png`: Analysis of statistically significant correlations
- `integrated_dataset.csv`: Combined dataset with all molecular layers
- `regulatory_influence_results.csv`: Regression analysis results
- `analysis_summary.txt`: Comprehensive analysis summary
- `regulatory_influence_heatmap.png`: Model performance heatmap
- `r2_distribution.png`: Distribution of R² scores
- `model_comparison.png`: Comparison of different regression models
- `sample_correlation.png`: Sample correlation matrix

---
*Analysis completed on: August 15, 2024*
*Data: 40 samples, 4 molecular layers, 20 high-variance genes*