Personal Genomics Skill v4.2.0
Comprehensive local DNA analysis with 1600+ markers across 30 categories. Privacy-first genetic analysis for AI agents.
Quick Start
python comprehensive_analysis.py /path/to/dna_file.txt
Triggers
Activate this skill when user mentions:
- DNA analysis, genetic analysis, genome analysis
- 23andMe, AncestryDNA, MyHeritage results
- Pharmacogenomics, drug-gene interactions
- Medication interactions, drug safety
- Genetic risk, disease risk, health risk
- Carrier status, carrier testing
- VCF file analysis
- APOE, MTHFR, CYP2D6, BRCA, or other gene names
- Polygenic risk scores
- Haplogroups, maternal lineage, paternal lineage
- Ancestry composition, ethnicity
- Hereditary cancer, Lynch syndrome
- Autoimmune genetics, HLA, celiac
- Pain sensitivity, opioid response
- Sleep optimization, chronotype, caffeine metabolism
- Dietary genetics, lactose intolerance, celiac
- Athletic genetics, sports performance
- UV sensitivity, skin type, melanoma risk
- Telomere length, longevity genetics
Supported Files
- 23andMe, AncestryDNA, MyHeritage, FTDNA
- VCF files (whole genome/exome, .vcf or .vcf.gz)
- Any tab-delimited rsid format
Output Location
~/dna-analysis/reports/
agent_summary.json - AI-optimized, priority-sorted
full_analysis.json - Complete data
report.txt - Human-readable
genetic_report.pdf - Professional PDF report
New v4.0 Features
Haplogroup Analysis
- Mitochondrial DNA (mtDNA) - maternal lineage
- Y-chromosome - paternal lineage (males only)
- Migration history context
- PhyloTree/ISOGG standards
Ancestry Composition
- Population comparisons (EUR, AFR, EAS, SAS, AMR)
- Admixture detection
- Ancestry informative markers
Hereditary Cancer Panel
- BRCA1/BRCA2 comprehensive
- Lynch syndrome (MLH1, MSH2, MSH6, PMS2)
- Other genes (APC, TP53, CHEK2, PALB2, ATM)
- ACMG-style classification
Autoimmune HLA
- Celiac (DQ2/DQ8) - can rule out if negative
- Type 1 Diabetes
- Ankylosing spondylitis (HLA-B27)
- Rheumatoid arthritis, lupus, MS
Pain Sensitivity
- COMT Val158Met
- OPRM1 opioid receptor
- SCN9A pain signaling
- TRPV1 capsaicin sensitivity
- Migraine susceptibility
PDF Reports
- Professional format
- Physician-shareable
- Executive summary
- Detailed findings
- Disclaimers included
New v4.1.0 Features
Medication Interaction Checker
from markers.medication_interactions import check_medication_interactions
result = check_medication_interactions(
medications=["warfarin", "clopidogrel", "omeprazole"],
genotypes=user_genotypes
)
# Returns critical/serious/moderate interactions with alternatives
- Accepts brand or generic names
- CPIC guidelines integrated
- PubMed citations included
- FDA warning flags
Sleep Optimization Profile
from markers.sleep_optimization import generate_sleep_profile
profile = generate_sleep_profile(genotypes)
# Returns ideal wake/sleep times, coffee cutoff, etc.
- Chronotype (morning/evening preference)
- Caffeine metabolism speed
- Personalized timing recommendations
Dietary Interaction Matrix
from markers.dietary_interactions import analyze_dietary_interactions
diet = analyze_dietary_interactions(genotypes)
# Returns food-specific guidance
- Caffeine, alcohol, saturated fat, lactose, gluten
- APOE-specific diet recommendations
- Bitter taste perception
Athletic Performance Profile
from markers.athletic_profile import calculate_athletic_profile
profile = calculate_athletic_profile(genotypes)
# Returns power/endurance type, recovery profile, injury risk
- Sport suitability scoring
- Training recommendations
- Injury prevention guidance
UV Sensitivity Calculator
from markers.uv_sensitivity import generate_uv_sensitivity_report
uv = generate_uv_sensitivity_report(genotypes)
# Returns skin type, SPF recommendation, melanoma risk
- Fitzpatrick skin type estimation
- Vitamin D synthesis capacity
- Melanoma risk factors
Natural Language Explanations
from markers.explanations import generate_plain_english_explanation
explanation = generate_plain_english_explanation(
rsid="rs3892097", gene="CYP2D6", genotype="GA",
trait="Drug metabolism", finding="Poor metabolizer carrier"
)
- Plain-English summaries
- Research variant flagging
- PubMed links
Telomere & Longevity
from markers.advanced_genetics import estimate_telomere_length
telomere = estimate_telomere_length(genotypes)
# Returns relative estimate with appropriate caveats
- TERT, TERC, OBFC1 variants
- Longevity associations (FOXO3, APOE)
Data Quality
- Call rate analysis
- Platform detection
- Confidence scoring
- Quality warnings
Export Formats
- Genetic counselor clinical export
- Apple Health compatible
- API-ready JSON
- Integration hooks
Marker Categories (21 total)
- Pharmacogenomics (159) - Drug metabolism
- Polygenic Risk Scores (277) - Disease risk
- Carrier Status (181) - Recessive carriers
- Health Risks (233) - Disease susceptibility
- Traits (163) - Physical/behavioral
- Haplogroups (44) - Lineage markers
- Ancestry (124) - Population informative
- Hereditary Cancer (41) - BRCA, Lynch, etc.
- Autoimmune HLA (31) - HLA associations
- Pain Sensitivity (20) - Pain/opioid response
- Rare Diseases (29) - Rare conditions
- Mental Health (25) - Psychiatric genetics
- Dermatology (37) - Skin and hair
- Vision & Hearing (33) - Sensory genetics
- Fertility (31) - Reproductive health
- Nutrition (34) - Nutrigenomics
- Fitness (30) - Athletic performance
- Neurogenetics (28) - Cognition/behavior
- Longevity (30) - Aging markers
- Immunity (43) - HLA and immune
- Ancestry AIMs (24) - Admixture markers
Agent Integration
The agent_summary.json provides:
{
"critical_alerts": [],
"high_priority": [],
"medium_priority": [],
"pharmacogenomics_alerts": [],
"apoe_status": {},
"polygenic_risk_scores": {},
"haplogroups": {
"mtDNA": {"haplogroup": "H", "lineage": "maternal"},
"Y_DNA": {"haplogroup": "R1b", "lineage": "paternal"}
},
"ancestry": {
"composition": {},
"admixture": {}
},
"hereditary_cancer": {},
"autoimmune_risk": {},
"pain_sensitivity": {},
"lifestyle_recommendations": {
"diet": [],
"exercise": [],
"supplements": [],
"avoid": []
},
"drug_interaction_matrix": {},
"data_quality": {}
}
Critical Findings (Always Alert User)
Pharmacogenomics
- DPYD variants - 5-FU/capecitabine FATAL toxicity risk
- HLA-B*5701 - Abacavir hypersensitivity
- HLA-B*1502 - Carbamazepine SJS (certain populations)
- MT-RNR1 - Aminoglycoside-induced deafness
Hereditary Cancer
- BRCA1/BRCA2 pathogenic - Breast/ovarian cancer syndrome
- Lynch syndrome genes - Colorectal/endometrial cancer
- TP53 pathogenic - Li-Fraumeni syndrome (multi-cancer)
Disease Risk
- APOE ε4/ε4 - ~12x Alzheimer's risk
- Factor V Leiden - Thrombosis risk, contraceptive implications
- HLA-B27 - Ankylosing spondylitis susceptibility (OR ~70)
Carrier Status
- CFTR - Cystic fibrosis (1 in 25 Europeans)
- HBB - Sickle cell (1 in 12 African Americans)
- HEXA - Tay-Sachs (1 in 30 Ashkenazi Jews)
Usage Examples
Basic Analysis
from comprehensive_analysis import main
main() # Uses command line args
Haplogroup Analysis
from markers.haplogroups import analyze_haplogroups
result = analyze_haplogroups(genotypes)
print(result["mtDNA"]["haplogroup"]) # e.g., "H"
Ancestry
from markers.ancestry_composition import get_ancestry_summary
ancestry = get_ancestry_summary(genotypes)
Cancer Panel
from markers.cancer_panel import analyze_cancer_panel
cancer = analyze_cancer_panel(genotypes)
if cancer["pathogenic_variants"]:
print("ALERT: Pathogenic variants detected")
Generate PDF
from pdf_report import generate_pdf_report
pdf_path = generate_pdf_report(analysis_results)
Export for Genetic Counselor
from exports import generate_genetic_counselor_export
clinical = generate_genetic_counselor_export(results, "clinical.json")
Privacy
- All analysis runs locally
- Zero network requests
- No data leaves the machine
Limitations
- Consumer arrays miss rare variants (~0.1% of genome)
- Results are probabilistic, not deterministic
- Not a medical diagnosis
- Most conditions 50-80% non-genetic
- Consult healthcare providers for medical decisions
- Negative hereditary cancer result does NOT rule out cancer syndrome
- Haplogroup resolution limited without WGS
When to Recommend Genetic Counseling
- Any pathogenic hereditary cancer variant
- APOE ε4/ε4 genotype
- Multiple critical pharmacogenomic findings
- Carrier status with reproduction implications
- High-risk autoimmune HLA types with symptoms
- Results causing significant user distress
1---2name: personal-genomics-skill-v4-2-03description: Comprehensive local DNA analysis with 1600+ markers across 30 categories. Privacy-first genetic analysis for AI agents.4---5
6# Personal Genomics Skill v4.2.0
7
8Comprehensive local DNA analysis with **1600+ markers** across **30 categories**. Privacy-first genetic analysis for AI agents.
9
10## Quick Start
11
12```bash
13python comprehensive_analysis.py /path/to/dna_file.txt
14```
15
16## Triggers
17
18Activate this skill when user mentions:
19- DNA analysis, genetic analysis, genome analysis
20- 23andMe, AncestryDNA, MyHeritage results
21- Pharmacogenomics, drug-gene interactions
22- Medication interactions, drug safety
23- Genetic risk, disease risk, health risk
24- Carrier status, carrier testing
25- VCF file analysis
26- APOE, MTHFR, CYP2D6, BRCA, or other gene names
27- Polygenic risk scores
28- Haplogroups, maternal lineage, paternal lineage
29- Ancestry composition, ethnicity
30- Hereditary cancer, Lynch syndrome
31- Autoimmune genetics, HLA, celiac
32- Pain sensitivity, opioid response
33- Sleep optimization, chronotype, caffeine metabolism
34- Dietary genetics, lactose intolerance, celiac
35- Athletic genetics, sports performance
36- UV sensitivity, skin type, melanoma risk
37- Telomere length, longevity genetics
38
39## Supported Files
40
41- 23andMe, AncestryDNA, MyHeritage, FTDNA
42- VCF files (whole genome/exome, .vcf or .vcf.gz)
43- Any tab-delimited rsid format
44
45## Output Location
46
47`~/dna-analysis/reports/`
48
49- `agent_summary.json` - AI-optimized, priority-sorted
50- `full_analysis.json` - Complete data
51- `report.txt` - Human-readable
52- `genetic_report.pdf` - Professional PDF report
53
54## New v4.0 Features
55
56### Haplogroup Analysis
57- Mitochondrial DNA (mtDNA) - maternal lineage
58- Y-chromosome - paternal lineage (males only)
59- Migration history context
60- PhyloTree/ISOGG standards
61
62### Ancestry Composition
63- Population comparisons (EUR, AFR, EAS, SAS, AMR)
64- Admixture detection
65- Ancestry informative markers
66
67### Hereditary Cancer Panel
68- BRCA1/BRCA2 comprehensive
69- Lynch syndrome (MLH1, MSH2, MSH6, PMS2)
70- Other genes (APC, TP53, CHEK2, PALB2, ATM)
71- ACMG-style classification
72
73### Autoimmune HLA
74- Celiac (DQ2/DQ8) - can rule out if negative
75- Type 1 Diabetes
76- Ankylosing spondylitis (HLA-B27)
77- Rheumatoid arthritis, lupus, MS
78
79### Pain Sensitivity
80- COMT Val158Met
81- OPRM1 opioid receptor
82- SCN9A pain signaling
83- TRPV1 capsaicin sensitivity
84- Migraine susceptibility
85
86### PDF Reports
87- Professional format
88- Physician-shareable
89- Executive summary
90- Detailed findings
91- Disclaimers included
92
93## New v4.1.0 Features
94
95### Medication Interaction Checker
96```python
97from markers.medication_interactions import check_medication_interactions
98
99result = check_medication_interactions(
100 medications=["warfarin", "clopidogrel", "omeprazole"],
101 genotypes=user_genotypes
102)
103# Returns critical/serious/moderate interactions with alternatives
104```
105- Accepts brand or generic names
106- CPIC guidelines integrated
107- PubMed citations included
108- FDA warning flags
109
110### Sleep Optimization Profile
111```python
112from markers.sleep_optimization import generate_sleep_profile
113
114profile = generate_sleep_profile(genotypes)
115# Returns ideal wake/sleep times, coffee cutoff, etc.
116```
117- Chronotype (morning/evening preference)
118- Caffeine metabolism speed
119- Personalized timing recommendations
120
121### Dietary Interaction Matrix
122```python
123from markers.dietary_interactions import analyze_dietary_interactions
124
125diet = analyze_dietary_interactions(genotypes)
126# Returns food-specific guidance
127```
128- Caffeine, alcohol, saturated fat, lactose, gluten
129- APOE-specific diet recommendations
130- Bitter taste perception
131
132### Athletic Performance Profile
133```python
134from markers.athletic_profile import calculate_athletic_profile
135
136profile = calculate_athletic_profile(genotypes)
137# Returns power/endurance type, recovery profile, injury risk
138```
139- Sport suitability scoring
140- Training recommendations
141- Injury prevention guidance
142
143### UV Sensitivity Calculator
144```python
145from markers.uv_sensitivity import generate_uv_sensitivity_report
146
147uv = generate_uv_sensitivity_report(genotypes)
148# Returns skin type, SPF recommendation, melanoma risk
149```
150- Fitzpatrick skin type estimation
151- Vitamin D synthesis capacity
152- Melanoma risk factors
153
154### Natural Language Explanations
155```python
156from markers.explanations import generate_plain_english_explanation
157
158explanation = generate_plain_english_explanation(
159 rsid="rs3892097", gene="CYP2D6", genotype="GA",
160 trait="Drug metabolism", finding="Poor metabolizer carrier"
161)
162```
163- Plain-English summaries
164- Research variant flagging
165- PubMed links
166
167### Telomere & Longevity
168```python
169from markers.advanced_genetics import estimate_telomere_length
170
171telomere = estimate_telomere_length(genotypes)
172# Returns relative estimate with appropriate caveats
173```
174- TERT, TERC, OBFC1 variants
175- Longevity associations (FOXO3, APOE)
176
177### Data Quality
178- Call rate analysis
179- Platform detection
180- Confidence scoring
181- Quality warnings
182
183### Export Formats
184- Genetic counselor clinical export
185- Apple Health compatible
186- API-ready JSON
187- Integration hooks
188
189## Marker Categories (21 total)
190
1911. **Pharmacogenomics** (159) - Drug metabolism
1922. **Polygenic Risk Scores** (277) - Disease risk
1933. **Carrier Status** (181) - Recessive carriers
1944. **Health Risks** (233) - Disease susceptibility
1955. **Traits** (163) - Physical/behavioral
1966. **Haplogroups** (44) - Lineage markers
1977. **Ancestry** (124) - Population informative
1988. **Hereditary Cancer** (41) - BRCA, Lynch, etc.
1999. **Autoimmune HLA** (31) - HLA associations
20010. **Pain Sensitivity** (20) - Pain/opioid response
20111. **Rare Diseases** (29) - Rare conditions
20212. **Mental Health** (25) - Psychiatric genetics
20313. **Dermatology** (37) - Skin and hair
20414. **Vision & Hearing** (33) - Sensory genetics
20515. **Fertility** (31) - Reproductive health
20616. **Nutrition** (34) - Nutrigenomics
20717. **Fitness** (30) - Athletic performance
20818. **Neurogenetics** (28) - Cognition/behavior
20919. **Longevity** (30) - Aging markers
21020. **Immunity** (43) - HLA and immune
21121. **Ancestry AIMs** (24) - Admixture markers
212
213## Agent Integration
214
215The `agent_summary.json` provides:
216
217```json
218{
219 "critical_alerts": [],
220 "high_priority": [],
221 "medium_priority": [],
222 "pharmacogenomics_alerts": [],
223 "apoe_status": {},
224 "polygenic_risk_scores": {},
225 "haplogroups": {
226 "mtDNA": {"haplogroup": "H", "lineage": "maternal"},
227 "Y_DNA": {"haplogroup": "R1b", "lineage": "paternal"}
228 },
229 "ancestry": {
230 "composition": {},
231 "admixture": {}
232 },
233 "hereditary_cancer": {},
234 "autoimmune_risk": {},
235 "pain_sensitivity": {},
236 "lifestyle_recommendations": {
237 "diet": [],
238 "exercise": [],
239 "supplements": [],
240 "avoid": []
241 },
242 "drug_interaction_matrix": {},
243 "data_quality": {}
244}
245```
246
247## Critical Findings (Always Alert User)
248
249### Pharmacogenomics
250- **DPYD** variants - 5-FU/capecitabine FATAL toxicity risk
251- **HLA-B*5701** - Abacavir hypersensitivity
252- **HLA-B*1502** - Carbamazepine SJS (certain populations)
253- **MT-RNR1** - Aminoglycoside-induced deafness
254
255### Hereditary Cancer
256- **BRCA1/BRCA2** pathogenic - Breast/ovarian cancer syndrome
257- **Lynch syndrome** genes - Colorectal/endometrial cancer
258- **TP53** pathogenic - Li-Fraumeni syndrome (multi-cancer)
259
260### Disease Risk
261- **APOE ε4/ε4** - ~12x Alzheimer's risk
262- **Factor V Leiden** - Thrombosis risk, contraceptive implications
263- **HLA-B27** - Ankylosing spondylitis susceptibility (OR ~70)
264
265### Carrier Status
266- **CFTR** - Cystic fibrosis (1 in 25 Europeans)
267- **HBB** - Sickle cell (1 in 12 African Americans)
268- **HEXA** - Tay-Sachs (1 in 30 Ashkenazi Jews)
269
270## Usage Examples
271
272### Basic Analysis
273```python
274from comprehensive_analysis import main
275main() # Uses command line args
276```
277
278### Haplogroup Analysis
279```python
280from markers.haplogroups import analyze_haplogroups
281result = analyze_haplogroups(genotypes)
282print(result["mtDNA"]["haplogroup"]) # e.g., "H"
283```
284
285### Ancestry
286```python
287from markers.ancestry_composition import get_ancestry_summary
288ancestry = get_ancestry_summary(genotypes)
289```
290
291### Cancer Panel
292```python
293from markers.cancer_panel import analyze_cancer_panel
294cancer = analyze_cancer_panel(genotypes)
295if cancer["pathogenic_variants"]:
296 print("ALERT: Pathogenic variants detected")
297```
298
299### Generate PDF
300```python
301from pdf_report import generate_pdf_report
302pdf_path = generate_pdf_report(analysis_results)
303```
304
305### Export for Genetic Counselor
306```python
307from exports import generate_genetic_counselor_export
308clinical = generate_genetic_counselor_export(results, "clinical.json")
309```
310
311## Privacy
312
313- All analysis runs locally
314- Zero network requests
315- No data leaves the machine
316
317## Limitations
318
319- Consumer arrays miss rare variants (~0.1% of genome)
320- Results are probabilistic, not deterministic
321- Not a medical diagnosis
322- Most conditions 50-80% non-genetic
323- Consult healthcare providers for medical decisions
324- Negative hereditary cancer result does NOT rule out cancer syndrome
325- Haplogroup resolution limited without WGS
326
327## When to Recommend Genetic Counseling
328
329- Any pathogenic hereditary cancer variant
330- APOE ε4/ε4 genotype
331- Multiple critical pharmacogenomic findings
332- Carrier status with reproduction implications
333- High-risk autoimmune HLA types with symptoms
334- Results causing significant user distress