Evaluation of computational methods for the reconstruction of HLA haplotypes

Human leukocyte antigen (HLA) haplotypes are frequently evaluated for population history inferences and association studies. However, the available typing techniques for the main HLA loci usually do not allow the determination of the allele phase and the constitution of a haplotype, which may be obtained by a very time-consuming and expensive family-based segregation study. Without the family-based study, computational inference by probabilistic models is necessary to obtain haplotypes. Several authors have used the expectation-maximization (EM) algorithm to determine HLA haplotypes, but high levels of erroneous inferences are expected because of the genetic distance among the main HLA loci and the presence of several recombination hotspots. In order to evaluate the efficiency of computational inference methods, 763 unrelated individuals stratified into three different datasets had their haplotypes manually defined in a family-based study of HLA-A, -B, -DRB1 and -DQB1 segregation, and these haplotypes were compared with the data obtained by the following three methods: the Expectation-Maximization (EM) and Excoffier-Laval-Balding (ELB) algorithms using the arlequin 3.11 software, and the PHASE method. When comparing the methods, we observed that all algorithms showed a poor performance for haplotype reconstruction with distant loci, estimating incorrect haplotypes for 38%-57% of the samples considering all algorithms and datasets. We suggest that computational haplotype inferences involving low-resolution HLA-A, HLA-B, HLA-DRB1 and HLA-DQB1 haplotypes should be considered with caution.

Facies architecture of the CRP-3 drillhole, Victoria Land Basin, Antarctica

The Cenozoic Victoria Land Basin (VLB) stratigraphic section penetrated by CRP-3 is mostly of Early Oligocene age. It contains an array of lithofacies comprising fine-grained mudrocks, interlaminated and interbedded mudrocks/sandstones, mud-rich and mud-poor sandstones, conglomerates and diametites that are together interpreted as the products of shallow marine to possibly non-marine environments of deposition, affected by the periodic advance and retreat of tidewater glaciers. This lithofacies assemblage can be readily rationalised using the facies scheme designed originally for CRP-2/2A, and published previously. The uppermost 330 metres below sea floor (mbsf) shows a cyclical arrangement of lithofacies also similar to that recognised throughout CRP-2/2A, and interpreted to reflect cyclical variations in relative sea-level driven by ice volume fluctuations ("Motif A"). Between 330 and 480 mbsf, a series of less clearly cyclical units, generally fining-upward but nonetheless incorporating a significant subset of the facies assemblage, has been identified and noted in the Initial Report as "Motif B. Below 480 mbsf, the section is arranged into a repetitive succession of fining-upward units, each of which comprises dolerite clast conglomerate at the base passing upward into relatively thick intervals of sandstones. The cycles present down 480 mbsf are defined as sequences, each interpreted to record cyclical variation of relative sea-level. The thickness distribution of sequences in CRP-3 provides some insights into the geological variables controlling sediment accumulation in the Early Oligocene section. The uppermost part of the section in CRP-3 comprises two or three thick, complete sequences that show a broadly symmetrical arrangement of lithofacies (similar to Sequences 9-11 in CRP-2/2A). This suggests a period of relatively rapid tectonic subsidence, which allowed preservation of the complete facies cycle. Below Sequence 3, however, is a considerable interval of thin, incomplete and erosionally truncated sequences (4-23), which incorporates both the remainder of Motif A sequences and all Motif B sequences recognised. The thinner and more truncated sequences suggest sediment accumulation under conditions of reduced accommodation, and given the lack of evidence for glacial conditions (see Powell et al., this volume) tends to argue for a period of reduced tectonic subsidence. The section below 480 mbsf consists of a series of fining-upward, conglomerate to sandstone intervals which cannot be readily interpreted in terms of relative sea-level change. A relatively mudrock-rich interval above the basal conglomerate/breccia (782-762 mbsf) may record initial flooding of the basin during early rift subsidence. The lithostratigraphy summarised above has been linked to seismic reflection data using depth conversion techniques (Henrys et al., this volume). The three uppermost reflectors ("o", "p" and "q") correlate to the package of thick sequences 1-3, and several deeper reflectors can also be correlated to sequence boundaries. The package of thick Sequences 1-3 shows a sheet-like cross-sectional geometry on seismic reflection lines, unlike the similar package recognised in CRP-2/2A.

Combining the analyses of introgressive hybridisation and linkage mapping to investigate the genetic architecture of population divergence in the lake whitefish (Coregonus clupeaformis, Mitchill)

Adaptation and reproductive isolation, the engines of biological diversity, are still elusive when discussing the genetic bases of speciation. Namely, the number of genes and magnitude of selection acting positively or negatively on genomic traits implicated in speciation is contentious. Here, we describe the first steps of an ongoing research program aimed at understanding the genetic bases of population divergence and reproductive isolation in the lake whitefish (Coregonus clupeaformis). A preliminary linkage map originating from a hybrid cross between dwarf and normal ecotypes is presented, whereby some of the segregating AFLP markers were found to be conserved among natural populations. Maximum-likelihood was used to estimate hybrid indices from non-diagnostic markers at 998 AFLP loci. This allowed identification of the most likely candidate loci that have been under the influence of selection during the natural hybridisation of whitefish originating from different glacial races. As some of these loci could be identified on the linkage map, the possibility that selection of traits in natural populations may eventually be correlated to specific chromosomal regions was demonstrated. The future prospects and potential of these approaches to elucidate the genetic bases of adaptation and reproductive isolation among sympatric ecotypes of lake whitefish is discussed.

Branching processes and computational collapse of discretized unimodal mappings

In computer simulations of smooth dynamical systems, the original phase space is replaced by machine arithmetic, which is a finite set. The resulting spatially discretized dynamical systems do not inherit all functional properties of the original systems, such as surjectivity and existence of absolutely continuous invariant measures. This can lead to computational collapse to fixed points or short cycles. The paper studies loss of such properties in spatial discretizations of dynamical systems induced by unimodal mappings of the unit interval. The problem reduces to studying set-valued negative semitrajectories of the discretized system. As the grid is refined, the asymptotic behavior of the cardinality structure of the semitrajectories follows probabilistic laws corresponding to a branching process. The transition probabilities of this process are explicitly calculated. These results are illustrated by the example of the discretized logistic mapping.

Simulation of insect movement with respect to plant architecture and morphogenesis

Developments in computer and three dimensional (3D) digitiser technologies have made it possible to keep track of the broad range of data required to simulate an insect moving around or over the highly heterogeneous habitat of a plant's surface. Properties of plant parts vary within a complex canopy architecture, and insect damage can induce further changes that affect an animal's movements, development and likelihood of survival. Models of plant architectural development based on Lindenmayer systems (L-systems) serve as dynamic platforms for simulation of insect movement, providing ail explicit model of the developing 3D structure of a plant as well as allowing physiological processes associated with plant growth and responses to damage to be described and Simulated. Simple examples of the use of the L-system formalism to model insect movement, operating Lit different spatial scales-from insects foraging on an individual plant to insects flying around plants in a field-are presented. Such models can be used to explore questions about the consequences of changes in environmental architecture and configuration on host finding, exploitation and its population consequences. In effect this model is a 'virtual ecosystem' laboratory to address local as well as landscape-level questions pertinent to plant-insect interactions, taking plant architecture into account. (C) 2002 Elsevier Science B.V. All rights reserved.

In the vernacular: On the architecture of the National Museum of Australia

This Article does not have an abstract.