title: Evaluation of linkage disequilibrium in wheat with an L1-regularized sparse Markov network.
abstract: Linkage disequilibrium (LD) is defined as a stochastic dependence between alleles at two or more loci. Although understanding LD is important in the study of the genetics of many species, little attention has been paid on how a covariance structure between many loci distributed across the genome should be represented. Given that biological systems at the cellular level often involve gene networks, it is appealing to evaluate LD from a network perspective, i.e., as a set of associated loci involved in a complex system. We applied a Markov network (MN) to study LD using data on 1,279 markers derived from 599 wheat inbred lines. The MN attempts to account for association between two markers, conditionally on the remaining markers in the network model. In this study, the recovery of the structure of a LD network was done through two variants of pseudo-likelihoods subject to an L1 penalty on the MN parameters. It is shown that, while the L1-regularized Markov network preserves features of a Bayesian network (BN), the nodes in the resulting networks have fewer links. The resulting sparse network, encoding conditional independencies, provides a clearer picture of association than marginal LD metrics, and a sparse graph eases interpretation markedly, since it includes a smaller number of edges than a BN. Thus, an L1-regularized sparse Markov network seems appealing for representing conditional LD with high-dimensional genomic data, where variables, e.g., single nucleotide polymorphism markers, are expected to be sparsely connected.
authors: Department of Animal Sciences, University of Wisconsin, Madison, WI 53706, USA. morota@ansci.wisc.edu
authors: Gota
authors: G
authors: Morota
authors: Daniel
authors: D
authors: Gianola
mesh_terms: D000483:Alleles; D001499:Bayes Theorem; D005787:Gene Frequency; D053263:Gene Regulatory Networks / Q000235:genetics*; D005819:Genetic Markers; D005838:Genotype; D015810:Linkage Disequilibrium / Q000235:genetics*; D008390:Markov Chains; D008957:Models, Genetic; D014908:Triticum / Q000235:genetics*
chemical_list: D005819:Genetic Markers
references: J Stat Softw. 2010;33(1):1-22
references: 20808728
references: J Mach Learn Res. 2009 Apr 1;10:883-906
references: 21857799
references: Biostatistics. 2008 Jul;9(3):432-41
references: 18079126
references: Plant Genome. 2010;3(2):106-116
references: 21566722
references: Theor Appl Genet. 1968 Jun;38(6):226-31
references: 24442307
references: Genetics. 1921 Mar;6(2):111-23
references: 17245958
references: Nat Biotechnol. 2006 Apr;24(4):427-33
references: 16601728
references: Genetics. 2001 Apr;157(4):1819-29
references: 11290733
references: Anim Genet. 2012 Jul;43 Suppl 1:19-35
references: 22742500
references: PLoS One. 2010 Dec 20;5(12):e14147
references: 21188141
references: Biometrics. 2011 Mar;67(1):164-73
references: 20377578
references: Genetics. 2000 May;155(1):421-30
references: 10790414
references: Biometrics. 1975 Jun;31(2):423-47
references: 1174616
references: Genetics. 2009 May;182(1):375-85
references: 19293140
references: J Anim Breed Genet. 2012 Dec;129(6):474-87
references: 23148973
references: Am J Hum Genet. 2004 Jun;74(6):1088-101
references: 15114533
references: Proc Natl Acad Sci U S A. 2000 Oct 24;97(22):12182-6
references: 11027309
references: Genetics. 2007 Nov;177(3):1889-913
references: 17947425
references: J Dairy Sci. 2008 Nov;91(11):4414-23
references: 18946147
references: BMC Bioinformatics. 2009 Nov 24;10:384
references: 19930695
references: J Am Stat Assoc. 2009 Jun 1;104(486):735-746
references: 19881892
references: Genetics. 1988 Nov;120(3):849-52
references: 3224810
references: Genetics. 2010 Oct;186(2):713-24
references: 20813882