Rare copy number variants in males and females with childhood attention-deficit/hyperactivity disorder.

While childhood attention-deficit/hyperactivity disorder (ADHD) is more prevalent in males than females, genetic contributors to this effect have not been established. Here, we explore sex differences in the contribution of common and/or rare genetic variants to ADHD. Participants were from the Adolescent Brain and Cognitive Development study (N = 1253 youth meeting DSM-5 criteria for ADHD [mean age = 11.46 years [SD = 0.87]; 31% female] and 5577 unaffected individuals [mean age = 11.42 years [SD = 0.89]; 50% female], overall 66% White, non-Hispanic (WNH), 19% Black/African American, and 15% other races. Logistic regression tested for interactions between sex (defined genotypically) and both rare copy number variants (CNV) and polygenic (common variant) risk in association with ADHD. There was a significant interaction between sex and the presence of a CNV deletion larger than 200 kb, both in the entire cohort (ß = -0.74, CI = [-1.27 to -0.20], FDR-corrected p = 0.048) and, at nominal significance levels in the WNH ancestry subcohort (ß = -0.86, CI = [-1.51 to -0.20], p = 0.010). Additionally, the number of deleted genes interacted with sex in association with ADHD (whole cohort. ß = -0.13, CI = [-0.23 to -0.029], FDR-corrected p = 0.048; WNH. ß = -0.17, CI = [-0.29 to -0.050], FDR-corrected p = 0.044) as did the total length of CNV deletions (whole cohort. ß = -0.12, CI = [-0.19 to -0.044], FDR-corrected p = 0.028; WNH. ß = -0.17, CI = [-0.28 to -0.061], FDR-corrected p = 0.034). This sex effect was driven by increased odds of childhood ADHD for females but not males in the presence of CNV deletions. No similar sex effect was found for CNV duplications or polygenic risk scores. The association between CNV deletions and ADHD was partially mediated by measures of cognitive flexibility. In summary, CNV deletions were associated with increased odds for childhood ADHD in females, but not males.

Genome-wide significant risk loci for mood disorders in the Old Order Amish founder population.

Genome-wide association studies (GWAS) of mood disorders in large case-control cohorts have identified numerous risk loci, yet pathophysiological mechanisms remain elusive, primarily due to the very small effects of common variants. We sought to discover risk variants with larger effects by conducting a genome-wide association study of mood disorders in a founder population, the Old Order Amish (OOA, n = 1,672). Our analysis revealed four genome-wide significant risk loci, all of which were associated with >2-fold relative risk. Quantitative behavioral and neurocognitive assessments (n = 314) revealed effects of risk variants on sub-clinical depressive symptoms and information processing speed. Network analysis suggested that OOA-specific risk loci harbor novel risk-associated genes that interact with known neuropsychiatry-associated genes via gene interaction networks. Annotation of the variants at these risk loci revealed population-enriched, non-synonymous variants in two genes encoding neurodevelopmental transcription factors, CUX1 and CNOT1. Our findings provide insight into the genetic architecture of mood disorders and a substrate for mechanistic and clinical studies.

Mendelian Randomization Analysis Reveals a Complex Genetic Interplay among Atopic Dermatitis, Asthma, and Gastroesophageal Reflux Disease.

Rationale: Gastroesophageal reflux disease (GERD) is commonly associated with atopic disorders, but cause-effect relationships remain unclear. Objectives: We applied Mendelian randomization analysis to explore whether GERD is causally related to atopic disorders of the lung (asthma) and/or skin (atopic dermatitis [AD]). Methods: We conducted two-sample bidirectional Mendelian randomization to infer the magnitude and direction of causality between asthma and GERD, using summary statistics from the largest genome-wide association studies conducted on asthma (Ncases = 56,167) and GERD (Ncases = 71,522). In addition, we generated instrumental variables for AD from the latest population-level genome-wide association study meta-analysis (Ncases = 22,474) and assessed their fidelity and confidence of predicting the likely causal pathway(s) leading to asthma and/or GERD. Measurements and Main Results: Applying three different methods, each method revealed similar magnitude of causal estimates that were directionally consistent across the sensitivity analyses. Using an inverse variance-weighted method, the largest effect size was detected for asthma predisposition to AD (odds ratio [OR], 1.46; 95% confidence interval [CI], 1.34-1.59), followed by AD to asthma (OR, 1.34; 95% CI, 1.24-1.45). A significant association was detected for genetically determined asthma on risk of GERD (OR, 1.06; 95% CI, 1.03-1.09) but not genetically determined AD on GERD. In contrast, GERD equally increased risks of asthma (OR, 1.21; 95% CI, 1.09-1.35) and AD (OR, 1.21; 95% CI, 1.07-1.37). Conclusions: This study uncovers previously unrecognized causal pathways that have clinical implications in European-ancestry populations: 1) asthma is a causal risk for AD, and 2) the predisposition to AD, including asthma, can arise from specific pathogenic mechanisms manifested by GERD.

Transforming Growth Factor-ß1 Decreases ß2-Agonist-induced Relaxation in Human Airway Smooth Muscle.

Helper T effector cytokines implicated in asthma modulate the contractility of human airway smooth muscle (HASM) cells. We have reported recently that a profibrotic cytokine, transforming growth factor (TGF)-ß1, induces HASM cell shortening and airway hyperresponsiveness. Here, we assessed whether TGF-ß1 affects the ability of HASM cells to relax in response to ß2-agonists, a mainstay treatment for airway hyperresponsiveness in asthma. Overnight TGF-ß1 treatment significantly impaired isoproterenol (ISO)-induced relaxation of carbachol-stimulated, isolated HASM cells. This single-cell mechanical hyporesponsiveness to ISO was corroborated by sustained increases in myosin light chain phosphorylation. In TGF-ß1-treated HASM cells, ISO evoked markedly lower levels of intracellular cAMP. These attenuated cAMP levels were, in turn, restored with pharmacological and siRNA inhibition of phosphodiesterase 4 and Smad3, respectively. Most strikingly, TGF-ß1 selectively induced phosphodiesterase 4D gene expression in HASM cells in a Smad2/3-dependent manner. Together, these data suggest that TGF-ß1 decreases HASM cell ß2-agonist relaxation responses by modulating intracellular cAMP levels via a Smad2/3-dependent mechanism. Our findings further define the mechanisms underlying ß2-agonist hyporesponsiveness in asthma, and suggest TGF-ß1 as a potential therapeutic target to decrease asthma exacerbations in severe and treatment-resistant asthma.

15q13.3 duplication in two patients with childhood-onset schizophrenia.

We report two cases of paternally inherited 15q13.3 duplications in carriers diagnosed with childhood-onset schizophrenia (COS), a rare neurodevelopmental disorder of proposed polygenic origin with onset in children before age 13. This study documents that the 15q13.3 deletion and duplication exhibit pathogenicity for COS, with both copy number variants (CNVs) sharing a disrupted CHRNA7 gene. CHRNA7 encodes the neuronal alpha7 nicotinic acetylcholine receptor (α7nAChR) and is a candidate gene that has been suggested as a pathophysiological process mediating adult-onset schizophrenia (AOS) and other neurodevelopmental disorders. These results support the incomplete penetrance and variable expressivity of this CNV and represent the first report of 15q13.3 duplication carriers exhibiting COS. Published 2016. This article is a U.S. Government work and is in the public domain in the USA. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics published by Wiley Periodicals, Inc.

P2X7-regulated protection from exacerbations and loss of control is independent of asthma maintenance therapy.

RATIONALE: The function of the P2X(7) nucleotide receptor protects against exacerbation in people with mild-intermittent asthma during viral illnesses, but the impact of disease severity and maintenance therapy has not been studied. OBJECTIVES: To evaluate the association between P2X(7), asthma exacerbations, and incomplete symptom control in a more diverse population. METHODS: A matched P2RX7 genetic case-control was performed with samples from Asthma Clinical Research Network trial participants enrolled before July 2006, and P2X(7) pore activity was determined in whole blood samples as an ancillary study to two trials completed subsequently. MEASUREMENTS AND MAIN RESULTS: A total of 187 exacerbations were studied in 742 subjects, and the change in asthma symptom burden was studied in an additional 110 subjects during a trial of inhaled corticosteroids (ICS) dose optimization. African American carriers of the minor G allele of the rs2230911 loss-of-function single nucleotide polymorphism were more likely to have a history of prednisone use in the previous 12 months, with adjustment for ICS and long-acting ß(2)-agonists use (odds ratio, 2.7; 95% confidence interval, 1.2-6.2; P = 0.018). Despite medium-dose ICS, attenuated pore function predicted earlier exacerbations in incompletely controlled patients with moderate asthma (hazard ratio, 3.2; confidence interval, 1.1-9.3; P = 0.033). After establishing control with low-dose ICS in patients with mild asthma, those with attenuated pore function had more asthma symptoms, rescue albuterol use, and FEV(1) reversal (P < 0.001, 0.03, and 0.03, respectively) during the ICS adjustment phase. CONCLUSIONS: P2X(7) pore function protects against exacerbations of asthma and loss of control, independent of baseline severity and the maintenance therapy.

Cortico-striatal differences in the epigenome in attention-deficit/ hyperactivity disorder.

While epigenetic modifications have been implicated in ADHD through studies of peripheral tissue, to date there has been no examination of the epigenome of the brain in the disorder. To address this gap, we mapped the methylome of the caudate nucleus and anterior cingulate cortex in post-mortem tissue from fifty-eight individuals with or without ADHD. While no single probe showed adjusted significance in differential methylation, several differentially methylated regions emerged. These regions implicated genes involved in developmental processes including neurogenesis and the differentiation of oligodendrocytes and glial cells. We demonstrate a significant association between differentially methylated genes in the caudate and genes implicated by GWAS not only in ADHD but also in autistic spectrum, obsessive compulsive and bipolar affective disorders through GWAS. Using transcriptomic data available on the same subjects, we found modest correlations between the methylation and expression of genes. In conclusion, this study of the cortico-striatal methylome points to gene and gene pathways involved in neurodevelopment, consistent with studies of common and rare genetic variation, as well as the post-mortem transcriptome in ADHD.

Mechanisms of oral carcinogenesis induced by dibenzo[a,l]pyrene: an environmental pollutant and a tobacco smoke constituent.

We previously reported that dibenzo[a,l]pyrene (DB[a,l]P), the most potent known environmental carcinogen among polycyclic aromatic hydrocarbons (PAH) congeners, is carcinogenic in the oral tissues of mice. We have now developed a new mouse model which employs the oral application of the fjord region diol epoxide, (±)-anti-11,12-dihydroxy-13,14-epoxy-11,12,13,14-tetrahydrodibenzo[a,l]pyrene (DB[a,l]PDE), a metabolite of the tobacco smoke constituent DB[a,l]P, and we show its specific induction of oral squamous cell carcinoma (OSCC) in both tongue and other oral tissues. Groups of B6C3F1 mice (20/group) received 6 or 3 nmol of (±)-anti-DB[a,l]PDE administered into the oral cavity; 3 times per week for 38 weeks. Additional groups received the vehicle alone or were left untreated. Mice were sacrificed 42 weeks after the first carcinogen administration. The high dose induced 74 and 100% OSCC in the tongue and other oral tissues, respectively; the corresponding values at the lower dose were 45 and 89%. Using immunohistochemistry, we showed that DB[a,l]PDE resulted in overexpression of p53 and COX-2 proteins in malignant tissues when compared to normal oral tissues and tongues. Consistent with the carcinogenicity, we demonstrated powerful mutagenicity in cII gene in B6C3F1 (Big Blue) mouse tongue. The mutational profile in lacI reporter gene is similar to those detected in human head and neck cancer, and p53 mutations were observed in mouse oral tumor tissues. Taken together, we conclude that the formation of diol epoxides plays a major role among the mechanisms by which DB[a,l]P exerts its oral mutagenicity and tumorigenicity.

Longitudinal trajectories of childhood and adolescent attention deficit hyperactivity disorder diagnoses in three cohorts.

Background: Attention deficit/hyperactivity disorder (ADHD) is usually conceptualized as a childhood-onset neurodevelopmental disorder, in which symptoms either decrease steadily into adulthood or remain stable. A recent study challenged this view, reporting that for most with ADHD, diagnostic status fluctuates with age. We ask if such a 'fluctuating' ADHD symptom trajectory subgroup is present in other population-based and clinic-based cohorts, centered on childhood and adolescence. Methods: Cohorts were the population-based Adolescent Brain Cognitive Development (ABCD: N = 9735), Neurobehavioral Clinical Research (NCR: N = 258), and the Nathan Kline Institute-Rockland (NKI-Rockland: N = 149). All participants had three or more assessments spanning different age windows. Participants were categorized into developmental diagnostic subgroups: fluctuant ADHD (defined by two or more switches between meeting and not meeting ADHD criteria), remitting ADHD, persisting ADHD, emerging ADHD and never affected. Data were collected between 2011 and 2022. Analyses were performed between May 2022 and April 2023. Findings: A subgroup with fluctuant child and adolescent ADHD diagnoses was found in all cohorts (29.3% of participants with ADHD in ABCD, 26.6% in NCR and 17% in NKI-Rockland). While the proportion of those with fluctuant ADHD increased with the number of assessments, it never constituted the dominant subgroup. Interpretation: We provide further evidence in three cohorts for the existence of a fluctuant ADHD diagnostic subgroup during childhood and adolescence, albeit in a minority of cases. Such fluctuant child and adolescent ADHD diagnoses may suggest a natural history more akin to relapsing-remitting mood disorders and/or a marked sensitivity to environmental shifts that occur across development. Funding: Intramural programs of the NHGRI and NIMH.

TAS2R activation promotes airway smooth muscle relaxation despite ß(2)-adrenergic receptor tachyphylaxis.

Recently, bitter taste receptors (TAS2Rs) were found in the lung and act to relax airway smooth muscle (ASM) via intracellular Ca(2+) concentration signaling generated from restricted phospholipase C activation. As potential therapy, TAS2R agonists could be add-on treatment when patients fail to achieve adequate bronchodilation with chronic ß-agonists. The ß(2)-adrenergic receptor (ß(2)AR) of ASM undergoes extensive functional desensitization. It remains unknown whether this desensitization affects TAS2R function, by cross talk at the receptors or distal common components in the relaxation machinery. We studied intracellular signaling and cell mechanics using isolated human ASM, mouse tracheal responses, and human bronchial responses to characterize TAS2R relaxation in the context of ß(2)AR desensitization. In isolated human ASM, magnetic twisting cytometry revealed >90% loss of isoproterenol-promoted decrease in cell stiffness after 18-h exposure to albuterol. Under these same conditions of ß(2)AR desensitization, the TAS2R agonist chloroquine relaxation response was unaffected. TAS2R-mediated stimulation of intracellular Ca(2+) concentration in human ASM was unaltered by albuterol pretreatment, in contrast to cAMP signaling, which was desensitized by >90%. In mouse trachea, ß(2)AR desensitization by ß-agonist amounted to 92 ± 6.0% (P < 0.001), while, under these same conditions, TAS2R desensitization was not significant (11 ± 3.5%). In human lung slices, chronic ß-agonist exposure culminated in 64 ± 5.7% (P < 0.001) desensitization of ß(2)AR-mediated dilation of carbachol-constricted airways that was reversed by chloroquine. We conclude that there is no evidence for physiologically relevant cross-desensitization of TAS2R-mediated ASM relaxation from chronic ß-agonist treatment. These findings portend a favorable therapeutic profile for TAS2R agonists for the treatment of bronchospasm in asthma or chronic obstructive lung disease.

Mutagenesis and carcinogenesis induced by dibenzo[a,l]pyrene in the mouse oral cavity: a potential new model for oral cancer.

Cancer of the oral cavity is a serious disease, affecting about 30,000 individuals in US annually. There are several animal models of oral cancer, but each has certain disadvantages. As a new model, we investigated whether topical application of the tobacco smoke carcinogen, dibenzo[a,l]pyrene (DB[a,l]P) is mutagenic and carcinogenic in the oral cavity of the B6C3F1 lacI and B6C3F1 mouse, respectively. B6C3F1 lacI mice received DB[a,l]P (0, 3, 6, 12 nmol) 3× per week. B6C3F1 mice received the same doses and also 24 nmol. At 38 weeks mutagenesis was measured in oral tissues in lacI mice. For the high dose group, the mutant fraction (MF) in upper mucosa and tongue increased about twofold relative to that in vehicle-alone. The increases were statistically significant. The mutational profile in the DB[a,l]P-induced mutants was compared with that induced by benzo[a]pyrene (BaP) in oral tissue. BaP is mutagenic in many tissues when administered by gavage. The mutational profile for DB[a,l]P was more similar to that reported for p53 mutations in head and neck cancers than was that of BaP. At 47 weeks, oral squamous cell carcinomas (OSCC) were found in 31% of the high-dose B6C3F1 group. Elevations of p53 and COX-2 protein were observed in tumor and dysplastic tissue. As DB[a,l]P induces mutations and tumors in the oral cavity, and has a mutational profile in oral tissue similar to that found in p53 in human OSCC, the treatment protocol described here may represent a new and relevant model for cancer of the oral cavity.

ADHD and its neurocognitive substrates: A two sample Mendelian randomization study.

Attention-deficit/hyperactivity disorder (ADHD) is associated with a wide array of neural and cognitive features, and other psychiatric disorders, identified mainly through cross-sectional associations studies. However, it is unclear if the disorder is causally associated with these neurocognitive features. Here, we applied a two-sample bidirectional Mendelian randomization (MR) study to summary GWAS data to explore the presence and direction of a causal effect between ADHD and a range of neurocognitive features and other psychiatric disorders. The inverse variance weighted method was used in the main analysis, and two MR methods (MR-Egger, weighted median) were used for robustness checks. We found that genetic risk for ADHD was causally associated with a decreased area of lateral orbitofrontal cortex. Conversely, we found that brain volume and some features of intrinsic functional connectivity had causal effects on ADHD risk. Bidirectional causal links were found between ADHD and adult general intelligence, as well as depression and autistic spectrum disorders. Such work highlights the important ties between ADHD and general cognitive ability, and suggest some neural features, previously merely associated with the disorder, may play a causal role in its pathogenesis.

Nicotine dependence phenotype, time to first cigarette, and risk of head and neck cancer.

BACKGROUND: A behavioral phenotype that characterizes nicotine dependence, the time to first cigarette after waking, is hypothesized to increase the risk of head and neck cancer. METHODS: A case-control study of histologically confirmed head and neck cancer was conducted that included 1055 cases and 795 controls with a history of cigarette smoking. RESULTS: The pack-years-adjusted odds ratio was 1.42 (95% confidence interval [95% CI], 1.02-1.99) for an interval of 31 minutes to 60 minutes to first cigarette after waking and 1.59 (95% CI, 1.19-2.11) for an interval of 1 minute to 30 minutes. The risk estimates were similar when smoking was modeled as total years, smoking status (current vs former), number of cigarettes smoked per day, years since quitting, and excess odds ratio. Findings were consistent for cancers of the floor of the mouth, palate, and pharynx. CONCLUSIONS: Time to first cigarette is an indicator of increased nicotine dependence, smoke uptake, and risk of head and neck cancer. This high-risk group of individuals would benefit from targeted smoking interventions.

Nicotine dependence phenotype and lung cancer risk.

BACKGROUND: A behavioral phenotype that characterizes nicotine dependence, the time to first cigarette after waking, is hypothesized to increase the risk of lung cancer. METHODS: A case-control study of histologically confirmed lung cancer was conducted. The current analysis included 4775 lung cancer cases and 2835 controls who were regular cigarette smokers. RESULTS: Compared with subjects who smoked their first cigarette > 60 minutes after waking, the pack-years-adjusted odds ratio was 1.31 (95% confidence interval [95% CI], 1.11-1.54) for subjects who smoked 31 minutes to 60 minutes after waking and 1.79 (95% CI, 1.56-2.07) for subjects who smoked within 30 minutes of waking. The risk estimates were similar when smoking was modeled as total years, smoking status (current vs former), number of cigarettes smoked per day, years since quitting, and excess odds ratio. The findings were consistent for all histologic types of lung cancer. CONCLUSIONS: The findings of the current study indicate that a specific nicotine dependence phenotype that is associated with the amount of smoke uptake per cigarette is independently associated with lung cancer risk. These findings may help to identify high-risk individuals who would benefit from targeted interventions.

A dynamic model for genome-wide association studies.

Although genome-wide association studies (GWAS) are widely used to identify the genetic and environmental etiology of a trait, several key issues related to their statistical power and biological relevance have remained unexplored. Here, we describe a novel statistical approach, called functional GWAS or fGWAS, to analyze the genetic control of traits by integrating biological principles of trait formation into the GWAS framework through mathematical and statistical bridges. fGWAS can address many fundamental questions, such as the patterns of genetic control over development, the duration of genetic effects, as well as what causes developmental trajectories to change or stop changing. In statistics, fGWAS displays increased power for gene detection by capitalizing on cumulative phenotypic variation in a longitudinal trait over time and increased robustness for manipulating sparse longitudinal data.

A novel small molecule target in human airway smooth muscle for potential treatment of obstructive lung diseases: a staged high-throughput biophysical screening.

BACKGROUND: A newly identified mechanism of smooth muscle relaxation is the interaction between the small heat shock protein 20 (HSP20) and 14-3-3 proteins. Focusing upon this class of interactions, we describe here a novel drug target screening approach for treating airflow obstruction in asthma. METHODS: Using a high-throughput fluorescence polarization (FP) assay, we screened a library of compounds that could act as small molecule modulators of HSP20 signals. We then applied two quantitative, cell-based biophysical methods to assess the functional efficacy of these molecules and rank-ordered their abilities to relax isolated human airway smooth muscle (ASM). Scaling up to the level of an intact tissue, we confirmed in a concentration-responsive manner the potency of the cell-based hit compounds. RESULTS: Among 58,019 compound tested, 268 compounds caused 20% or more reduction of the polarized emission in the FP assay. A small subset of these primary screen hits, belonging to two scaffolds, caused relaxation of isolated ASM cell in vitro and attenuated active force development of intact tissue ex vivo. CONCLUSIONS: This staged biophysical screening paradigm provides proof-of-principle for high-throughput and cost-effective discovery of new small molecule therapeutic agents for obstructive lung diseases.

A disequilibrium model for detecting genetic mutations for cancer.

It has been recognized that genetic mutations in specific nucleotides may give rise to cancer via the alteration of signaling pathways. Thus, the detection of those cancer-causing mutations has received considerable interest in cancer genetic research. Here, we propose a statistical model for characterizing genes that lead to cancer through point mutations using genome-wide single nucleotide polymorphism (SNP) data. The basic idea of the model is that mutated genes may be in high association with their nearby SNPs because of evolutionary forces. By genotyping SNPs in both normal and cancer cells, we formulate a polynomial likelihood to estimate the population genetic parameters related to cancer, such as allele frequencies of cancer-causing alleles, mutation rates of alleles derived from maternal or paternal parents, and zygotic linkage disequilibria between different loci after the mutation occurs. We implement the EM algorithm to estimate some of these parameters because of the missing information in the likelihood construction. The model allows the elegant tests of the significant associations between mutated cancer genes and genome-wide SNPs, thus providing a way for predicting the occurrence and formation of cancer with genetic information. The model, validated through computer simulation, may help cancer geneticists design efficient experiments and formulate hypotheses for cancer gene identification.

Increase of rejection rate in case-control studies with the differential genotyping error rates.

Genotyping error adversely affects the statistical power of case-control association studies and introduces bias in the estimated parameters when the same error mechanism and probabilities apply to both affected and unaffected individuals; that is, when there is non-differential genotype misclassification. Simulation studies have shown that differential genotype misclassification leads to a rejection rate that is higher than the nominal significance level (type I error rate) for some tests of association. This study extends previous work by examining this issue analytically using the non-centrality parameter of the asymptotic distribution of the chi-squared test and linear trend test (LTT) when there is no difference between case and control genotype frequencies, but there is differential misclassification with SNP data. The parameters examined are the minor allele frequency (MAF) and sample size. When MAF is less than 0.2, differential genotyping errors lead to a rejection rate much larger than the nominal significance level. As the MAF decreases to zero, the increase in the rejection rate becomes larger. The errors that most increase the rejection rate are differential recording of the more common homozygote as the other homozygote and differential recording of the more common homozygote as the heterozygote. The rejection rate increases as the sample size increases for fixed differential genotyping error rates and nominal significance level for each test.

A statistical model for genetic mapping of viral infection by integrating epidemiological behavior.

Large-scale studies of genetic variation may be helpful for understanding the genetic control mechanisms of viral infection and, ultimately, predicting and eliminating infectious disease outbreaks. We propose a new statistical model for detecting specific DNA sequence variants that are responsible for viral infection. This model considers additive, dominance and epistatic effects of haplotypes from three different genomes, recipient, transmitter and virus, through an epidemiological process. The model is constructed within the maximum likelihood framework and implemented with the EM algorithm. A number of hypothesis tests about population genetic structure and diversity and the pattern of genetic control are formulated. A series of closed forms for the EM algorithm to estimate haplotype frequencies and haplotype effects in a network of genetic interactions among three genomes are derived. Simulation studies were performed to test the statistical properties of the model, recommending necessary sample sizes for obtaining reasonably good accuracy and precision of parameter estimation. By integrating, for the first time, the epidemiological principle of viral infection into genetic mapping, the new model shall find an immediate application to studying the genetic architecture of viral infection.