Many advances in the understanding of meiosis have been made by measuring how often errors in chromosome segregation occur. This process of nondisjunction can be studied by counting experimental progeny, but direct measurement of nondisjunction rates is complicated by not all classes of nondisjunctional progeny being viable. For
X chromosome nondisjunction in Drosophila female meiosis, all of the normal progeny survive, while nondisjunctional eggs produce viable progeny only if fertilized by sperm that carry the appropriate sex chromosome. The rate of nondisjunction has traditionally been estimated by assuming a binomial process and doubling the number of observed nondisjunctional progeny, to account for the inviable classes. However, the correct way to derive statistics (such as confidence intervals or hypothesis testing) by this approach is far from clear. Instead, we use the multinomial-Poisson hierarchy model and demonstrate that the old estimator is in fact the maximum-likelihood estimator (MLE). Under more general assumptions, we derive asymptotic normality of this estimator and construct confidence interval and hypothesis testing formulae. Confidence intervals under this framework are always larger than under the binomial framework, and application to published data shows that use of the multinomial approach can avoid an apparent type 1 error made by use of the binomial assumption. The current study provides guidance for researchers designing genetic experiments on nondisjunction and improves several methods for the analysis of genetic data.MEIOSIS is a specialized cell division, where a diploid cell undergoes a single round of replication followed by two rounds of segregation to produce four haploid gametes. During this segregation, chromosomes must correctly separate (or disjoin) from their homologs at meiosis I, followed by sister chromatids disjoining at meiosis II. When chromosomes fail to disjoin from their partners, the resultant nondisjunction produces aneuploid gametes with the wrong number of chromosomes. The study of meiotic nondisjunction in Drosophila has a long and distinguished history of publication in
genetics, with the inaugural article published in this journal being Calvin Bridges'' use of nondisjunction to prove the chromosome theory of heredity (
Bridges 1916). The first study that screened variants isolated from natural populations used nondisjunction to identify meiotic mutants (
Sandler et al. 1968), as did the first EMS-induced mutant screen (
Baker and Carpenter 1972). Subsequent screens using new mutagens or techniques have also relied on measuring nondisjunction to identify mutants of interest (
Sekelsky et al. 1999). Indeed, much of the progress that has been made in the study of meiosis would not have been possible without the use of nondisjunction to identify new mutations that are defective at some step in chromosome segregation.However, one difficulty in estimating nondisjunction rates is that in most instances the resulting aneuploid progeny cannot survive. Fortunately, in Drosophila it is possible to design crosses to recover them. Sex determination in flies is based on the number of
X chromosomes, rather than a masculinizing
Y chromosome as in mammals. This means that
XO flies are viable (but sterile) males, while
XXY flies are viable females. Therefore, it is possible to recover both normal and nondisjunctional progeny, as a nullo-
X egg fertilized by an
X-bearing sperm will survive as an
XO male, while a diplo-
X egg fertilized by a sperm lacking an
X will be female (
XXY). By using visible markers on the sex chromosomes, these exceptional progeny are straightforward to identify. However, if those eggs are fertilized by the other class of sperm, the resulting
OY or
XXX progeny are inviable. Therefore, the nondisjunction rate that occurs during meiosis is not equal to the proportion of nondisjunctional progeny, as only 50% of nondisjunctional eggs receive sperm compatible with viability, while all normal eggs are viable.Given this experimental limitation, what is the correct method to calculate the error rate during meiosis? For this discussion, let
N be the total number of progeny produced in an experiment, let
X1 be the number of inviable nondisjunctional progeny (
OY and
XXX), let
X2 be the number of viable nondisjunctional progeny (
XO and
XXY), and let
X3 be the number of normal progeny (
XY and
XX), such that
N =
X1 +
X2 +
X3. If all progeny could be counted, then the nondisjunction rate would simply be (
X1 +
X2)/
N.However, only flies that survive to adulthood can be counted, and therefore both
X1 and
N are unknown. As
X- and
Y-bearing sperm are produced in equal numbers, live and dead nondisjunctional progeny are also expected in equal numbers. Therefore, K.W. Cooper (
Cooper 1948) proposed the widely used estimator for the
X chromosome nondisjunction rate, where
X2 is substituted for
X1 in the above formula, giving the rate as:(1)While this estimator works, the statistical properties of this estimator are not clear. Instead of following the early literature to combine
X1 and
X2 and use a binomial distribution, we go back to the three original categories and model the process as a multinomial distribution with latent number of progeny
N, considering all three possible phenotypes for each progeny (nondisjunctional dead, nondisjunctional living, and normal). Whether a nondisjunctional oocyte becomes a nondisjunctional dead or nondisjunctional living progeny depends on the sex chromosome content of the sperm that fertilized it. As
X- and
Y-bearing sperm are produced in equal numbers during male meiosis, the usual genetic expectation for the rates of nondisjunctional dead and living progeny will be . However, even assuming that the rates of nondisjunctional dead and living progeny are different, with a Poisson assumption of
N, we can derive the maximum-likelihood estimators (MLEs) for the nondisjunctional dead and nondisjunctional living rates. Under the usual genetic expectation of equality, the MLE of the nondisjunctional rate coincides with Cooper''s estimator, and we furthermore derive the exact distribution of . Under another set of reasonable assumptions, we show the consistency and asymptotic normality of Cooper''s estimator, and derive asymptotic results when comparing two nondisjunction rates. All these distributional results enable us to develop confidence interval and hypothesis testing related to
p, or
px −
py in the case of comparing two nondisjunction rates from populations
x and
y.
相似文献