Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

Genetic constraints on the evolution of mate recognition under natural selection

Field populations of Drosophila serrata display reproductive character displacement in cuticular hydrocarbons (CHCs) when sympatric with Drosophila birchii. We have previously shown that the naturally occurring pattern of reproductive character displacement can be experimentally replicated by exposing field allopatric populations of D. serrata to experimental sympatry with D. birchii. Here, we tested whether the repeated evolution of reproductive character displacement in natural and experimental populations was a consequence of genetic constraints on the evolution of CHCs. The genetic variance-covariance (G) matrices for CHCs were determined for populations of D. serrata that had evolved in either the presence or absence of D. birchii under field and experimental conditions. Natural selection on mate recognition under both field and experimental sympatric conditions increased the genetic variance in CHCs consistent with a response to selection based on rare alleles. A close association between G eigenstructure and the eigenstructure of the phenotypic divergence (D) matrix in natural and experimental populations suggested that G matrix eigenstructure may have determined the direction in which reproductive character displacement evolved during the reinforcement of mate recognition.

Identification of the enamelin (g.8344delG) mutation in a new kindred and presentation of a standardized ENAM nomenclature

The amelogenesis imperfectas (Al) area geneticatly heterogeneous group of diseases that result in defective development of tooth enamel. Although X-linked, autosomal. dominant and autosomal. recessive forms of Al have been clinically characterized, only two genes (AMELX and ENAM) have been associated with Al. To date, three enamelin (ENAM) mutations have been identified. These mutations cause phenotypically diverse forms of autosomal. dominant Al. Detailed phenotype-genotype correlations have not been performed for autosomal. dominant Al due to ENAM mutations. We identified a previously unreported kindred segregating for the ENAM mutation, g.8344delG. Light and electron microscopy analyses of unerupted permanent teeth show the enamel is markedly reduced in thickness, Lacks a prismatic structure and has a laminated appearance. Taken together these histological features support the enamelin protein as being critical for the development of a normal. enamel. thickness and that it Likely has a role in regulating c-axis crystallite growth. Because there is growing molecular and phenotypic diversity in the enamelin defects, it is critical to have a nomenclature and numbering system for characterizing these conditions. We present a standardized nomenclature for ENAM mutations that will allow consistent reporting and communication. (C) 2003 Elsevier Science Ltd. All rights reserved.

Quantitative trait loci for agronomic and seed quality traits in an F-2 and F-4 : 6 soybean population

Molecular breeding is becoming more practical as better technology emerges. The use of molecular markers in plant breeding for indirect selection of important traits can favorably impact breeding efficiency. The purpose of this research is to identify quantitative trait loci (QTL) on molecular linkage groups (MLG) which are associated with seed protein concentration, seed oil concentration, seed size, plant height, lodging, and maturity, in a population from a cross between the soybean cultivars 'Essex' and 'Williams.' DNA was extracted from F-2 generation soybean leaves and amplified via polymerase chain reaction (PCR) using simple sequence repeat (SSR) markers. Markers that were polymorphic between the parents were analyzed against phenotypic trait data from the F-2 and F-4:6 generation. For the F-2 population, significant additive QTL were Satt540 (MLG M, maturity, r(2)=0.11; height, r(2)=0.04, seed size, r(2)=0.061, Satt373 (MLG L, seed size, r(2)=0.04; height, r(2)=0.14), Satt50 (MLG A1, maturity r(2)=0.07), Satt14 (MLG D2, oil, r(2)=0.05), and Satt251 (protein r(2)=0.03, oil, r(2)=0.04). Significant dominant QTL for the F-2 population were Satt540 (MLG M, height, r(2)=0.04; seed size, r(2)=0.06) and Satt14 (MLG D2, oil, r(2)=0.05). In the F-4:6 generation significant additive QTL were Satt239 (MLG I, height, r(2)=0.02 at Knoxville, TN and r(2)=0.03 at Springfield, TN), Satt14 (MLG D2, seed size, r(2)=0.14 at Knoxville, TN), Satt373 (MLG L, protein, r(2)=0.04 at Knoxville, TN) and Satt251 (MLG B I, lodging r(2)=0.04 at Springfield, TN). Averaged over both environments in the F-4:6 generation, significant additive QTL were identified as Satt251 (MLG B 1, protein, r(2)=0.03), and Satt239 (MLG 1, height, r(2)=0.03). The results found in this study indicate that selections based solely on these QTL would produce limited gains (based on low r(2) values). Few QTL were detected to be stable across environments. Further research to identify stable QTL over environments is needed to make marker-assisted approaches more widely adopted by soybean breeders.

Identification of residues and domains of Raf important for function in vivo and in vitro

Random mutagenesis and genetic screens for impaired Raf function in Caenorhabditis elegans were used to identify six loss-of-function alleles of lin-45 raf that result in a substitution of a single amino acid. The mutations were classified as weak, intermediate, and strong based on phenotypic severity. We engineered these mutations into the homologous residues of vertebrate Raf-1 and analyzed the mutant proteins for their underlying biochemical defects. Surprisingly, phenotype strength did not correlate with the catalytic activity of the mutant proteins. Amino acid substitutions Val-589 and Ser-619 severely compromised Raf kinase activity, yet these mutants displayed weak phenotypes in the genetic screen. Interestingly, this is because these mutant Raf proteins efficiently activate the MAPK (mitogen-activated protein kinase) cascade in living cells, a result that may inform the analysis of knockout mice. Equally intriguing was the observation that mutant proteins with non-functional Ras-binding domains, and thereby deficient in Ras-mediated membrane recruitment, displayed only intermediate strength phenotypes. This confirms that secondary mechanisms exist to couple Ras to Raf in vivo. The strongest phenotype in the genetic screens was displayed by a S508N mutation that again did not correlate with a significant loss of kinase activity or membrane recruitment by oncogenic Ras in biochemical assays. Ser-508 lies within the Raf-1 activation loop, and mutation of this residue in Raf-1 and the equivalent Ser-615 in B-Raf revealed that this residue regulates Raf binding to MEK. Further characterization revealed that in response to activation by epidermal growth factor, the Raf-S508N mutant protein displayed both reduced catalytic activity and aberrant activation kinetics: characteristics that may explain the C. elegans phenotype.

The role of melanocortin-1 receptor polymorphism in skin cancer risk phenotypes

We have examined melanocortin-1 receptor (MC1R) variant allele frequencies in the general population and in a collection of adolescent dizygotic and monozygotic twins to determine statistical associations of pigmentation phenotypes with increased skin cancer risk. This included hair and skin color, freckling, mole count and sun exposed skin reflectance. Nine variants were studied and designated as either strong R (OR = 63; 95% CI 32-140) or weak r (OR = 5; 95% CI 3-11) red hair alleles. Penetrance of each MC1R variant allele was consistent with an allelic model where effects were multiplicative for red hair but additive for skin reflectance. To assess the interaction of the brown eye color gene BEY2/OCA2 on the phenotypic effects of variant MC1R alleles we imputed OCA2 genotype in the twin collection. A modifying effect of OCA2 on MC1R variant alleles was seen on constitutive skin color, freckling and mole count. In order to study the individual effects of these variants on pigmentation phenotype we have established a series of human primary melanocyte strains genotyped for the MC1R receptor. These include strains which are MC1R wild-type consensus, variant heterozygotes, and homozygotes for strong R alleles Arg151Cys and Arg160Trp. Ultrastructural analysis demonstrated that only consensus strains contained stage III and IV melanosomes in their terminal dendrites whereas Arg151Cys and Arg160Trp homozygous strains contained only immature stage I and II melanosomes. Such genetic association studies combined with the functional analysis of MC1R variant alleles in melanocytic cells should provide a link in understanding the association between pigmentary phototypes and skin cancer risk.

The t-SNARE syntaxin 4 is regulated during macrophage activation to function in membrane traffic and cytokine secretion

Activation of macrophages with lipopolysaccharide (LPS) induces the rapid synthesis and secretion of proinflammatory cytokines, such as tumor necrosis factor (TNFalpha), for priming the immune response [1, 2]. TNFalpha plays a key role in inflammatory disease [3]; yet, little is known of the intracellular trafficking events leading to its secretion. In order to identify molecules involved in this secretory pathway, we asked whether any of the known trafficking proteins are regulated by LPS. We found that the levels of SNARE proteins were rapidly and significantly up- or downregulated during macrophage activation. A subset of t-SNAREs (Syntaxin 4/SNAP23/Munc18c) known to control regulated exocytosis in other cell types [4, 5] was substantially increased by LPS in a temporal pattern coinciding with peak TNFalpha secretion. Syntaxin 4 formed a complex with Munc18c at the cell surface of macrophages. Functional studies involving the introduction of Syntaxin 4 cDNA or peptides into macrophages implicate this t-SNARE in a rate-limiting step of TNFalpha secretion and in membrane ruffling during macrophage activation. We conclude that in macrophages, SNAREs are regulated in order to accommodate the rapid onset of cytokine secretion and for membrane traffic associated with the phenotypic changes of immune activation. This represents a novel regulatory role for SNAREs in regulated secretion and in macrophage-mediated host defense.

The use of 16s rDNA restriction fragment length polymorphism analysis for the identification of non-fermenting gram negative bacilli in a cystic fibrosis microbiology referral service

The utility of 16s rDNA restriction fragment length polymorphism (RFLP) analysis for the partial genomovar differentiation of Burkholderia cepacia complex bacterium is well documented. We compared the 16s rDNA RFLP signatures for a number of non-fermenting gram negative bacilli (NF GNB) LMG control strains and clinical isolates pertaining to the genera Burkholderia, Pseudomonas, Achromobacter (Alcaligenes), Ralstonia, Stenotrophomonas and Pandoraea. A collection of 24 control strain (LMG) and 25 clinical isolates were included in the study. Using conventional PCR, a 1.2 kbp 16s rDNA fragment was generated for each organism. Following restriction digestion and electrophoresis, each clinical isolate RFLP signature was compared to those of the control strain panel. Nineteen different RFLP signatures were detected from the 28 control strains included in the study. TwentyoneyTwenty- five of the clinical isolates could be classified by RFLP analysis into a single genus and species when compared to the patterns produced by the control strain panel. Four clinical B. pseudomallei isolates produced RFLP signatures which were indistinguishable from B. cepacia genomovars I, III and VIII. The identity of these four isolates were confirmed using B. pseudomallei specific PCR. 16s rDNA RFLP analysis can be a useful identification strategy when applied to NF GNB, particularly for those which exhibit colistin sulfate resistance. The use of this molecular based methodology has proved very useful in the setting of a CF referral laboratory particularly when utilised in conjunction with B. cepacia complex and genomovar specific PCR techniques. Species specific PCR or sequence analysis should be considered for selected isolates; especially where discrepancies between epidemiology, phenotypic and genotypic characteristics occur.

The pyramidal neuron in occipital, temporal and prefrontal cortex of the owl monkey (Aotus trivirgatus): regional specialization in cell structure

Recent studies have revealed marked regional variation in pyramidal cell morphology in primate cortex. In particular, pyramidal cells in human and macaque prefrontal cortex (PFC) are considerably more spinous than those in other cortical regions. PFC pyramidal cells in the New World marmoset monkey, however, are less spinous than those in man and macaques. Taken together, these data suggest that the pyramidal cell has become more branched and more spinous during the evolution of PFC in only some primate lineages. This specialization may be of fundamental importance in determining the cognitive styles of the different species. However, these data are preliminary, with only one New World and two Old World species having been studied. Moreover, the marmoset data were obtained from different cases. In the present study we investigated PFC pyramidal cells in another New World monkey, the owl monkey, to extend the basis for comparison. As in the New World marmoset monkey, prefrontal pyramidal cells in owl monkeys have relatively few spines. These species differences appear to reflect variation in the extent to which PFC circuitry has become specialized during evolution. Highly complex pyramidal cells in PFC appear not to have been a feature of a common prosimian ancestor, but have evolved with the dramatic expansion of PFC in some anthropoid lineages.