A genome-wide association study identifies novel gene associations with facial skin wrinkling and mole count in Latin Americans.

BACKGROUND: Genome-wide association studies (GWASs) have identified genes influencing skin ageing and mole count in Europeans, but little is known about the relevance of these (or other genes) in non-Europeans. OBJECTIVES: To conduct a GWAS for facial skin ageing and mole count in adults < 40 years old, of mixed European, Native American and African ancestry, recruited in Latin America. METHODS: Skin ageing and mole count scores were obtained from facial photographs of over 6000 individuals. After quality control checks, three wrinkling traits and mole count were retained for genetic analyses. DNA samples were genotyped with Illumina's HumanOmniExpress chip. Association testing was performed on around 8 703 729 single-nucleotide polymorphisms (SNPs) across the autosomal genome. RESULTS: Genome-wide significant association was observed at four genome regions: two were associated with wrinkling (in 1p13·3 and 21q21·2), one with mole count (in 1q32·3) and one with both wrinkling and mole count (in 5p13·2). Associated SNPs in 5p13·2 and in 1p13·3 are intronic within SLC45A2 and VAV3, respectively, while SNPs in 1q32·3 are near the SLC30A1 gene, and those in 21q21·2 occur in a gene desert. Analyses of SNPs in IRF4 and MC1R are consistent with a role of these genes in skin ageing. CONCLUSIONS: We replicate the association of wrinkling with variants in SLC45A2, IRF4 and MC1R reported in Europeans. We identify VAV3 and SLC30A1 as two novel candidate genes impacting on wrinkling and mole count, respectively. We provide the first evidence that SLC45A2 influences mole count, in addition to variants in this gene affecting melanoma risk in Europeans.

MRPL53, a New Candidate Gene for Orofacial Clefting, Identified Using an eQTL Approach.

A valuable approach to understand how individual and population genetic differences can predispose to disease is to assess the impact of genetic variants on cellular functions (e.g., gene expression) of cell and tissue types related to pathological states. To understand the genetic basis of nonsyndromic cleft lip with or without cleft palate (NSCL/P) susceptibility, a complex and highly prevalent congenital malformation, we searched for genetic variants with a regulatory role in a disease-related tissue, the lip muscle (orbicularis oris muscle [OOM]), of affected individuals. From 46 OOM samples, which are frequently discarded during routine corrective surgeries on patients with orofacial clefts, we derived mesenchymal stem cells and correlated the individual genetic variants with gene expression from these cultured cells. Through this strategy, we detected significant cis-eQTLs (i.e., DNA variants affecting gene expression) and selected a few candidates to conduct an association study in a large Brazilian cohort (624 patients and 668 controls). This resulted in the discovery of a novel susceptibility locus for NSCL/P, rs1063588, the best eQTL for the MRPL53 gene, where evidence for association was mostly driven by the Native American ancestry component of our Brazilian sample. MRPL53 (2p13.1) encodes a 39S protein subunit of mitochondrial ribosomes and interacts with MYC, a transcription factor required for normal facial morphogenesis. Our study illustrates not only the importance of sampling admixed populations but also the relevance of measuring the functional effects of genetic variants over gene expression to dissect the complexity of disease phenotypes.

Nasal airflow simulations suggest convergent adaptation in Neanderthals and modern humans.

Both modern humans (MHs) and Neanderthals successfully settled across western Eurasian cold-climate landscapes. Among the many adaptations considered as essential to survival in such landscapes, changes in the nasal morphology and/or function aimed to humidify and warm the air before it reaches the lungs are of key importance. Unfortunately, the lack of soft-tissue evidence in the fossil record turns difficult any comparative study of respiratory performance. Here, we reconstruct the internal nasal cavity of a Neanderthal plus two representatives of climatically divergent MH populations (southwestern Europeans and northeastern Asians). The reconstruction includes mucosa distribution enabling a realistic simulation of the breathing cycle in different climatic conditions via computational fluid dynamics. Striking across-specimens differences in fluid residence times affecting humidification and warming performance at the anterior tract were found under cold/dry climate simulations. Specifically, the Asian model achieves a rapid air conditioning, followed by the Neanderthals, whereas the European model attains a proper conditioning only around the medium-posterior tract. In addition, quantitative-genetic evolutionary analyses of nasal morphology provided signals of stabilizing selection for MH populations, with the removal of Arctic populations turning covariation patterns compatible with evolution by genetic drift. Both results indicate that, departing from important craniofacial differences existing among Neanderthals and MHs, an advantageous species-specific respiratory performance in cold climates may have occurred in both species. Fluid dynamics and evolutionary biology independently provided evidence of nasal evolution, suggesting that adaptive explanations regarding complex functional phenotypes require interdisciplinary approaches aimed to quantify both performance and evolutionary signals on covariation patterns.

NAT2 gene diversity and its evolutionary trajectory in the Americas.

N-acetyltransferase 2 (NAT2) is responsible for metabolizing xenobiotics; NAT2 polymorphisms lead to three phenotypes: rapid, intermediate and slow acetylators. We aimed to investigate NAT2 diversity in Native Americans. NAT2 exon 2 was sequenced for 286 individuals from 21 populations (Native American and American Mestizos). Excluding the basal/rapid haplotype NAT2*4, the most frequent haplotypes are NAT2*5B (35.95%) in hunter-gatherers and NAT2*7B (20.61%) and NAT2*5B (19.08%) in agriculturalists that were related to the slow phenotype. A new haplotype was identified in two Amerindians. Data from the ~44 kb region surrounding NAT2 in 819 individuals from Africa, East-Asia, Europe and America were used in additional analyses. No significant differences in the acetylator NAT2 haplotype and phenotype distributions were found between Native American populations practicing farming and/or herding and those practicing hunting and gathering, probably because of the absence or weakness of selection pressures and presence of demographic and random processes preventing detection of any selection signal.

Ocular and craniofacial phenotypes in a large Brazilian family with congenital aniridia.

Congenital aniridia is a rare genetic disorder characterized by varying degrees of iris hypoplasia that are associated with additional ocular abnormalities. More than 90% of the causal mutations identified are found in the PAX6 gene, a transcription factor of critical importance in the process of neurogenesis and ocular development. Here, we investigate clinical, molecular, and craniofacial features of a large Brazilian family with congenital aniridia. Among the 56 eyes evaluated, phenotype variation encompassed bilateral total aniridia to mild iris defects with extensive variation between eyes of the same individual. PAX6 molecular screening indicated a heterozygous splice mutation (c.141 + 1G>A). Thus, we hypothesize that this splicing event may cause variation in the expression of the wild-type transcript, which may lead to the observed variation in phenotype. Affected individuals were more brachycephalic, even though their face height and cephalic circumference were not significantly different when compared to those of non-affected relatives. From this, we infer that the head shape of affected subjects may also be a result of the PAX6 splice-site mutation. Our data summarize the clinical variability associated with the ocular phenotype in a large family with aniridia, and help shed light on the role of PAX6 in neurocranial development.

T-cell receptor BV gene segment polymorphisms in healthy elderly individuals from the south of Brazil.

The ageing of the immune system (immunosenescence) is believed to be involved in both morbidity and mortality in elderly humans due to a higher incidence of infections, autoimmune diseases, cancers and other pathological situations. As any specific immune response involves recognition of antigens by T cells, the ability to develop a given immune response is also dependent on the T-cell repertoire available at a given time point. Different T-cell receptor beta variable segment (BV) (TCRBV) gene segment alleles have been associated with diseases in various human populations. In the present work we analysed the allelic frequencies of four biallelic polymorphisms in TCRBV gene segments (TCRBV3S1, TCRBV13S5, TCRBV13S6 and TCRBV18) in healthy elderly human subjects (80 years old or more) from the south of Brazil, where life expectancies reach similar levels to those observed in developed countries. Except for allele 2 of the TCRBV13S6 polymorphism, which was more frequent in elderly than in young individuals (P = 0.0105), there were no differences in allele or genotype frequencies between young and elderly individuals. The data suggest that there is no direct correlation between the TCRBV3S1, TCRBV13S5 and TCRBV18 polymorphisms analysed and healthy senescence in this particular group of elderly individuals. The higher frequency of TCRBV13S6 allele 2 in healthy elderly individuals should be confirmed in other samples to establish the significance of this finding.