The genetics of skin, hair, and eye color variation and its relevance to forensic pigmentation predictive tests

This review examines the potential application of single nucleotide polymorphism (SNP)-based predictive tests for skin, hair, and eye color to forensic analysis in support of police investigations lacking DNA database matches or eyewitness testimony. Brief descriptions of the biology of melanogenesis and the main genes involved are presented in order to understand the basis of common pigmentation variation in humans. We outline the most recently developed forensically sensitive multiplex tests that can be applied to investigative analyses. The review also describes the biology of the SNPs with the closest associations to, and therefore the best predictors for, common variation in eye, hair, and skin pigmentation. Because pigmentation pathways are complex in their patterns, many of the better-studied human albinism traits provide insight into how pigmentation SNPs interact, control, or modify gene expression and show varying degrees of association with the key genes identified to date. These aspects of SNP action are discussed in an overview of each of the functional groups of pigmentation genes.

Population pharmacokinetic/pharmacogenetic model of lopinavir/ritonavir in HIV-infected patients

AIM

This study aims to develop a population pharmacokinetic/pharmacogenetic model for lopinavir/ritonavir (LPV/r) in European HIV-infected patients.

MATERIALS & METHODS

A total of 693 LPV/r plasma concentrations were assessed and 15 single-nucleotide polymorphisms were genotyped. The population pharmacokinetic/pharmacogenetic model was created using a nonlinear mixed-effect approach (NONMEM^® v.7.2.0., ICON Development Solutions, Dublin, Ireland).

RESULTS

Covariates significantly related to LPV/r apparent clearance (CL/F) were ritonavir trough concentration (RTC), BMI, high-density lipoprotein cholesterol (HDL-C) and certain single-nucleotide polymorphisms in genes encoding for metabolizing enzymes, which are representable as follows: CL/F = (0.216BMI + 0.0125HDL-C) × 0.713^RTC × 1.26^{rs28371764[C/T]} × 0.528^{rs6945984[C/C]} × 0.302 ^{CYP3A4[1461insA/del]} Conclusion: The LPV/r standard dose appears to be appropriate for the rs28371764[C/T] genotype. However, lower doses should be recommended for the rs6945984[C/C] and CYP3A4[1461insA/del] genotypes and even for those patients without any of these variants, as the standard dose seems to be higher than that which is required in order to achieve therapeutic levels.

A three-stage genome-wide association study identifies a susceptibility locus for late radiotherapy toxicity at 2q24.1

There is increasing evidence supporting the role of genetic variants in the development of radiation-induced toxicity. However, previous candidate gene association studies failed to elucidate the common genetic variation underlying this phenotype, which could emerge years after the completion of treatment. We performed a genome-wide association study on a Spanish cohort of 741 individuals with prostate cancer treated with external beam radiotherapy (EBRT). The replication cohorts consisted of 633 cases from the UK and 368 cases from North America. One locus comprising TANC1 (lowest unadjusted P value for overall late toxicity=6.85×10(-9), odds ratio (OR)=6.61, 95% confidence interval (CI)=2.23-19.63) was replicated in the second stage (lowest unadjusted P value for overall late toxicity=2.08×10(-4), OR=6.17, 95% CI=2.25-16.95; Pcombined=4.16×10(-10)). The inclusion of the third cohort gave unadjusted Pcombined=4.64×10(-11). These results, together with the role of TANC1 in regenerating damaged muscle, suggest that the TANC1 locus influences the development of late radiation-induced damage.

Transcriptome and genome sequencing uncovers functional variation in humans

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project--the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

Supercomplex Assembly Determines Electron Flux in the Mitochondrial Electron Transport Chain

The textbook description of mitochondrial respiratory complexes (RCs) views them as free-moving entities linked by the mobile carriers coenzyme Q (CoQ) and cytochrome c (cyt c). This model (known as the fluid model) is challenged by the proposal that all RCs except complex II can associate in supercomplexes (SCs). The proposed SCs are the respirasome (complexes I, III, and IV), complexes I and III, and complexes III and IV. The role of SCs is unclear, and their existence is debated. By genetic modulation of interactions between complexes I and III and III and IV, we show that these associations define dedicated CoQ and cyt c pools and that SC assembly is dynamic and organizes electron flux to optimize the use of available substrates.

Enhanced localization of genetic samples through linkage-disequilibrium correction

Characterizing the spatial patterns of genetic diversity in human populations has a wide range of applications, from detecting genetic mutations associated with disease to inferring human history. Current approaches, including the widely used principal-component analysis, are not suited for the analysis of linked markers, and local and long-range linkage disequilibrium (LD) can dramatically reduce the accuracy of spatial localization when unaccounted for. To overcome this, we have introduced an approach that performs spatial localization of individuals on the basis of their genetic data and explicitly models LD among markers by using a multivariate normal distribution. By leveraging external reference panels, we derive closed-form solutions to the optimization procedure to achieve a computationally efficient method that can handle large data sets. We validate the method on empirical data from a large sample of European individuals from the POPRES data set, as well as on a large sample of individuals of Spanish ancestry. First, we show that by modeling LD, we achieve accuracy superior to that of existing methods. Importantly, whereas other methods show decreased performance when dense marker panels are used in the inference, our approach improves in accuracy as more markers become available. Second, we show that accurate localization of genetic data can be achieved with only a part of the genome, and this could potentially enable the spatial localization of admixed samples that have a fraction of their genome originating from a given continent. Finally, we demonstrate that our approach is resistant to distortions resulting from long-range LD regions; such distortions can dramatically bias the results when unaccounted for.

The genetic legacy of the pre-colonial period in contemporary Bolivians

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1-2%) of a sub-Saharan component in Bolivians.

Uniparental markers of contemporary Italian population reveals details on its pre-Roman heritage

Background

According to archaeological records and historical documentation, Italy has been a melting point for populations of different geographical and ethnic matrices. Although Italy has been a favorite subject for numerous population genetic studies, genetic patterns have never been analyzed comprehensively, including uniparental and autosomal markers throughout the country.

Methods/Principal Findings

A total of 583 individuals were sampled from across the Italian Peninsula, from ten distant (if homogeneous by language) ethnic communities — and from two linguistic isolates (Ladins, Grecani Salentini). All samples were first typed for the mitochondrial DNA (mtDNA) control region and selected coding region SNPs (mtSNPs). This data was pooled for analysis with 3,778 mtDNA control-region profiles collected from the literature. Secondly, a set of Y-chromosome SNPs and STRs were also analyzed in 479 individuals together with a panel of autosomal ancestry informative markers (AIMs) from 441 samples. The resulting genetic record reveals clines of genetic frequencies laid according to the latitude slant along continental Italy – probably generated by demographical events dating back to the Neolithic. The Ladins showed distinctive, if more recent structure. The Neolithic contribution was estimated for the Y-chromosome as 14.5% and for mtDNA as 10.5%. Y-chromosome data showed larger differentiation between North, Center and South than mtDNA. AIMs detected a minor sub-Saharan component; this is however higher than for other European non-Mediterranean populations. The same signal of sub-Saharan heritage was also evident in uniparental markers.

Conclusions/Significance

Italy shows patterns of molecular variation mirroring other European countries, although some heterogeneity exists based on different analysis and molecular markers. From North to South, Italy shows clinal patterns that were most likely modulated during Neolithic times.

Genetic epistasis in female suicide attempters

BACKGROUND

Complex behaviors such as suicidal behavior likely exhibit gene-gene interactions. The main aim of this study is to explore potential single nucleotide polymorphisms combinations with epistatic effect in suicidal behavior using a data mining tool (Multifactor Dimensionality Reduction).

METHODS

Genomic DNA from peripheral blood samples was analyzed using SNPlex Technology. Multifactor Dimensionality Reduction was used to detect epistatic interactions between single nucleotide polymorphisms from the main central nervous system (CNS) neurotransmitters (dopamine: 9; noradrenaline: 19; serotonin: 23; inhibitory neurotransmitters: 60) in 889 individuals (417 men and 472 women) aged 18 years or older (585 psychiatric controls without a history of suicide attempts, and 304 patients with a history of suicide attempts). Individual analysis of association between single nucleotide polymorphisms and suicide attempts was estimated using logistic regression models.

RESULTS

Multifactor Dimensionality Reduction showed significant epistatic interactions involving four single nucleotide polymorphisms in female suicide attempters with a classification test accuracy of 60.7% (59.1%-62.4%, 95% CI): rs1522296, phenylalanine hydroxylase gene (PAH); rs7655090, dopamine receptor D5 gene (DRD5); rs11888528, chromosome 2 open reading frame 76, close to diazepam binding inhibitor gene (DBI); and rs2376481, GABA-A receptor subunit γ3 gene (GABRG3). The multivariate logistic regression model confirmed the relevance of the epistatic interaction [OR(95% CI)=7.74(4.60-13.37)] in females.

CONCLUSIONS

Our results suggest an epistatic interaction between genes of all monoamines and GABA in female suicide attempters.