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. ^

Molecular tools to improve chestnut management: El Bierzo as a case study

The European chestnut (Castanea sativa Mill.) is a multipurpose species that has been widely cultivated around the Mediterranean basin since ancient times. New varieties were brought to the Iberian Peninsula during the Roman Empire, which coexist since then with native populations that survived the last glaciation. The relevance of chestnut cultivation has being steadily growing since the Middle Ages, until the rural decline of the past century put a stop to this trend. Forest fires and diseases were also major factors. Chestnut cultivation is gaining momentum again due to its economic (wood, fruits) and ecologic relevance, and represents currently an important asset in many rural areas of Europe. In this Thesis we apply different molecular tools to help improve current management strategies. For this study we have chosen El Bierzo (Castile and Leon, NW Spain), which has a centenary tradition of chestnut cultivation and management, and also presents several unique features from a genetic perspective (next paragraph). Moreover, its nuts are widely appreciated in Spain and abroad for their organoleptic properties. We have focused our experimental work on two major problems faced by breeders and the industry: the lack of a fine-grained genetic characterization and the need for new strategies to control blight disease. To characterize with sufficient detail the genetic diversity and structure of El Bierzo orchards, we analyzed DNA from 169 trees grafted for nut production covering the entire region. We also analyzed 62 nuts from all traditional varieties. El Bierzo constitutes an outstanding scenario to study chestnut genetics and the influence of human management because: (i) it is located at one extreme of the distribution area; (ii) it is a major glacial refuge for the native species; (iii) it has a long tradition of human management (since Roman times, at least); and (iv) its geographical setting ensures an unusual degree of genetic isolation. Thirteen microsatellite markers provided enough informativeness and discrimination power to genotype at the individual level. Together with an unexpected level of genetic variability, we found evidence of genetic structure, with three major gene pools giving rise to the current population. High levels of genetic differentiation between groups supported this organization. Interestingly, genetic structure does not match with spatial boundaries, suggesting that the exchange of material and cultivation practices have strongly influenced natural gene flow. The microsatellite markers selected for this study were also used to classify a set of 62 samples belonging to all traditional varieties. We identified several cases of synonymies and homonymies, evidencing the need to substitute traditional classification systems with new tools for genetic profiling. Management and conservation strategies should also benefit from these tools. The avenue of high-throughput sequencing technologies, combined with the development of bioinformatics tools, have paved the way to study transcriptomes without the need for a reference genome. We took advantage of RNA sequencing and de novo assembly tools to determine the transcriptional landscape of chestnut in response to blight disease. In addition, we have selected a set of candidate genes with high potential for developing resistant varieties via genetic engineering. Our results evidenced a deep transcriptional reprogramming upon fungal infection. The plant hormones ET and JA appear to orchestrate the defensive response. Interestingly, our results also suggest a role for auxins in modulating such response. Many transcription factors were identified in this work that interact with promoters of genes involved in disease resistance. Among these genes, we have conducted a functional characterization of a two major thaumatin-like proteins (TLP) that belongs to the PR5 family. Two genes encoding chestnut cotyledon TLPs have been previously characterized, termed CsTL1 and CsTL2. We substantiate here their protective role against blight disease for the first time, including in silico, in vitro and in vivo evidence. The synergy between TLPs and other antifungal proteins, particularly endo-p-1,3-glucanases, bolsters their interest for future control strategies based on biotechnological approaches.

Evolution of foraging behavior in Drosophila by density-dependent selection

One of the rare examples of a single major gene underlying a naturally occurring behavioral polymorphism is the foraging locus of Drosophila melanogaster. Larvae with the rover allele, forR, have significantly longer foraging path lengths on a yeast paste than do those homozygous for the sitter allele, fors. These variants do not differ in general activity in the absence of food. The evolutionary significance of this polymorphism is not as yet understood. Here we examine the effect of high and low animal rearing densities on the larval foraging path-length phenotype and show that density-dependent natural selection produces changes in this trait. In three unrelated base populations the long path (rover) phenotype was selected for under high-density rearing conditions, whereas the short path (sitter) phenotype was selected for under low-density conditions. Genetic crosses suggested that these changes resulted from alterations in the frequency of the fors allele in the low-density-selected lines. Further experiments showed that density-dependent selection during the larval stage rather than the adult stage of development was sufficient to explain these results. Density-dependent mechanisms may be sufficient to maintain variation in rover and sitter behavior in laboratory populations.

Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Mérino ewes

Ewes from the Booroola strain of Australian Mérino sheep are characterized by high ovulation rate and litter size. This phenotype is due to the action of the FecBB allele of a major gene named FecB, as determined by statistical analysis of phenotypic data. By genetic analysis of 31 informative half-sib families from heterozygous sires, we showed that the FecB locus is situated in the region of ovine chromosome 6 corresponding to the human chromosome 4q22–23 that contains the bone morphogenetic protein receptor IB (BMPR-IB) gene encoding a member of the transforming growth factor-β (TGF-β) receptor family. A nonconservative substitution (Q249R) in the BMPR-IB coding sequence was found to be associated fully with the hyperprolificacy phenotype of Booroola ewes. In vitro, ovarian granulosa cells from FecBB/FecBB ewes were less responsive than granulosa cells from FecB+/FecB+ ewes to the inhibitory effect on steroidogenesis of GDF-5 and BMP-4, natural ligands of BMPR-IB. It is suggested that in FecBB/FecBB ewes, BMPR-IB would be inactivated partially, leading to an advanced differentiation of granulosa cells and an advanced maturation of ovulatory follicles.

Analysis of Melanoma Onset: Assessing Familial Aggregation by Using Estimating Equations and Fitting Variance Components via Bayesian Random Effects Models

We investigate whether relative contributions of genetic and shared environmental factors are associated with an increased risk in melanoma. Data from the Queensland Familial Melanoma Project comprising 15,907 subjects arising from 1912 families were analyzed to estimate the additive genetic, common and unique environmental contributions to variation in the age at onset of melanoma. Two complementary approaches for analyzing correlated time-to-onset family data were considered: the generalized estimating equations (GEE) method in which one can estimate relationship-specific dependence simultaneously with regression coefficients that describe the average population response to changing covariates; and a subject-specific Bayesian mixed model in which heterogeneity in regression parameters is explicitly modeled and the different components of variation may be estimated directly. The proportional hazards and Weibull models were utilized, as both produce natural frameworks for estimating relative risks while adjusting for simultaneous effects of other covariates. A simple Markov Chain Monte Carlo method for covariate imputation of missing data was used and the actual implementation of the Bayesian model was based on Gibbs sampling using the free ware package BUGS. In addition, we also used a Bayesian model to investigate the relative contribution of genetic and environmental effects on the expression of naevi and freckles, which are known risk factors for melanoma.

Proteomic analysis of k-casein micro-heterogeneity

The proteome of bovine milk is dominated by just six gene products that constitute approximately 95% of milk protein. Nonetheless, over 150 protein spots can be readily detected following two-dimensional electrophoresis of whole milk. Many of these represent isoforms of the major gene products produced through extensive posttranslational modification. Peptide mass fingerprinting of in-gel tryptic digests (using matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) in reflectron mode with alpha-cyano-4-hydroxycinnamic acid as the matrix) identified 10 forms of K-casein with isoelectric point (pl) values from 4.47 to 5.81, but could not distinguish between them. MALDI-TOF MS in linear mode, using sinapinic acid as the matrix, revealed a large tryptic peptide (mass > 5990 Da) derived from the C-terminus that contained all the known sites of genetic variance, phosphorylation and glycosylation. Two genetic variants present as singly or doubly phosphorylated forms could be distinguished using mass data alone. Glycoforms containing a single acidic tetrasaccharide were also identified. The differences in electrophoretic mobility of these isoforms were consistent with the addition of the acidic groups. While more extensively glycosylated forms were also observed, substantial loss of N-acetylneuraminic acid from the glycosyl group was evident in the MALDI spectra such that ions corresponding to the intact glycopeptide were not observed and assignment of the glycoforms was not possible. However, by analysing the pl shifts observed on the two-dimensional gels in conjunction with the MS data, the number of N-acetylneuraminic acid residues, and hence the glycoforms present, could be determined.

Familial aggregation of human infection with Schistosoma japonicum in the Poyang Lake region, China

Despite the success of extensive control measures that have been implemented in China for over 50 years, the number of individuals infected with Schistosoma japonicum remains high in the existing endemic areas. A variance components analysis was undertaken to estimate the heritable and environmental components that contribute to S. japonicum infection in the Poyang Lake region of Jiangxi Province, PR China. The total target population was 3148 from four separate administrative villages. Two thousand seven hundred and five of these comprised 400 families ranging in size from 3 to 188. After adjustments were made for gender, water contact and past history of having had schistosomiasis, the heritable component was estimated to account for as much as 58% of the phenotype variation under the polygenic model. Household was not shown to be an important environmental factor. Incorporating village effects indicated that the results were valid for the total population. We conclude that genetic heritability in this region is high and plays an important role in determining risk of infection with S. japonicum. (c) 2005 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Identification of membrane-bound and secreted proteins from Echinococcus granulosus by signal sequence trap

The signal sequence trap technique was applied to identify genes coding for secreted and membrane bound proteins from Echinococcus granulosus, the etiologic agent of cystic hydatid disease. An E. granulosus protoscolex cDNA library was constructed in the AP-PST vector such that randomly primed cDNAs were fused with a placental alkaline phosphatase reporter gene lacking its endogenous signal peptide. E. granulosus cDNAs encoding a functional signal peptide were selected by their ability to rescue secretion of alkaline phosphatase by COS-7 cells that had been transfected with the cDNA library. Eighteen positive clones were identified and sequenced. Their deduced amino acid sequences showed significant similarity with amino acid transporters, Krebs cycle intermediates transporters, presenilins and vacuolar protein sorter proteins. Other cDNAs encoded secreted proteins without homologues. Three sequences were transcribed antisense to E. granulosus expressed sequence tags. All the mRNAs were expressed in protoscoleces and adult worms, but some of them were not found in oncospheres. The putative E. granulosus secreted and membrane bound proteins identified are likely to play important roles in the metabolism, development and survival in the host and represent potential targets for diagnosis, drugs and vaccines against E. granulosus. (c) 2005 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Alternative splicing takes control of cytokine signaling

Thesis (Ph.D.)--University of Washington, 2016-08

Brain-derived neurotrophic factor (BDNF) gene : no major impact on antidepressant treatment response

The brain-derived neurotrophic factor (BDNF) has been suggested to play a pivotal role in the aetiology of affective disorders. In order to further clarify the impact of BDNF gene variation on major depression as well as antidepressant treatment response, association of three BDNF polymorphisms [rs7103411, Val66Met (rs6265) and rs7124442] with major depression and antidepressant treatment response was investigated in an overall sample of 268 German patients with major depression and 424 healthy controls. False discovery rate (FDR) was applied to control for multiple testing. Additionally, ten markers in BDNF were tested for association with citalopram outcome in the STAR*D sample. While BDNF was not associated with major depression as a categorical diagnosis, the BDNF rs7124442 TT genotype was significantly related to worse treatment outcome over 6 wk in major depression (p=0.01) particularly in anxious depression (p=0.003) in the German sample. However, BDNF rs7103411 and rs6265 similarly predicted worse treatment response over 6 wk in clinical subtypes of depression such as melancholic depression only (rs7103411: TTgenetic variation in BDNF and antidepressant treatment response or remission. Post-hoc analyses provide some preliminary support for a potential minor role of genetic variation in BDNF and antidepressant treatment outcome in the context of melancholic depression.

A deletion and point mutation study of the human papillomavirus type 16 major capsid gene

Recombinant human papillomavirus (HPV) virus-like particles (VLPs) made from the major capsid protein L1 are promising vaccine candidates for use as vaccines against genital and other HPV infections, and particularly against HPV-16. However, HPV-16 genotype variants have different binding affinities for neutralising mouse Mabs raised against HPV-16 L1 VLPs. This paper analyses, using a panel of well-characterised Mabs, the effects on the antigenicity of various C- and N-terminal deletants of HPV-16 L1 made in insect cells via recombinant baculovirus, of an A → T mutation at residue 266 (A266T), and of a C → G mutation at conserved position 428 (C428G). The effects of these changes on assembly of the variant L1s were studied by electron microscopy. Binding of Mab H16:E70 to A266T was reduced by almost half in comparison to wild type L1. Retention of the C-terminal region 428-483 was critical for the binding of conformation-specific Mabs (H16:V5, H16:E70, H16:U4 and H16:9A) whereas deletion of the nuclear localisation signal (NLS) or the C428G mutation or an N-terminal deletion (residues 2-9) did not affect the antigenicity. The N-terminal deletion resulted in a mixed population of 30 and 55 nm VLPs, which differs from the same construct expressed in Escherichia coli, whereas pentamer aggregates resulted from deletion of the 428-465 region or the C428G mutation. The results have implications both for considering use of single-genotype HPV vaccines, and for design of novel second-generation vaccines. © 2006 Elsevier B.V. All rights reserved.

The interleukin 1 gene cluster contains a major susceptibility locus for ankylosing spondylitis

Ankylosing spondylitis (AS) is a common and highly heritable inflammatory arthropathy. Although the gene HLA-B27 is almost essential for the inheritance of the condition, it alone is not sufficient to explain the pattern of familial recurrence of the disease. We have previously demonstrated suggestive linkage of AS to chromosome 2q13, a region containing the interleukin 1 (IL-1) family gene cluster, which includes several strong candidates for involvement in the disease. In the current study, we describe strong association and transmission of IL-1 family gene cluster single-nucleotide polymorphisms and haplotypes with AS.

Homologues of the maize rust resistance gene Rp1-D are genetically associated with a major rust resistance QTL in sorghum

As part of a comparative mapping study between sugarcane and sorghum, a sugarcane cDNA clone with homology to the maize Rp1-D rust resistance gene was mapped in sorghum. The cDNA probe hybridised to multiple loci, including one on sorghum linkage group (LG) E in a region where a major rust resistance QTL had been previously mapped. Partial sorghum Rp1-D homologues were isolated from genomic DNA of rust-resistant and -susceptible progeny selected from a sorghum mapping population. Sequencing of the Rp1-D homologues revealed five discrete sequence classes: three from resistant progeny and two from susceptible progeny. PCR primers specific to each sequence class were used to amplify products from the progeny and confirmed that the five sequence classes mapped to the same locus on LG E. Cluster analysis of these sorghum sequences and available sugarcane, maize and sorghum Rp1-D homologue sequences showed that the maize Rp1-D sequence and the partial sugarcane Rp1-D homologue were clustered with one of the sorghum resistant progeny sequence classes, while previously published sorghum Rp1-D homologue sequences clustered with the susceptible progeny sequence classes. Full-length sequence information was obtained for one member of a resistant progeny sequence class ( Rp1-SO) and compared with the maize Rp1-D sequence and a previously identified sorghum Rp1 homologue ( Rph1-2). There was considerable similarity between the two sorghum sequences and less similarity between the sorghum and maize sequences. These results suggest a conservation of function and gene sequence homology at the Rp1 loci of maize and sorghum and provide a basis for convenient PCR-based screening tools for putative rust resistance alleles in sorghum.