Genomic evidence for the Pleistocene and recent population history of Native Americans

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative “Paleoamerican” relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

Copy number variants on the X chromosome in women with primary ovarian insufficiency

Objective: To investigate whether submicroscopic copy number variants (CNVs) on the X chromosome can be identified in women with primary ovarian insufficiency (POI), defined as spontaneous secondary amenorrhea before 40 years of age accompanied by follicle-stimulating hormone levels above 40 IU/L on at least two occasions. Design: Analysis of intensity data of single nucleotide polymorphism (SNP) probes generated by genomewide Illumina 370k CNV BeadChips, followed by the validation of identified loci using a custom designed ultra-high-density comparative genomic hybridization array containing 48,325 probes evenly distributed over the X chromosome. Setting: Multicenter genetic cohort study in the Netherlands. Patient(s): 108 Dutch Caucasian women with POI, 97 of whom passed quality control, who had a normal karyogram and absent fragile X premutation, and 235 healthy Dutch Caucasian women as controls. Intervention(s): None. Main Outcome Measure(s): Amount and locus of X chromosomal microdeletions or duplications. Result(s): Intensity differences between SNP probes identify microdeletions and duplications. The initial analysis identified an overrepresentation of deletions in POI patients. Moreover, CNVs in two genes on the Xq21.3 locus (i.e., PCDH11X and TGIF2LX) were statistically significantly associated with the POI phenotype. Mean size of identified CNVs was 262 kb. However, in the validation study the identified putative Xq21.3 deletions samples did not show deviations in intensities in consecutive probes. Conclusion(s): X chromosomal submicroscopic CNVs do not play a major role in Caucasian POI patients. We provide guidelines on how submicroscopic cytogenetic POI research should be conducted. (Fertil Steril (R) 2011;95:1584-8. (C) 2011 by American Society for Reproductive Medicine.)

Multilocus Analyses of Seven Candidate Genes Suggest Interacting Pathways for Obesity-Related Traits in Brazilian Populations

We investigated whether variants in major candidate genes for food intake and body weight regulation contribute to obesity-related traits under a multilocus perspective. We studied 375 Brazilian subjects from partially isolated African-derived populations (quilombos). Seven variants displaying conflicting results in previous reports and supposedly implicated in the susceptibility of obesity-related phenotypes were investigated: beta(2)-adrenergic receptor (ADRB2) (Arg16Gly), insulin induced gene 2 (INSIG2) (rs7566605), leptin (LEP) (A19G), LEP receptor (LEPR) (Gln223Arg), perilipin (PLIN) (6209T > C), peroxisome proliferator-activated receptor-gamma (PPARG) (Pro12Ala), and resistin (RETN) (-420C > G). Regression models as well as generalized multifactor dimensionality reduction (GMDR) were employed to test the contribution of individual effects and higher-order interactions to BMI and waist-hip ratio (WHR) variation and risk of overweight/obesity. The best multilocus association signal identified in the quilombos was further examined in an independent sample of 334 Brazilian subjects of European ancestry. In quilombos, only the PPARG polymorphism displayed significant individual effects (WHR variation, P = 0.028). No association was observed either with the risk of overweight/obesity (BMI >= 25 kg/m(2)), risk of obesity alone (BMI >= 30 kg/m(2)) or BMI variation. However, GMDR analyses revealed an interaction between the LEPR and ADRB2 polymorphisms (P = 0.009) as well as a third-order effect involving the latter two variants plus INSIG2 (P = 0.034) with overweight/obesity. Assessment of the LEPR-ADRB2 interaction in the second sample indicated a marginally significant association (P = 0.0724), which was further verified to be limited to men (P = 0.0118). Together, our findings suggest evidence for a two-locus interaction between the LEPR Gln223Arg and ADRB2 Arg16Gly variants in the risk of overweight/obesity, and highlight further the importance of multilocus effects in the genetic component of obesity.

Genomic Analysis of Wild Tomato Introgressions Determining Metabolism- and Yield-Associated Traits

With the aim of determining the genetic basis of metabolic regulation in tomato fruit, we constructed a detailed physical map of genomic regions spanning previously described metabolic quantitative trait loci of a Solanum pennellii introgression line population. Two genomic libraries from S. pennellii were screened with 104 colocated markers from five selected genomic regions, and a total of 614 bacterial artificial chromosome (BAC)/cosmids were identified as seed clones. Integration of sequence data with the genetic and physical maps of Solanum lycopersicum facilitated the anchoring of 374 of these BAC/cosmid clones. The analysis of this information resulted in a genome-wide map of a nondomesticated plant species and covers 10% of the physical distance of the selected regions corresponding to approximately 1% of the wild tomato genome. Comparative analyses revealed that S. pennellii and domesticated tomato genomes can be considered as largely colinear. A total of 1,238,705 bp from both BAC/cosmid ends and nine large insert clones were sequenced, annotated, and functionally categorized. The sequence data allowed the evaluation of the level of polymorphism between the wild and cultivated tomato species. An exhaustive microsynteny analysis allowed us to estimate the divergence date of S. pennellii and S. lycopersicum at 2.7 million years ago. The combined results serve as a reference for comparative studies both at the macrosyntenic and microsyntenic levels. They also provide a valuable tool for fine-mapping of quantitative trait loci in tomato. Furthermore, they will contribute to a deeper understanding of the regulatory factors underpinning metabolism and hence defining crop chemical composition.

ADCK3, an ancestral kinase, is mutated in a form of recessive ataxia associated with coenzyme Q(10) deficiency

Muscle coenzyme Q(10) (CoQ(10) or ubiquinone) deficiency has been identified in more than 20 patients with presumed autosomal-recessive ataxia. However, mutations in genes required for CoQ(10) biosynthetic pathway have been identified only in patients with infantile-onset multisystemic diseases or isolated nephropathy. Our SNP-based genome-wide scan in a large consanguineous family revealed a locus for autosomal-recessive ataxia at chromosome 1q41. The causative mutation is a homozygous splice-site mutation in the aarF-domain-containing kinase 3 gene (ADCK3). Five additional mutations in ADCK3 were found in three patients with sporadic ataxia, including one known to have CoQ(10) deficiency in muscle. All of the patients have childhood-onset cerebellar ataxia with slow progression, and three of six have mildly elevated lactate levels. ADCK3 is a mitochondrial protein homologous to the yeast COQ8 and the bacterial UbiB proteins, which are required for CoQ biosynthesis. Three out of four patients tested showed a low endogenous pool of CoQ(10) in their fibroblasts or lymphoblasts, and two out of three patients showed impaired ubiquinone synthesis, strongly suggesting that ADCK3 is also involved in CoQ(10) biosynthesis. The deleterious nature of the three identified missense changes was confirmed by the introduction of them at the corresponding positions of the yeast COQ8 gene. Finally, a phylogenetic analysis shows that ADCK3 belongs to the family of atypical kinases, which includes phosphomositide and choline kinases, suggesting that ADCK3 plays an indirect regulatory role in ubiquinone biosynthesis possibly as part of a feedback loop that regulates ATP production.

Genetic and Functional Evidence Implicating DLL1 as the Gene That Influences Susceptibility to Visceral Leishmaniasis at Chromosome 6q27

Background. Visceral leishmaniasis (VL) is caused by Leishmania donovani and Leishmania infantum chagasi. Genome-wide linkage studies from Sudan and Brazil identified a putative susceptibility locus on chromosome 6q27. Methods. Twenty-two single-nucleotide polymorphisms (SNPs) at genes PHF10, C6orf70, DLL1, FAM120B, PSMB1, and TBP were genotyped in 193 VL cases from 85 Sudanese families, and 8 SNPs at genes PHF10, C6orf70, DLL1, PSMB1, and TBP were genotyped in 194 VL cases from 80 Brazilian families. Family-based association, haplotype, and linkage disequilibrium analyses were performed. Multispecies comparative sequence analysis was used to identify conserved noncoding sequences carrying putative regulatory elements. Quantitative reverse-transcription polymerase chain reaction measured expression of candidate genes in splenic aspirates from Indian patients with VL compared with that in the control spleen sample. Results. Positive associations were observed at PHF10, C6orf70, DLL1, PSMB1, and TBP in Sudan, but only at DLL1 in Brazil (combined P = 3 x 10(-4) at DLL1 across Sudan and Brazil). No functional coding region variants were observed in resequencing of 22 Sudanese VL cases. DLL1 expression was significantly (P = 2 x 10(-7)) reduced (mean fold change, 3.5 [SEM, 0.7]) in splenic aspirates from patients with VL, whereas other 6q27 genes showed higher levels (1.27 x 10(-6) < P < .01) than did the control spleen sample. A cluster of conserved noncoding sequences with putative regulatory variants was identified in the distal promoter of DLL1. Conclusions. DLL1, which encodes Delta-like 1, the ligand for Notch3, is strongly implicated as the chromosome 6q27 VL susceptibility gene.

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

Modelling dominance in a flexible intercross analysis

The aim of this paper is to develop a flexible model for analysis of quantitative trait loci (QTL) in outbred line crosses, which includes both additive and dominance effects. Our flexible intercross analysis (FIA) model accounts for QTL that are not fixed within founder lines and is based on the variance component framework. Genome scans with FIA are performed using a score statistic, which does not require variance component estimation. RESULTS: Simulations of a pedigree with 800 F2 individuals showed that the power of FIA including both additive and dominance effects was almost 50% for a QTL with equal allele frequencies in both lines with complete dominance and a moderate effect, whereas the power of a traditional regression model was equal to the chosen significance value of 5%. The power of FIA without dominance effects included in the model was close to those obtained for FIA with dominance for all simulated cases except for QTL with overdominant effects. A genome-wide linkage analysis of experimental data from an F2 intercross between Red Jungle Fowl and White Leghorn was performed with both additive and dominance effects included in FIA. The score values for chicken body weight at 200 days of age were similar to those obtained in FIA analysis without dominance. CONCLUSION: We have extended FIA to include QTL dominance effects. The power of FIA was superior, or similar, to standard regression methods for QTL effects with dominance. The difference in power for FIA with or without dominance is expected to be small as long as the QTL effects are not overdominant. We suggest that FIA with only additive effects should be the standard model to be used, especially since it is more computationally efficient.

Genetic studies of body mass index yield new insights for obesity biology

Obesity is heritable and predisposes to many diseases. To understand the genetic basis of obesity better, here we conduct a genome-wide association study and Metabochip meta-analysis of body mass index (BMI), a measure commonly used to define obesity and assess adiposity, in upto 339,224 individuals. This analysis identifies 97 BMI-associated loci (P < 5 x 10(-8)), 56 of which are novel. Five loci demonstrate clear evidence of several independent association signals, and many loci have significant effects on other metabolic phenotypes. The 97 loci account for similar to 2.7% of BMI variation, and genome-wide estimates suggest that common variation accounts for >20% of BMI variation. Pathway analyses provide strong support for a role of the central nervous systemin obesity susceptibility and implicate new genes and pathways, including those related to synaptic function, glutamate signalling, insulin secretion/action, energy metabolism, lipid biology and adipogenesis.

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

Fatores de risco individuais e familiares no desenvolvimento da pré-eclâmpia

Preeclampsia is a spectral disease, with different clinical forms which can evolve with severe multisystemic complications. This present study aimed to determine the risk factors associated with preeclampsia (PE); to validate the existence of aggregation of hypertensive disease in families of women with preeclampsia and verify the existence of association between polymorphisms in the VEGF gene and level of VEGF and its soluble receptor (sFlt1). A case-control study was performed (n = 851). Genotyping of VEGF was performed and serum levels of VEGF and sFlt1 were measured by ELISA. It was observed that 38% of mothers (173, 455) of a case of preeclampsia and 30.8% (78 of 361) of controls had history of hypertension (p <0.0001). Similarly, when examining the history of maternal preeclampsia, we observed that 14.6% (48 of 328) of mothers of women with preeclampsia and 9.6% (12 of 294) of mothers of controls had a history of preeclampsia (p = 0.0001). As for maternal history of preeclampsia, we found that 5.1% (15 of 295) of cases and 3.6% (7 of 314) of controls had a history of preeclampsia (p = 0.0568). Sisters of women with preeclampsia also had a history of hypertensive disease in 9% (41 of 455) versus 6.6% (13 of 361), p = 0.002. Similarly when examining the history of preeclampsia in sisters, it was observed that 22.7% (57 of 251) of a sister of case versus 11.4% (26 of 228) of controls had a history of preeclampsia (P = 0.0011). We observed a decrease in free VEGF in the serum of patients (P <0.05) and increased soluble VEGF receptor. There was no association between polymorphisms in the VEGF gene and preeclampsia. The data obtained in this work validate that hypertensive disease in mothers and sisters with preeclampsia are risk factors for preeclampsia. The risk of illness in the family is higher according to disease severity. High incidence of preeclampsia can be assumed by the high incidence of this disease among the controls. Significant differences between the frequency of preeclampsia in mothers of cases and controls indicate familial factors. Work is being conducted with the to eventually perform genome wide association studies to identify susceptibility loci

Xylella fastidiosa gene expression analysis by DNA microarrays

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

Identification of a Microdeletion at the 7q33-q35 Disrupting the CNTNAP2 Gene in a Brazilian Stuttering Case

Speech and language disorders are some of the most common referral reasons to child development centers accounting for approximately 40% of cases. Stuttering is a disorder in which involuntary repetition, prolongation, or cessation of the sound precludes the flow of speech. About 5% of individuals in the general population have a stuttering problem, and about 80% of the affected children recover naturally. The causal factors of stuttering remain uncertain in most cases; studies suggest that genetic factors are responsible for 70% of the variance in liability for stuttering, whereas the remaining 30% is due to environmental effects supporting a complex cause of the disorder. The use of high-resolution genome wide array comparative genomic hybridization has proven to be a powerful strategy to narrow down candidate regions for complex disorders. We report on a case with a complex set of speech and language difficulties including stuttering who presented with a 10Mb deletion of chromosome region 7q33-35 causing the deletion of several genes and the disruption of CNTNAP2 by deleting the first three exons of the gene. CNTNAP2 is known to be involved in the cause of language and speech disorders and autism spectrum disorder and is in the same pathway as FOXP2, another important language gene, which makes it a candidate gene for causal studies speech and language disorders such as stuttering. (C) 2010 Wiley-Liss, Inc.

Microallelotyping defines novel regions of loss of heterozygosity in uterine leiomyomas

Uterine leiomyomas are extremely common, benign, smooth muscle tumors that represent a significant public health problem. Although there have been few molecular studies of uterine leiomyomas, most of them have reported a very low frequency of loss of heterozygosity (LOH) in different regions of the genome. The detection of LOH has been used to identify genomic regions that harbor tumor suppressor genes and to characterize different tumor types. We have used a set of 15 microsatellite polymorphism markers to examine the frequency of allele loss in a panel of 64 human uterine leiomyomas matched to normal DNAs. The markers were chosen from regions involved in losses identified by comparative genomic hybridization in a subset of uterine leiomyomas described in a previous report. DNA from tumors and normal tissue was amplified by the polymerase chain reaction and subsequently analyzed using an ABI Prism 377 DNA automated sequencer. The frequency of LOH observed was low, except for the markers D15S87 (15q26.3), D7S493 (7p15.3), and D7S517 (7p22.2). No changes in microsatellite size were detected in our samples. These results provide useful clues for identifying putative tumor suppressor genes associated with a subset of uterine leiomyomas. (C) 2004 Wiley-Liss, Inc.

QTL for percentage of carcass and carcass parts in a broiler x layer cross

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