Hardy-Weinberg Equilibrium and the Ternary Plot

The Hardy-Weinberg law, formulated about 100 years ago, states that under certain assumptions, the three genotypes AA, AB and BB at a bi-allelic locus are expected to occur in the proportions p2, 2pq, and q2 respectively, where p is the allele frequency of A, and q = 1-p. There are many statistical tests being used to check whether empirical marker data obeys the Hardy-Weinberg principle. Among these are the classical xi-square test (with or without continuity correction), the likelihood ratio test, Fisher's Exact test, and exact tests in combination with Monte Carlo and Markov Chain algorithms. Tests for Hardy-Weinberg equilibrium (HWE) are numerical in nature, requiring the computation of a test statistic and a p-value. There is however, ample space for the use of graphics in HWE tests, in particular for the ternary plot. Nowadays, many genetical studies are using genetical markers known as Single Nucleotide Polymorphisms (SNPs). SNP data comes in the form of counts, but from the counts one typically computes genotype frequencies and allele frequencies. These frequencies satisfy the unit-sum constraint, and their analysis therefore falls within the realm of compositional data analysis (Aitchison, 1986). SNPs are usually bi-allelic, which implies that the genotype frequencies can be adequately represented in a ternary plot. Compositions that are in exact HWE describe a parabola in the ternary plot. Compositions for which HWE cannot be rejected in a statistical test are typically “close" to the parabola, whereas compositions that differ significantly from HWE are “far". By rewriting the statistics used to test for HWE in terms of heterozygote frequencies, acceptance regions for HWE can be obtained that can be depicted in the ternary plot. This way, compositions can be tested for HWE purely on the basis of their position in the ternary plot (Graffelman & Morales, 2008). This leads to nice graphical representations where large numbers of SNPs can be tested for HWE in a single graph. Several examples of graphical tests for HWE (implemented in R software), will be shown, using SNP data from different human populations

Quantifying rock fabrics: a test of autocorrelation of the spatial distribution of cristals

A novel test of spatial independence of the distribution of crystals or phases in rocks based on compositional statistics is introduced. It improves and generalizes the common joins-count statistics known from map analysis in geographic information systems. Assigning phases independently to objects in RD is modelled by a single-trial multinomial random function Z(x), where the probabilities of phases add to one and are explicitly modelled as compositions in the K-part simplex SK. Thus, apparent inconsistencies of the tests based on the conventional joins{count statistics and their possibly contradictory interpretations are avoided. In practical applications we assume that the probabilities of phases do not depend on the location but are identical everywhere in the domain of de nition. Thus, the model involves the sum of r independent identical multinomial distributed 1-trial random variables which is an r-trial multinomial distributed random variable. The probabilities of the distribution of the r counts can be considered as a composition in the Q-part simplex SQ. They span the so called Hardy-Weinberg manifold H that is proved to be a K-1-affine subspace of SQ. This is a generalisation of the well-known Hardy-Weinberg law of genetics. If the assignment of phases accounts for some kind of spatial dependence, then the r-trial probabilities do not remain on H. This suggests the use of the Aitchison distance between observed probabilities to H to test dependence. Moreover, when there is a spatial uctuation of the multinomial probabilities, the observed r-trial probabilities move on H. This shift can be used as to check for these uctuations. A practical procedure and an algorithm to perform the test have been developed. Some cases applied to simulated and real data are presented. Key words: Spatial distribution of crystals in rocks, spatial distribution of phases, joins-count statistics, multinomial distribution, Hardy-Weinberg law, Hardy-Weinberg manifold, Aitchison geometry

Distribució espacial i ús de l'hàbitat de la comunitat de peixos a l'estany de Banyoles

En el marc d'un projecte més ampli sobre la comunitat de peixos de la conca lacustre de Banyoles, s'ha estudiat la distribució espacial de les diferents espècies així com l'ús que aquestes fan de l'hàbitat. El poblament piscícola de l'estany de Banyoles és el resultat d'un llarg historial d'introducció d'espècies exòtiques i extinció de les poblacions autòctones locals. S'ha revisat la seva composició actual detectant un total de 18 espècies (4 autòctones i 14 introduïdes) de les quals només 13 presenten una població estable. L'estudi de l'hàbitat s'ha centrat a l'Estany, l'element principal del sistema lacustre, analitzant per separat la zona litoral i la zona limnètica. En el primer cas s'han realitzat les captures d'individus mitjançant pesca elèctrica des d'una embarcació adaptada per aquesta pràctica. La totalitat del litoral ha estat dividida en trams de característiques homogènies on s'han obtingut les captures per unitat d'esforç per cada espècie. El mostreig s'ha desenvolupat entre l'estiu de 1997 i la primavera del 2000 realitzant un total de 10 campanyes de pesca. Les espècies més abundants al litoral són la perca americana (Micropterus salmoides) i el peix sol (Lepomis gibbosus), essent també presents la perca (Perca fluviatilis), carpa (Cyprinus carpio) i el gardí (Scardinius erythrophthalmus). S'han capturat altres espècies com ara anguila (Anguilla anguilla), bagra (Squalius cephalus), sandra (Sander lucioperca), carpí (Carassius auratus) i madrilleta vera (Rutilus rutilus), però són molt menys abundants en nombre. S'ha examinat, per cadascuna de les espècies, si existeix selecció de l'hàbitat i en cas afirmatiu, quin és el preferent en base a la classificació del litoral en sis tipus de vegetació predominant. Les espècies més abundants, perca americana i peix sol, ocupen tots els hàbitats disponibles però amb una densitat diferent. La perca mostra també una clara selecció de l'hàbitat a favor de les zones molt estructurades amb abundant presència de jonca litoral. Carpa i gardí seleccionen els ambients més fondos amb major presència de matèria orgànica d'origen vegetal procedent del bogar. En general els individus ocupen les zones amb una densitat de vegetació intermèdia, majoritàriament zones de jonca a l'estiu i zones amb mansega a l'hivern, on troben refugi i els recursos tròfics necessaris. La perca americana, a més, presenta una elevada fidelitat a un mateix punt del litoral al llarg de la seva vida. La zona limnètica ha estat prospectada mensualment mitjançant ecosondació, realitzant transectes perpendiculars a l'eix principal de l'Estany, cobrint la seva totalitat. La composició d'espècies s'ha obtingut a partir de les captures fetes amb xarxes (tresmalls) amb periodicitat estacional. L'anàlisi geoestadística de la densitat de peixos ha permès descriure l'estructura espacial d'aquesta a partir dels variogrames, així com la seva variabilitat tant espacial com temporal, i obtenir els mapes de densitat. A l'hivern, la densitat de peixos a la zona limnètica assoleix els seus valors mínims i els individus es troben formant agregats dispersos, pels diferents estrats de fondària. A partir de la primavera la densitat augmenta, pel reclutament i la major freqüència d'individus que abandonen el litoral; la densitat esdevé més homogènia a les primer capes de fondària. A l'estiu la densitat és màxima i l'estrat més homogeni coincideix amb la posició de la termoclina. Aquest estructura varia en disminuir la temperatura i barrejar-se la columna d'aigua, tornant a la situació hivernal. La perca i la madrilleta vera són les espècies predominants en aquest ambient, juntament amb la carpa. La seva distribució no és homogènia i respon a les característiques limnològiques de les diferents cubetes de l'Estany. Una particularitat d'aquest, relacionada amb el seu origen càrstic, es la formació d'una ploma hidrotermal que afecta la distribució dels peixos, probablement en augmentar la terbolesa. S'ha integrat l'ús de l'hàbitat de les espècies que ocupen tant la zona limnètica com la litoral a partir del seguiment d'individus, concretament de perca i bagra. S'ha utilitzat un sistema automàtic de posicionament que estima la localització dels individus marcats amb transmissor de telemetria acústica. L'anàlisi dels desplaçaments mostra un rang superior per la bagra en comparació amb la perca. Ambdues espècies mostren una orientació en els seus desplaçaments. La perca ocupa el litoral a la nit i es desplaça a la zona limnètica de dia, amb un ritme d'activitat marcat per dos màxims coincidint amb la sortida i posta de sol; en canvi la bagra mostra una major activitat nocturna amb zones de repòs properes al litoral. S'ha estimat igualment els dominis vital de cada individu marcat.