Novel Statistical Methods in Quantitative Genetics : Modeling Genetic Variance for Quantitative Trait Loci Mapping and Genomic Evaluation

This thesis develops and evaluates statistical methods for different types of genetic analyses, including quantitative trait loci (QTL) analysis, genome-wide association study (GWAS), and genomic evaluation. The main contribution of the thesis is to provide novel insights in modeling genetic variance, especially via random effects models. In variance component QTL analysis, a full likelihood model accounting for uncertainty in the identity-by-descent (IBD) matrix was developed. It was found to be able to correctly adjust the bias in genetic variance component estimation and gain power in QTL mapping in terms of precision.  Double hierarchical generalized linear models, and a non-iterative simplified version, were implemented and applied to fit data of an entire genome. These whole genome models were shown to have good performance in both QTL mapping and genomic prediction. A re-analysis of a publicly available GWAS data set identified significant loci in Arabidopsis that control phenotypic variance instead of mean, which validated the idea of variance-controlling genes.  The works in the thesis are accompanied by R packages available online, including a general statistical tool for fitting random effects models (hglm), an efficient generalized ridge regression for high-dimensional data (bigRR), a double-layer mixed model for genomic data analysis (iQTL), a stochastic IBD matrix calculator (MCIBD), a computational interface for QTL mapping (qtl.outbred), and a GWAS analysis tool for mapping variance-controlling loci (vGWAS).

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.2014

Brazilian coffee genome project: an EST-based genomic resource

O café é um dos principais produtos agrícolas, sendo considerado o segundo item em importância do comércio internacional de commodities. O gênero Coffea pertence à família Rubiaceae que também inclui outras plantas importantes. Este gênero contém aproximadamente 100 espécies, mas a produção comercial é baseada somente em duas espécies, Coffea arabica e Coffea canephora, que representam aproximadamente 70 % e 30 % do mercado total de café, respectivamente. O Projeto Genoma Café Brasileiro foi desenvolvido com o objetivo de disponibilizar os modernos recursos da genômica à comunidade científica e aos diferentes segmentos da cadeia produtiva do café. Para isso, foram seqüenciados 214.964 clones escolhidos aleatoriamente de 37 bibliotecas de cDNA de C. arabica, C. canephora e C. racemosa representando estádios específicos do desenvolvimento de células e de tecidos do cafeeiro, resultando em 130.792, 12.381 e 10.566 seqüências de cada espécie, respectivamente, após processo de trimagem. Os ESTs foram agrupados em 17.982 contigs e em 32.155 singletons. A comparação destas seqüências pelo programa BLAST revelou que 22 % não tiveram nenhuma similaridade significativa às seqüências no banco de dados do National Center for Biotechnology Information (de função conhecida ou desconhecida). A base de dados de ESTs do cafeeiro resultou na identificação de cerca de 33.000 unigenes diferentes. Os resultados de anotação das seqüências foram armazenados em base de dados online em http://www.lge.ibi.unicamp.br/cafe. Os recursos desenvolvidos por este projeto disponibilizam ferramentas genéticas e genômicas que podem ser decisivas para a sustentabilidade, a competitividade e a futura viabilidade da agroindústria cafeeira nos mercados interno e externo.

Caracterização molecular e diversidade genética de diferentes variedades de abacate por marcadores microssatélites

Pertencente à família Lauraceae, o abacateiro compreende três raças hortícolas: antilhana, guatemalense e mexicana. Os marcadores moleculares são uma ferramenta rápida e eficaz para estudos genômicos, uma vez que detectam o polimorfismo diretamente ao nível do DNA e não sofrem qualquer tipo de influência ambiental. Com base nesse polimorfismo, é possível fazer inferências sobre as relações entre o genótipo e o fenótipo dos indivíduos, o que, em última análise, permite aumentar a eficiência dos programas de melhoramento. Diante o exposto, o objetivo foi investigar a diversidade genética entre sete variedades de abacate a partir de 5 lócus de marcadores moleculares microssatélites (SSR). Nas amostras de abacateiros avaliadas, encontrou-se um total de 18 alelos, com uma média de 3,6 alelos por lócus. O dendrograma gerado a partir de análise de agrupamento UPGMA agrupou, separadamente do resto dos genótipos, a cultivar Geada da raça Antilhana, possivelmente por esta variedade ser uma raça pura, e o restante foi agrupado em dois grandes grupos das raças, a Guatemalense e a Mexicana. Os genótipos das sete variedades de abacate apresentam diversidade genética nos cinco lócus de marcadores moleculares microssatélites (SSR) avaliados, o que indica que são materiais promissores para utilização em futuros programas de melhoramento.

A comparison of statistical methods for genomic selection in a mice population

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estudos ecogenômicos e bioprospectivos de Shewanella ssp

Bacteria trom Shewanella and Geobacter ganera are the most studied iron-reducing microorganisms particularly due to their electron transport systems and contribution to some industrial and environmental problems, including steel corrosion, bioenergy and bioremediation of petroleum-impacted sites. The present study was focused in two ways: the first is an in silico comparative ecogenomic study of Shewanella spp. with sequenced genomes, and the second is an experimental metagenomic work to detect iron-reducing Shewanella through PCR-DGGE of a metabolic gene. The in silico study resulted in positive correIation between copy number of 16S rDNA and genome size in Shewanella spp., with clusters of rrn near lhe origin of replication. This way, the genus is inferred as opportunist. There are no compact genomes and their sequences length varied, ranging from 4306142 nt in S. amazonensis SB2B to 5935403 nt in S. woodyi ATCC 51908, without correIation to temperature range characteristic of each specie. Intragenomic 16S rDNA sequences possess little divergence, but reasonable to resuIt in different phyIogenetic trees, depending on the sequence that is chosen to compare. For moIecuIar detection of iron-reducing Shewanella, it is proposed the mtrB gene as new biomarker. because it codes to a fundamental protein at Fe (III)-reduction. The specific primers were designed and evaluated in silico and resulted in a fragment of 360 pb. In the second study, these primers were tested in a genomic sample from S. oneidensis MR-1, amplifying the expected region. After this successfuI resuIt, the primer set was used as a tool to assess the iron-reducing communities of ShewaneIla genus under an environmental stress, i.e. crude oil contamination in mangrove sediment in Rio Grande do Norte State (Brazil). The primers presented high specificity and the reactions performed resulted in one single band of ampIification in the metagenomic samples. The fingerprinting obtained at DGGE reveaIed temporal variation of Shewanella spp. in analyzed samples. The resuIts presented show the detection of a biotechnological important group of microorganisms, the iron-reducing Shewanella spp. using a metabolic gane as target. It is concluded there are eight or more 16S rDNA sequences in Shewanella genus, with little divergence among them that affects the phylogeny; the pair of primers designed to ampIify mtrB sequences is a viable alternative to detect iron-reducing ShewanelIa in metagenomic approaches; such bacteria are present in the mangrove sediment anaIyzed, with temporal variations in the samples. This is the first experimental study that screened the iron-reducing Shewanella genus in a metagenomic experiment of mangrove sediments subjected to oil contamination through a key metabolic gene

Development and characterization of 14 microsatellite loci from an enriched genomic library of Eucalyptus camaldulensis Dehnh

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genomic instability in non-neoplastic oral mucosa cells can predict risk during 4-nitroquinoline 1-oxide-induced rat tongue carcinogenesis

4-Nitroquinotine 1-oxide (4NQO)-induced rat tongue carcinogenesis is a useful model for studying oral squamous cell carcinoma. The aim of this study was to investigate the level of DNA damage induced by 4NQO in oral mucosa cells by the single cell get (comet) assay. Mate Wistar rats were distributed into three groups of 10 animals each and treated with 50 ppm 4NQO solution by drinking water for 4, 12 or 20 weeks. Ten animals were used as negative control. Statistically significant increase of DNA damage was observed in non-neoplastic oral cells at four weeks of 4NQO administration when compared with control (P < 0.05). The level of DNA damage was directly associated with the severity of histological changes. The results suggest that histologically normal tissue is able to harbor genetically unstable cells contributing to the initiation of oral carcinogenesis. Genomic instability appears to be associated with the risk and progression of oral cancer. (C) 2004 Elsevier Ltd. All rights reserved.

Tumor microenvironmental genomic alterations in juvenile nasopharyngeal angiofibroma

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genomic instability in blood cells is able to predict the oral cancer risk: an experimental study in rats

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Novel insights into the genomic basis of citrus canker based on the genome sequences of two strains of Xanthomonas fuscans subsp aurantifolii

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Acheiropodia is caused by a genomic deletion in C7orf2, the human orthologue of the Lmbr1 gene

Acheiropodia is an autosomal recessive developmental disorder presenting with bilateral congenital amputations of the upper and lower extremities and aplasia of the hands and feet. This severely handicapping condition appears to affect only the extremities, with no other systemic manifestations reported. Recently, a locus for acheiropodia was mapped on chromosome 7q36. Herein we report the narrowing of the critical region for the acheiropodia gene and the subsequent identification of a common mutation in C7orf2-the human orthologue of the mouse Lmbr1 gene-that is responsible for the disease. Analysis of five families with acheiropodia, by means of 15 polymorphic markers, narrowed the critical region to 1.3 cM, on the basis of identity by descent, and to <0.5 Mb, on the basis of physical mapping. Analysis of C7orf2, the human orthologue of the mouse Lmbr1 gene, identified a deletion in all five families, thus identifying a common acheiropodia mutation. The deletion was identified at both the genomic-DNA and mRNA level. It leads to the production of a C7orf2 transcript lacking exon 4 and introduces a premature stop codon downstream of exon 3. Given the nature of the acheiropodia phenotype, it appears likely that the Lmbr1 gene plays an important role in limb development.

Comparative genomic hybridization detects novel amplifications in fibroadenomas of the breast

Comparative genomic hybridization analysis was performed for identification of chromosomal imbalances in 23 samples of fibroadenomas of the breast. Chromosomal gains rather than losses were a feature of these lesions. Only two cases with a familial and/or previous history of breast lesions had gain of 1q or 16q as the sole abnormality. The most frequently overrepresented segments were 5p14 (10/23 cases), 5q34-qter (6/23 cases), 13q32-qter (6/23 cases), 10q25-qter (5/23 cases), and 18q22 (4/23 cases). Some of these regions have previously been associated with breast carcinoma, but this study indicates that gain of these regions can also occur in benign breast lesions. Our findings may provide a basis for conducting further investigations to locate and identify genes associated with proliferation that may be involved in the early steps of tumorigenesis of the breast. (C) 2001 Wiley-Liss, Inc.