Power and sample size calculations for case-control genetic association tests when errors are present: application to single nucleotide polymorphisms

The purpose of this work is to quantify the effects that errors in genotyping have on power and the sample size necessary to maintain constant asymptotic Type I and Type II error rates (SSN) for case-control genetic association studies between a disease phenotype and a di-allelic marker locus, for example a single nucleotide polymorphism (SNP) locus. We consider the effects of three published models of genotyping errors on the chi-square test for independence in the 2 x 3 table. After specifying genotype frequencies for the marker locus conditional on disease status and error model in both a genetic model-based and a genetic model-free framework, we compute the asymptotic power to detect association through specification of the test's non-centrality parameter. This parameter determines the functional dependence of SSN on the genotyping error rates. Additionally, we study the dependence of SSN on linkage disequilibrium (LD), marker allele frequencies, and genotyping error rates for a dominant disease model. Increased genotyping error rate requires a larger SSN. Every 1% increase in sum of genotyping error rates requires that both case and control SSN be increased by 2-8%, with the extent of increase dependent upon the error model. For the dominant disease model, SSN is a nonlinear function of LD and genotyping error rate, with greater SSN for lower LD and higher genotyping error rate. The combination of lower LD and higher genotyping error rates requires a larger SSN than the sum of the SSN for the lower LD and for the higher genotyping error rate.     

Gene for selenium-dependent glutathione peroxidase maps to human chromosomes 3, 21 and X

A human glutathione peroxidase cDNA has been used as a probe to hybridize to DNAs isolated from human - rodent somatic cell hybrids that have segregated human chromosomes. A 609 bp probe which contains the entire coding region hybridizes to human chromosomes 3, 21 and Xp. Fragments of the cDNA coding sequence and of the 3' untranslated region were also used as probes. These fragments hybridized to each of the three chromosomes with the same efficiency, suggesting similarity between the loci, whereas an intronic probe detected only the gene on chromosome 3. The general organization of each gene was determined from the hybridization data. The data suggest that the locus on chromosome 3 is a functional gene containing a single intron and a pattern of restriction sites identical to those found in the cDNA coding sequence. The data also suggest that the sequences on chromosomes X and 21 have equal conservation of the 3' untranslated and coding sequences but do not contain introns, providing evidence that the latter two sequences are processed pseudogenes. A simple two allele polymorphism in PvuII digests was detected at the locus on chromosome 21.     

Comparison of three methods to estimate genetic ancestry and control for stratification in genetic association studies among admixed populations

Population stratification may confound the results of genetic association studies among unrelated individuals from admixed populations. Several methods have been proposed to estimate the ancestral information in admixed populations and used to adjust the population stratification in genetic association tests. We evaluate the performances of three different methods: maximum likelihood estimation, ADMIXMAP and Structure through various simulated data sets and real data from Latino subjects participating in a genetic study of asthma. All three methods provide similar information on the accuracy of ancestral estimates and control type I error rate at an approximately similar rate. The most important factor in determining accuracy of the ancestry estimate and in minimizing type I error rate is the number of markers used to estimate ancestry. We demonstrate that approximately 100 ancestry informative markers (AIMs) are required to obtain estimates of ancestry that correlate with correlation coefficients more than 0.9 with the true individual ancestral proportions. In addition, after accounting for the ancestry information in association tests, the excess of type I error rate is controlled at the 5% level when 100 markers are used to estimate ancestry. However, since the effect of admixture on the type I error rate worsens with sample size, the accuracy of ancestry estimates also needs to increase to make the appropriate correction. Using data from the Latino subjects, we also apply these methods to an association study between body mass index and 44 AIMs. These simulations are meant to provide some practical guidelines for investigators conducting association studies in admixed populations.     

Trend tests for case-control studies of genetic markers: power,sample size and robustness

The Cochran-Armitage trend test is commonly used as a genotype-based test for candidate gene association. Corresponding to each underlying genetic model there is a particular set of scores assigned to the genotypes that maximizes its power. When the variance of the test statistic is known, the formulas for approximate power and associated sample size are readily obtained. In practice, however, the variance of the test statistic needs to be estimated. We present formulas for the required sample size to achieve a prespecified power that account for the need to estimate the variance of the test statistic. When the underlying genetic model is unknown one can incur a substantial loss of power when a test suitable for one mode of inheritance is used where another mode is the true one. Thus, tests having good power properties relative to the optimal tests for each model are useful. These tests are called efficiency robust and we study two of them: the maximin efficiency robust test is a linear combination of the standardized optimal tests that has high efficiency and the MAX test, the maximum of the standardized optimal tests. Simulation results of the robustness of these two tests indicate that the more computationally involved MAX test is preferable.     

Structural organization of the human genome: distribution of nucleotides, Alu-repeats and exons in chromosomes 21 and 22]

Analysis of DNA sequences of the human chromosomes 21 and 22 performed using a specially designed MegaGene software allowed us to obtain the following results. Purine and pyrimidine nucleotide residues are unevenly distributed along both chromosomes, displaying maxima and minima (Y waves phi) with a period of about 3 Mbp. Distribution of G + C along both chromosomes has no distinct maxima and minima, however, chromosome 21 contains considerably less G + C than chromosome 22. Both exons and Alu repeats are unevenly distributed along chromosome 21: they are scarce in its left part and abundant in the right part, while MIR elements are quite monotonously spread along this chromosome. The Alu repeats show a wave-like distribution pattern similar for both repeat orientations. The number of the Alu repeats of opposite orientations was equal for both studied chromosomes, and this may be considered a new property of the human genome. The positive correlation between the exon and Alu distribution patterns along the chromosome, the concurrent distribution of Alu repeats in both orientations along the chromosome, and the equal copy numbers for Alu in direct and inverted orientations within an individual chromosome point to their important role in the human genome, and do not fit the notion that Alu repeats belong to parasitic (junk) DNA.     

Algorithmic co-optimization of genetic constructs and growth conditions: application to 6-ACA,a potential nylon-6 precursor

Optimizing bio-production involves strain and process improvements performed as discrete steps. However, environment impacts genotype and a strain that is optimal under one set of conditions may not be under different conditions. We present a methodology to simultaneously vary genetic and process factors, so that both can be guided by design of experiments (DOE). Advances in DNA assembly and gene insulation facilitate this approach by accelerating multi-gene pathway construction and the statistical interpretation of screening data. This is applied to a 6-aminocaproic acid (6-ACA) pathway in Escherichia coli consisting of six heterologous enzymes. A 32-member fraction factorial library is designed that simultaneously perturbs expression and media composition. This is compared to a 64-member full factorial library just varying expression (0.64 Mb of DNA assembly). Statistical analysis of the screening data from these libraries leads to different predictions as to whether the expression of enzymes needs to increase or decrease. Therefore, if genotype and media were varied separately this would lead to a suboptimal combination. This is applied to the design of a strain and media composition that increases 6-ACA from 9 to 48 mg/l in a single optimization step. This work introduces a generalizable platform to co-optimize genetic and non-genetic factors.     

A simple fitness proxy for structured populations with continuous traits, with case studies on the evolution of haplo-diploids and genetic dimorphisms

For structured populations in equilibrium with everybody born equal, ln(R (0)) is a useful fitness proxy for evolutionarily steady strategy (ESS) and most adaptive dynamics calculations, with R (0) the average lifetime number of offspring in the clonal and haploid cases, and half the average lifetime number of offspring fathered or mothered for Mendelian diploids. When individuals have variable birth states, as is, for example, the case in spatial models, R (0) is itself an eigenvalue, which usually cannot be expressed explicitly in the trait vectors under consideration. In that case, Q(Y| X):=-det (I-L(Y| X)) can often be used as fitness proxy, with L the next-generation matrix for a potential mutant characterized by the trait vector Y in the (constant) environment engendered by a resident characterized by X. If the trait space is connected, global uninvadability can be determined from it. Moreover, it can be used in all the usual local calculations like the determination of evolutionarily singular trait vectors and their local invadability and attractivity. We conclude with three extended case studies demonstrating the usefulness of Q: the calculation of ESSs under haplo-diploid genetics (I), of evolutionarily steady genetic dimorphisms (ESDs) with a priori proportionality of macro- and micro-gametic outputs (an assumption that is generally made but the fulfilment of which is a priori highly exceptional) (II), and of ESDs without such proportionality (III). These case studies should also have some interest in their own right for the spelled out calculation recipes and their underlying modelling methodology.     

The human "interferon-beta 2/hepatocyte stimulating factor/interleukin-6" gene: DNA polymorphism studies and localization to chromosome 7p21

The human interferon-beta 2 gene (IFNB2) product is identical to that for the B-cell stimulation factor-2 (BSF-2), the hybridoma growth factor (HGF) ("interleukin-6"), and the hepatocyte stimulating factor (HSF). Proteins derived from this gene mediate the plasma protein response to tissue injury (acute-phase response) and regulate the growth and differentiation of both B and T cells. By using the enzymes MspI, BstNI, and BglI, three polymorphic systems were detected with probes for the IFNB2 gene. The MspI and BglI polymorphisms are likely to be due to base pair substitutions; the BstNI polymorphism was revealed by nine other enzymes and is likely to be due to DNA insertions within 1 kb of the 3' flanking region of the gene. This region is rich in AT dinucleotides, and slippage at DNA replication may generate the insertions of between 0.07 and 0.23 kb that were observed. The polymorphic MspI site also lies within the vicinity of the fifth exon. The BglI polymorphic site is likely to lie in 5' flanking DNA. The three polymorphisms are separate, and a variety of haplotypes was observed. A low level of linkage disequilibrium exists between the MspI and the BglI alleles. MspI and BstNI polymorphisms were observed in Caucasoids, CAR Pygmies, Zaire Pygmies, Melanesians, and Chinese but at differing frequencies, and not all alleles were present in all populations. The BglI polymorphism was observed in Caucasoids and Africans only. Linkage studies involving the IFNB2 gene and 27 other chromosome 7 markers have localized it to between D7S135 and D7S370 at 7p22-p21.(ABSTRACT TRUNCATED AT 250 WORDS)     

Folic acid and chromosome breakage. IV. Variance estimates from longitudinal studies of normal individuals and their implications for sample size and power to detect differences between populations

The frequency of chromosome aberrations was studied in minimal essential medium (MEM) with and without folic acid (FA) in lymphocytes of 4 normal individuals, each sampled 12 times over a 1-year period. The cells cultured without FA had significantly more breaks and gaps. In both media about 75% of aberrations were classified as gaps. Calculations based on variance estimates suggest that the use of medium without FA could enhance the statistical power to distinguish differences in proportions of chromosome breakage between groups in the same study.     

Extended cytogenetic maps of sheep chromosome 1 and their cattle and river buffalo homoeologues: comparison with the OAR1 RH map and human chromosomes 2, 3, 21 and 1q

Cytogenetic maps are useful tools for several applications, such as the physical anchoring of linkage and RH maps or genome sequence contigs to specific chromosome regions or the analysis of chromosome rearrangements. Recently, a detailed RH map was reported in OAR1. In the present study, we selected 38 markers equally distributed in this RH map for identification of ovine genomic DNA clones within the ovine BAC library CHORI-243 using the virtual sheep genome browser and performed FISH mapping for both comparison of OAR1 and homoeologous chromosomes BBU1q-BBU6 and BTA1-BTA3 and considerably extending the cytogenetic maps of the involved species-specific chromosomes. Comparison of the resulting maps with human-identified homology with HSA2q, HSA3, HSA21 and HSA1q reveals complex chromosome rearrangements differentiating human and bovid chromosomes. In addition, we identified 2 new small human segments from HSA2q and HSA3q conserved in the telomeric regions of OAR1p and homoeologous chromosome regions of BTA3 and BBU6, and OAR1q, respectively. Evaluation of the present OAR1 cytogenetic map and the OAR1 RH map supports previous RH assignments with 2 main exceptions. The 2 loci BMS4011 and CL638002 occupy inverted positions in these 2 maps.     

Assignment of human teratocarcinoma derived growth factor (TDGF) sequences to chromosomes 2q37, 3q22, 6p25 and 19q13.1.

The human teratocarcinoma derived growth factor 1 (TDGF1) gene maps on chromosome (Chr) 3p21.3. One pseudogene (TDGF3) maps on Chr Xq21-->q22. We now report the nucleotide sequence and chromosome location of three additional TDGF pseudogenes. The three new sequences (TDGF2, TDGF4 and TDGF5) are truncated at the 5' end and have accumulated several point mutations, deletions and insertions. TDGF2, TDGF4 and TDGF6 map on Chrs 2q37, 6p25 and 3q22, respectively. Finally, Southern blot analysis of DNA from normal individuals shows a highly variable restriction pattern of the TDGF sequences.     

Heritable fragile sites on human chromosomes. X. New folate-sensitive fragile sites: 6p23, 9p21, 9q32, and 11q23

Four new folate-sensitive fragile sites are documented at 6p23, 9p21, 9q32, and 11q23. These have all been shown to be heritable except for the one at 9p21, which has been seen only in a single individual. As with the other autosomal fragile sites, these appear to be innocuous in heterozygotes.     

RFLP-based maps of the homoeologous group-6 chromosomes of wheat and their application in the tagging of Pm12, a powdery mildew resistance gene transferred from Aegilops speltoides to wheat

Genetic maps of the homoeologous group-6 chromosomes of bread wheat, Triticum aestivum, have been constructed spanning 103 cM on 6A, 90 cM on 6B and 124 cM on 6D. These maps were transferred to a Chinese Spring (CS) x line #31 cross to locate a dominant powdery mildew resistance gene, Pm12, introgressed into line #31 from Aegilops speltoides. Pm12 was shown to lie on the short arm of translocation chromosome 6BS-6SS.6SL in line #31, but could not be mapped more precisely due to the lack of recombination between the 6S Ae. speltoides segment and chromosome 6B. Possible strategies to reduce the size of the alien segment, which probably encompasses the complete long arm and more than 82% of the short arm of chromosome 6B, are discussed.     

Two human genes encoding zinc finger proteins,ZNF12 (KOX 3) and ZNF 26 (KOX 20), map to chromosomes 7p22-p21 and 12q24.33, respectively

Summary Two members of the human zinc finger Krüppel family, ZNF 12 (KOX 3) and ZNF 26 (KOX 20), have been localized by somatic cell hybrid analysis and in situ chromosomal hybridization. The presence of individual human zinc finger genes in mouse-human hybrid DNAs was correlated with the presence of specific human chromosomes or regions of chromosomes in the corresponding cell hybrids. Analysis of such mouse-human hybrid DNAs allowed the assignment of the ZNF 12 (KOX 3) gene to chromosome region 7p. The ZNF 26 (KOX 20) gene segregated with chromosome region 12q13-qter. The zinc finger genes ZNF 12 (KOX 3) and ZNF 26 (KOX 20) were localized by in situ chromosomal hybridization to human chromosome regions 7p22-21 and 12q24.33, respectively. These genes and the previously mapped ZNF 24 (KOX 17) and ZNF 29 (KOX 26) genes, are found near fragile sites.     

Population genetic structure among populations of three predaceous nabid species: Nabis alternatus Parshley, Nabis roseipennis Reuter and Nabis americoferous Carayon (Hemiptera: Nabidae)

An electrophoretic survey of 28 presumptive loci distributed among 15 enzyme systems was conducted to determine the genetic structure among local populations of three predaceous nabid species: Nabis alternatus, N. roseipennis and N.americoferous in human-disturbed habitat. For all three species, F_IS values indicated a deficiency of heterozygotes at most loci in most populations. The amount of genetic differentiation based on the standardized variance, F_ST when combined across all poulations was 0.100 for N. alternatus, 0.189 for N. roseipennis, and 0.335 for N. americoferous. When local populations were grouped according to a set of defined geographical limits, based on the naturally occurring spatial arrangement of the sampled localities, some unexpected patterns of genetic variation were revealed. The degree of genetic differentiation detected among local populations of nabids is discussed in light of the frequent disruptions and ecological changes in their ephemeral habitats.     

Limb mammary syndrome: a new genetic disorder with mammary hypoplasia, ectrodactyly, and other Hand/Foot anomalies maps to human chromosome 3q27.

We report on a large Dutch family with a syndrome characterized by severe hand and/or foot anomalies, and hypoplasia/aplasia of the mammary gland and nipple. Less frequent findings include lacrimal-duct atresia, nail dysplasia, hypohydrosis, hypodontia, and cleft palate with or without bifid uvula. This combination of symptoms has not been reported previously, although there is overlap with the ulnar mammary syndrome (UMS) and with ectrodactyly, ectodermal dysplasia, and clefting syndrome. Allelism with UMS and other related syndromes was excluded by linkage studies with markers from the relevant chromosomal regions. A genomewide screening with polymorphic markers allowed the localization of the genetic defect to the subtelomeric region of chromosome 3q. Haplotype analysis reduced the critical region to a 3-cM interval of chromosome 3q27. This chromosomal segment has not been implicated previously in disorders with defective development of limbs and/or mammary tissue. Therefore, we propose to call this apparently new disorder "limb mammary syndrome" (LMS). The SOX2 gene at 3q27 might be considered an excellent candidate gene for LMS because the corresponding protein stimulates expression of FGF4, an important signaling molecule during limb outgrowth and development. However, no mutations were found in the SOX2 open reading frame, thus excluding its involvement in LMS.     

A new methodology combining PCR, cloning, and sequencing of clones discriminated by RFLP for the study of microbial populations: application to an UASB reactor sample

This work describes a methodology combining DNA extraction, polymerase chain reaction amplification with primers targeting 16S ribosomal RNA genes, cloning, and sequencing of clones previously analyzed by restriction fragment length polymorphism (RFLP), which can be applied to study the microbial diversity in a given habitat. The methodology allows the minimization of the sequencing effort, which is particularly relevant when analyzing large numbers of clones. The methodology does not require particularly skilled personnel and can easily be adaptable to the molecular characterization of virtually any particular microbial population, provided that both adequate primers and suitable restriction enzymes for RFLP analysis of the clone library have been chosen. An example of application is presented, in which a sample taken from a continuously operating upflow anaerobic sludge blanket reactor was analyzed. RFLP analysis of the initial 162 clones with HaeIII allowed the identification of only 28 distinct profiles. As expected, identical RFLP profiles corresponded to identical nucleotide sequences.