Population genetics of VNTR loci and their applications in evolutionary studies

Variable number of tandem repeats (VNTR) are genetic loci at which short sequence motifs are found repeated different numbers of times among chromosomes. To explore the potential utility of VNTR loci in evolutionary studies, I have conducted a series of studies to address the following questions: (1) What are the population genetic properties of these loci? (2) What are the mutational mechanisms of repeat number change at these loci? (3) Can DNA profiles be used to measure the relatedness between a pair of individuals? (4) Can DNA fingerprint be used to measure the relatedness between populations in evolutionary studies? (5) Can microsatellite and short tandem repeat (STR) loci which mutate stepwisely be used in evolutionary analyses?^ A large number of VNTR loci typed in many populations were studied by means of statistical methods developed recently. The results of this work indicate that there is no significant departure from Hardy-Weinberg expectation (HWE) at VNTR loci in most of the human populations examined, and the departure from HWE in some VNTR loci are not solely caused by the presence of population sub-structure.^ A statistical procedure is developed to investigate the mutational mechanisms of VNTR loci by studying the allele frequency distributions of these loci. Comparisons of frequency distribution data on several hundreds VNTR loci with the predictions of two mutation models demonstrated that there are differences among VNTR loci grouped by repeat unit sizes.^ By extending the ITO method, I derived the distribution of the number of shared bands between individuals with any kinship relationship. A maximum likelihood estimation procedure is proposed to estimate the relatedness between individuals from the observed number of shared bands between them.^ It was believed that classical measures of genetic distance are not applicable to analysis of DNA fingerprints which reveal many minisatellite loci simultaneously in the genome, because the information regarding underlying alleles and loci is not available. I proposed a new measure of genetic distance based on band sharing between individuals that is applicable to DNA fingerprint data.^ To address the concern that microsatellite and STR loci may not be useful for evolutionary studies because of the convergent nature of their mutation mechanisms, by a theoretical study as well as by computer simulation, I conclude that the possible bias caused by the convergent mutations can be corrected, and a novel measure of genetic distance that makes the correction is suggested. In summary, I conclude that hypervariable VNTR loci are useful in evolutionary studies of closely related populations or species, especially in the study of human evolution and the history of geographic dispersal of Homo sapiens. (Abstract shortened by UMI.) ^

Analysis of recurrent failure times: A time-dependent Yule process approach

Analysis of recurrent events has been widely discussed in medical, health services, insurance, and engineering areas in recent years. This research proposes to use a nonhomogeneous Yule process with the proportional intensity assumption to model the hazard function on recurrent events data and the associated risk factors. This method assumes that repeated events occur for each individual, with given covariates, according to a nonhomogeneous Yule process with intensity function λx(t) = λ 0(t) · exp( x′β). One of the advantages of using a non-homogeneous Yule process for recurrent events is that it assumes that the recurrent rate is proportional to the number of events that occur up to time t. Maximum likelihood estimation is used to provide estimates of the parameters in the model, and a generalized scoring iterative procedure is applied in numerical computation. ^ Model comparisons between the proposed method and other existing recurrent models are addressed by simulation. One example concerning recurrent myocardial infarction events compared between two distinct populations, Mexican-American and Non-Hispanic Whites in the Corpus Christi Heart Project is examined. ^

Outcome-based adaptive randomization for binary data in clinical trials---Compare and contrast

Monte Carlo simulation has been conducted to investigate parameter estimation and hypothesis testing in some well known adaptive randomization procedures. The four urn models studied are Randomized Play-the-Winner (RPW), Randomized Pôlya Urn (RPU), Birth and Death Urn with Immigration (BDUI), and Drop-the-Loses Urn (DL). Two sequential estimation methods, the sequential maximum likelihood estimation (SMLE) and the doubly adaptive biased coin design (DABC), are simulated at three optimal allocation targets that minimize the expected number of failures under the assumption of constant variance of simple difference (RSIHR), relative risk (ORR), and odds ratio (OOR) respectively. Log likelihood ratio test and three Wald-type tests (simple difference, log of relative risk, log of odds ratio) are compared in different adaptive procedures. ^ Simulation results indicates that although RPW is slightly better in assigning more patients to the superior treatment, the DL method is considerably less variable and the test statistics have better normality. When compared with SMLE, DABC has slightly higher overall response rate with lower variance, but has larger bias and variance in parameter estimation. Additionally, the test statistics in SMLE have better normality and lower type I error rate, and the power of hypothesis testing is more comparable with the equal randomization. Usually, RSIHR has the highest power among the 3 optimal allocation ratios. However, the ORR allocation has better power and lower type I error rate when the log of relative risk is the test statistics. The number of expected failures in ORR is smaller than RSIHR. It is also shown that the simple difference of response rates has the worst normality among all 4 test statistics. The power of hypothesis test is always inflated when simple difference is used. On the other hand, the normality of the log likelihood ratio test statistics is robust against the change of adaptive randomization procedures. ^

Logistic regression with Markov chains as covariates

The tobacco-specific nitrosamine 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) is an obvious carcinogen for lung cancer. Since CBMN (Cytokinesis-blocked micronucleus) has been found to be extremely sensitive to NNK-induced genetic damage, it is a potential important factor to predict the lung cancer risk. However, the association between lung cancer and NNK-induced genetic damage measured by CBMN assay has not been rigorously examined. ^ This research develops a methodology to model the chromosomal changes under NNK-induced genetic damage in a logistic regression framework in order to predict the occurrence of lung cancer. Since these chromosomal changes were usually not observed very long due to laboratory cost and time, a resampling technique was applied to generate the Markov chain of the normal and the damaged cell for each individual. A joint likelihood between the resampled Markov chains and the logistic regression model including transition probabilities of this chain as covariates was established. The Maximum likelihood estimation was applied to carry on the statistical test for comparison. The ability of this approach to increase discriminating power to predict lung cancer was compared to a baseline "non-genetic" model. ^ Our method offered an option to understand the association between the dynamic cell information and lung cancer. Our study indicated the extent of DNA damage/non-damage using the CBMN assay provides critical information that impacts public health studies of lung cancer risk. This novel statistical method could simultaneously estimate the process of DNA damage/non-damage and its relationship with lung cancer for each individual.^

Foraminifer counts in core JM01-1199, and review of radiocarbon ages in the Norwegian-Barents Sea-Svalbard margin

We compare six high-resolution Holocene, sediment cores along a S-N transect on the Norwegian-Svalbard continental margin from ca 60°N to 77.4°N, northern North Atlantic. Planktonic foraminifera in the cores were investigated to show the changes in upper surface and subsurface water mass distribution and properties, including summer sea-surface temperatures (SST). The cores are located below the axis of the Norwegian Current and the West Spitsbergen Current, which today transport warm Atlantic Water to the Arctic. Sediment accumulation rates are generally high at all the core sites, allowing for a temporal resolution of 10-102 years. SST is reconstructed using different types of transfer functions, resulting in very similar SST trends, with deviations of no more than +- 1.0/1.5 °C. A transfer function based on the maximum likelihood statistical approach is found to be most relevant. The reconstruction documents an abrupt change in planktonic foraminiferal faunal composition and an associated warming at the Younger Dryas-Preboreal transition. The earliest part of the Holocene was characterized by large temperature variability, including the Preboreal Oscillations and the 8.2 k event. In general, the early Holocene was characterized by SSTs similar to those of today in the south and warmer than today in the north, and a smaller S-N temperature gradient (0.23 °C/°N) compared to the present temperature gradient (0.46 °C/°N). The southern proxy records (60-69°N) were more strongly influenced by slightly cooler subsurface water probably due to the seasonality of the orbital forcing and increased stratification due to freshening. The northern records (72-77.4°N) display a millennial-scale change associated with reduced insolation and a gradual weakening of the North Atlantic thermohaline circulation (THC). The observed northwards amplification of the early Holocene warming is comparable to the pattern of recent global warming and future climate modelling, which predicts greater warming at higher latitudes. The overall trend during mid and late Holocene was a cooling in the north, stable or weak warming in the south, and a maximum S-N SST gradient of ca 0.7 °C/°N at 5000 cal. years BP. Superimposed on this trend were several abrupt temperature shifts. Four of these shifts, dated to 9000-8000, 5500-3000 and 1000 and ~400 cal. years BP, appear to be global, as they correlate with periods of global climate change. In general, there is a good correlation between the northern North Atlantic temperature records and climate records from Norway and Svalbard.

Detection histories for eight species of Amazona parrots in Venezuela during the NeoMaps bird surveys in 2010

Documenting changes in distribution is necessary for understanding species' response to environmental changes, but data on species distributions are heterogeneous in accuracy and resolution. Combining different data sources and methodological approaches can fill gaps in knowledge about the dynamic processes driving changes in species-rich, but data-poor regions. We combined recent bird survey data from the Neotropical Biodiversity Mapping Initiative (NeoMaps) with historical distribution records to estimate potential changes in the distribution of eight species of Amazon parrots in Venezuela. Using environmental covariates and presence-only data from museum collections and the literature, we first used maximum likelihood to fit a species distribution model (SDM) estimating a historical maximum probability of occurrence for each species. We then used recent, NeoMaps survey data to build single-season occupancy models (OM) with the same environmental covariates, as well as with time- and effort-dependent detectability, resulting in estimates of the current probability of occurrence. We finally calculated the disagreement between predictions as a matrix of probability of change in the state of occurrence. Our results suggested negative changes for the only restricted, threatened species, Amazona barbadensis, which has been independently confirmed with field studies. Two of the three remaining widespread species that were detected, Amazona amazonica, Amazona ochrocephala, also had a high probability of negative changes in northern Venezuela, but results were not conclusive for Amazona farinosa. The four remaining species were undetected in recent field surveys; three of these were most probably absent from the survey locations (Amazona autumnalis, Amazona mercenaria and Amazona festiva), while a fourth (Amazona dufresniana) requires more intensive targeted sampling to estimate its current status. Our approach is unique in taking full advantage of available, but limited data, and in detecting a high probability of change even for rare and patchily-distributed species. However, it is presently limited to species meeting the strong assumptions required for maximum-likelihood estimation with presence-only data, including very high detectability and representative sampling of its historical distribution.

Turbidity measurements in 96-wells microplates to determine the densities of copiotrophic and oligotrophic bacteria with the Most-Probable-Number-method during POLARSTERN cruise ANT-XI/2

A new microtiter-plate dilution method was applied during the expedition ANTARKTIS-XI/2 with RV Polarstern to determine the distribution of copiotrophic and oligotrophic bacteria in the water columns at polar fronts. Twofold serial dilutions were performed with an eight-channel Electrapette in 96-wells plates by mixing 150 µl of seawater with 150 µl of copiotrophic or olitrophic Trypticase-Broth, three times per well. After incubation of about 6 month at 2 °C, turbidities were measured with an eight-channel photometer at 405 nm and combinations of positive test results for three consecutive dilutions chosen and compared with a Most Probable Number table, calculated for 8 replicates and twofold serial dilutions. Densities of 12 to 661 cells/ml for copiotrophs, and 1 to 39 cells/ml for oligotrophs were found. Colony Forming Units on copiotrophic Trypticase-Agar were between 6 and 847 cells/ml, which is in the same range as determined with the MPN method.