Cost and power analysis of linkage disequilibrium methods for locating genes influencing quantitative traits in humans

Linkage disequilibrium methods can be used to find genes influencing quantitative trait variation in humans. Linkage disequilibrium methods can require smaller sample sizes than linkage equilibrium methods, such as the variance component approach to find loci with a specific effect size. The increase in power is at the expense of requiring more markers to be typed to scan the entire genome. This thesis compares different linkage disequilibrium methods to determine which factors influence the power to detect disequilibrium. The costs of disequilibrium and equilibrium tests were compared to determine whether the savings in phenotyping costs when using disequilibrium methods outweigh the additional genotyping costs.^ Nine linkage disequilibrium tests were examined by simulation. Five tests involve selecting isolated unrelated individuals while four involved the selection of parent child trios (TDT). All nine tests were found to be able to identify disequilibrium with the correct significance level in Hardy-Weinberg populations. Increasing linked genetic variance and trait allele frequency were found to increase the power to detect disequilibrium, while increasing the number of generations and distance between marker and trait loci decreased the power to detect disequilibrium. Discordant sampling was used for several of the tests. It was found that the more stringent the sampling, the greater the power to detect disequilibrium in a sample of given size. The power to detect disequilibrium was not affected by the presence of polygenic effects.^ When the trait locus had more than two trait alleles, the power of the tests maximized to less than one. For the simulation methods used here, when there were more than two-trait alleles there was a probability equal to 1-heterozygosity of the marker locus that both trait alleles were in disequilibrium with the same marker allele, resulting in the marker being uninformative for disequilibrium.^ The five tests using isolated unrelated individuals were found to have excess error rates when there was disequilibrium due to population admixture. Increased error rates also resulted from increased unlinked major gene effects, discordant trait allele frequency, and increased disequilibrium. Polygenic effects did not affect the error rates. The TDT, Transmission Disequilibrium Test, based tests were not liable to any increase in error rates.^ For all sample ascertainment costs, for recent mutations ($<$100 generations) linkage disequilibrium tests were less expensive than the variance component test to carry out. Candidate gene scans saved even more money. The use of recently admixed populations also decreased the cost of performing a linkage disequilibrium test. ^

De Novo Assembly And Transcriptome Analysis Of The Rubber Tree (hevea Brasiliensis) And Snp Markers Development For Rubber Biosynthesis Pathways.

Hevea brasiliensis (Willd. Ex Adr. Juss.) Muell.-Arg. is the primary source of natural rubber that is native to the Amazon rainforest. The singular properties of natural rubber make it superior to and competitive with synthetic rubber for use in several applications. Here, we performed RNA sequencing (RNA-seq) of H. brasiliensis bark on the Illumina GAIIx platform, which generated 179,326,804 raw reads on the Illumina GAIIx platform. A total of 50,384 contigs that were over 400 bp in size were obtained and subjected to further analyses. A similarity search against the non-redundant (nr) protein database returned 32,018 (63%) positive BLASTx hits. The transcriptome analysis was annotated using the clusters of orthologous groups (COG), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Pfam databases. A search for putative molecular marker was performed to identify simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs). In total, 17,927 SSRs and 404,114 SNPs were detected. Finally, we selected sequences that were identified as belonging to the mevalonate (MVA) and 2-C-methyl-D-erythritol 4-phosphate (MEP) pathways, which are involved in rubber biosynthesis, to validate the SNP markers. A total of 78 SNPs were validated in 36 genotypes of H. brasiliensis. This new dataset represents a powerful information source for rubber tree bark genes and will be an important tool for the development of microsatellites and SNP markers for use in future genetic analyses such as genetic linkage mapping, quantitative trait loci identification, investigations of linkage disequilibrium and marker-assisted selection.

Multidimensional cluster stability analysis from a Brazilian Bradyrhizobium sp RFLP/PCR data set

The taxonomy of the N(2)-fixing bacteria belonging to the genus Bradyrhizobium is still poorly refined, mainly due to conflicting results obtained by the analysis of the phenotypic and genotypic properties. This paper presents an application of a method aiming at the identification of possible new clusters within a Brazilian collection of 119 Bradryrhizobium strains showing phenotypic characteristics of B. japonicum and B. elkanii. The stability was studied as a function of the number of restriction enzymes used in the RFLP-PCR analysis of three ribosomal regions with three restriction enzymes per region. The method proposed here uses Clustering algorithms with distances calculated by average-linkage clustering. Introducing perturbations using sub-sampling techniques makes the stability analysis. The method showed efficacy in the grouping of the species B. japonicum and B. elkanii. Furthermore, two new clusters were clearly defined, indicating possible new species, and sub-clusters within each detected cluster. (C) 2008 Elsevier B.V. All rights reserved.

Analysis of genomic and functional RFLP derived markers associated with sucrose content, fiber and yield QTLs in a sugarcane (Saccharum spp.) commercial cross

Expressed sequence tags derived markers have a great potential to be used in functional map construction and QTL tagging. In the present work, sugarcane genomic probes and expressed sequence tags having homology to genes, mostly involved in carbohydrate metabolism were used in RFLP assays to identify putative QTLs as well as their epistatic interactions for fiber content, cane yield, pol and tones of sugar per hectare, at two crop cycles in a progeny derived from a bi-parental cross of sugarcane elite materials. A hundred and twenty marker trait associations were found, of which 26 at both crop cycle and 32 only at first ratoon cane. A sucrose synthase derived marker was associated with a putative QTL having a high negative effect on cane yield and also with a QTL having a positive effect on Pol at both crop cycles. Fifty digenic epistatic marker interactions were identified for the four traits evaluated. Of these, only two were observed at both crop cycles.

Mendelian inheritance, linkage and linkage disequilibrium in microsatellite loci of Copaifera langsdorffii Desf.

Copaifera langsdorffii is a Neotropical tree with wide distribution in the Brazilian Atlantic rain forest and savanna. Although eight microsatellite loci (SSR) were developed in 2000 and have been widely used since then, there is yet no information about their inheritance, linkage and linkage disequilibrium (LD). Through the analysis of 28 open-pollinated (OP) progenies, the SSR loci revealed Mendelian inheritance and independent assortment. Using these progenies, young and adult trees LD was mainly detected in OP progenies. Our results show clear evidence that the eight SSR loci can be used without restriction in genetic diversity, mating system and parentage analysis.

Analysis of five polymorphic DNA markers for indirect genetic diagnosis of haemophilia A in the Brazilian population

Hemophilia A is an X-linked, inherited, bleeding disorder caused by the partial or total inactivity of the coagulation factor VIII (FVIII). Due to difficulties in the direct recognition of the disease-associated mutation in the F8 gene, indirect diagnosis using polymorphic markers located inside or close to the gene is used as an alternative for determining the segregation of the mutant gene within families and thus for detecting carrier individuals and/or assisting in prenatal diagnosis. This study characterizes the allelic and haplotype frequencies, genetic diversity, population differentiation and linkage disequilibrium of five microsatellites (F8Int1, F8Int13, F8Int22, F8Int25.3 and IKBKG) in samples of healthy individuals from Sao Paulo, Rio Grande do Sul and Pernambuco and of patients from Sao Paulo with haemophilia A to determine the degree of informativeness of these microsatellites for diagnostic purposes. The interpopulational diversity parameters highlight the differences among the analyzed population samples. Regional differences in allelic frequencies must be taken into account when conducting indirect diagnosis of haemophilia A. With the exception of IKBKG, all of the microsatellites presented high heterozygosity levels. Using the markers described, diagnosis was possible in 10 of 11 families. The F8Int22, F8Int1, F8Int13, F8Int25.3 and IKBKG microsatellites were informative in seven, six, five and two of the cases, respectively, demonstrating the effectiveness of using these microsatellites in prenatal diagnosis and in carrier identification in the Brazilian population.

Mutation analysis of the CDKN2A promoter in Australian melanoma families

Approximately 50% of all melanoma families worldwide show linkage to 9p21-22, but only about half of these have been shown to contain germ line CDKN2A mutations. It has been hypothesized that a proportion of these families carry mutations in the noncoding regions of CDKN2A. Several Canadian families have been reported to carry a mutation in the 5' UTR, at position -34 relative to the start site, which gives rise to a novel AUG translation initiation codon that markedly decreases translation from the wild-type AUG (Liu et al., 1999). Haplotype sharing in these Canadian families suggested that this mutation is of British origin. We sequenced 1,327 base pairs (bp) of CDKN2A, making up 1,116 bp of the 5' UTR and promoter, all of exon 1, and 61 bp of intron 1, in at least one melanoma case from 110 Australian families with three or more affected members known not to carry mutations within the p16 coding region. In addition, 431 bp upstream of the start codon was sequenced in an additional 253 affected probands from two-case melanoma families for which the CDKN2A mutation status was unknown. Several known polymorphisms at positions -33, -191, -493, and -735 were detected, in addition to four novel variants at positions 120, -252, -347, and -981 relative to the start codon. One of the probands from a two-case family was found to have the previously reported Q50R mutation. No family member was found to carry the mutation at position -34 or any other disease-associated mutation. For further investigation of noncoding CDKN2A mutations that may affect transcription, allele-specific expression analysis was carried out in 31 of the families with at least three affected members who showed either complete or indeterminate 9p haplotype sharing without CDKN2A exonic mutations. Reverse transcription polymerase chain reaction and automated sequencing showed expression of both CDKN2A alleles in all family members tested. The lack of CDKN2A promoter mutations and the absence of transcriptional silencing in the germ line of this cohort of families suggest that mutations in the promoter and 5' UTR play a very limited role in melanoma predisposition. (C) 2001 Wiley-Liss, Inc.

A DNA fingerprinting procedure for ultra high-throughput genetic analysis of insects

Existing procedures for the generation of polymorphic DNA markers are not optimal for insect studies in which the organisms are often tiny and background molecular Information is often non-existent. We have used a new high throughput DNA marker generation protocol called randomly amplified DNA fingerprints (RAF) to analyse the genetic variability In three separate strains of the stored grain pest, Rhyzopertha dominica. This protocol is quick, robust and reliable even though it requires minimal sample preparation, minute amounts of DNA and no prior molecular analysis of the organism. Arbitrarily selected oligonucleotide primers routinely produced similar to 50 scoreable polymorphic DNA markers, between individuals of three Independent field isolates of R. dominica. Multivariate cluster analysis using forty-nine arbitrarily selected polymorphisms generated from a single primer reliably separated individuals into three clades corresponding to their geographical origin. The resulting clades were quite distinct, with an average genetic difference of 37.5 +/- 6.0% between clades and of 21.0 +/- 7.1% between individuals within clades. As a prelude to future gene mapping efforts, we have also assessed the performance of RAF under conditions commonly used in gene mapping. In this analysis, fingerprints from pooled DNA samples accurately and reproducibly reflected RAF profiles obtained from Individual DNA samples that had been combined to create the bulked samples.

The human melanocortin-1 receptor locus: analysis of transcription unit, locus polymorphism and haplotype evolution

The complete sequence of the MCIR locus has been assembled, the coding region of the gene is intronless and placed within a 12 kb region flanked by the NULP1 and TUBB4 genes. The immediate promoter region has an E-box site with homology to the M-box consensus known to bind the microphthalmia transcription factor (MITF), however, promoter deletion analysis and transactivation studies have failed to show activation through this element by MITF. Polymorphism within the coding region, immediate 5' promoter region and a variable number tandem repeat (VNTR) minisatellite within the locus have been examined in a collection of Caucasian families and African individuals. Haplotype analysis shows linkage disequilibrium between the VNTR and MCIR coding region red hair variant alleles which can be used to estimate the age of these missense changes. Assuming a mean VNTR mutation rate of 1% and a star phylogeny, we estimate the Arg151Cys variant arose 7500 years before the present day, suggesting these variants may have arisen in the Caucasian population more recently than previously thought. (C) 2001 Published by Elsevier Science B.V.

Diversity of Phytophthora clandestina isolated from subterranean clover in southern australia: Analysis of virulence and RAPD profiles

The virulence spectrum of 112 isolates of Phytophthora clandestina collected from 56 sites in four subterranean clover-growing states in southern Australia was determined using differential cultivars of subterranean clover. Five races were detected, with race 0 in all states except New South Wales, race 1 in all states, race 2 only in Victoria, race 3 only in New South Wales, and race 4 in Victoria and Western Australia. The level of genotypic diversity among the different P. clandestina populations was investigated using five RAPD primers. Among 30 bands amplified, only two were polymorphic. This enabled identification of four multilocus RAPD genotypes. Three of the four genotypes occurred in all four states. Races 2 and 3 occurred with RAPD genotypes 1 and 2 only whereas races 0 and 1 occurred in all four multilocus RAPD genotypes. These results indicate that the pathogenicity spectrum of P. clandestina can change rapidly.

Lucerne biology and genetic improvement - an analysis of past activities and future goals in Australia

Breeding methodologies for cultivated lucerne (Medicago sativa L.), an autotetraploid, have changed little over the last 50 years, with reliance on polycross methods and recurrent phenotypic selection. There has been, however, an increase in our understanding of lucerne biology, in particular the genetic relationships between members of the M. sativa complex, as deduced by DNA analysis. Also, the differences in breeding behaviour and vigour of diploids versus autotetraploids, and the underlying genetic causes, are discussed in relation to lucerne improvement. Medicago falcata, a member of the M. sativa complex, has contributed substantially to lucerne improvement in North America, and its diverse genetics would appear to have been under-utilised in Australian programs over the last two decades, despite the reduced need for tolerance to freezing injury in Australian environments. Breeding of lucerne in Australia only commenced on a large scale in 1977, driven by an urgent need to introgress aphid resistance into adapted backgrounds. The release in the early 1980s of lucernes with multiple pest and disease resistance (aphids, Phytophthora, Colletotrichum) had a significant effect on increasing lucerne productivity and persistence in eastern Australia, with yield increases under high disease pressure of up to 300% being recorded over the predominant Australian cultivar, up to 1977, Hunter River. Since that period, irrigated lucerne yields have plateaued, highlighting the need to identify breeding objectives, technologies, and the germplasm that will create new opportunities for increasing performance. This review discusses major goals for lucerne improvement programs in Australia, and provides indications of the germplasm sources and technologies that are likely to deliver the desired outcomes.

Identification of causal relationships among traits related to drought resistance in Stylosanthes scabra using QTL analysis

Previous studies have shown that a negative relationship exists between transpiration efficiency (TE) and carbon isotope discrimination (Delta) and between TE and specific leaf area (SLA) in Stylosanthes scabra, A glasshouse experiment was conducted to confirm these relationships in an F-2 population and to study the causal nature of these relationships through quantitative trait loci (QTL) analysis, One hundred and twenty F-2 genotypes from a cross between two genotypes within S. scabra were used. Three replications for each genotype were maintained through vegetative propagation, Water stress was imposed by maintaining plants at 40% of field capacity for about 45 d. To facilitate QTL analysis, a genetic linkage map consisting of 151 RAPD markers was developed, Results from this study show that Delta was significantly and negatively correlated with TE and biomass production. Similarly, SLA showed significant negative correlation with TE and biomass production, Most of the QTL for TE and Delta were present on linkage groups 5 and 11. Similarly, QTL for SLA, transpiration and biomass productivity traits were clustered on linkage groups 13 and 24, One unlinked marker was also associated with these traits, There were several markers coincident between different traits, At all the coincident QTL, the direction of QTL effects was consistent with phenotypic data, At the coincident markers between TE and Delta, high alleles of TE were associated with low alleles of Delta. Similarly, low alleles of SLA were associated with high alleles of biomass productivity traits and transpiration. At the coincident markers between trans-4-hydroxy-N-methyl proline (MHP) and relative water content (RWC), low alleles of MHP were associated with high alleles of RWC, This study suggests the causal nature of the relationship between TE and Delta. Phenotypic data and QTL, data show that SLA was more closely associated with biomass production than with TE, This study also shows that a cause-effect relationship may exist between SLA and biomass production.

Cluster Analysis of Business Data

This journal provides immediate open access to its content on the principle that making research freely available to the public supports a greater global exchange of knowledge.