PennSeq: accurate isoform-specific gene expression quantification in RNA-Seq by modeling non-uniform read distribution.

Correctly estimating isoform-specific gene expression is important for understanding complicated biological mechanisms and for mapping disease susceptibility genes. However, estimating isoform-specific gene expression is challenging because various biases present in RNA-Seq (RNA sequencing) data complicate the analysis, and if not appropriately corrected, can affect isoform expression estimation and downstream analysis. In this article, we present PennSeq, a statistical method that allows each isoform to have its own non-uniform read distribution. Instead of making parametric assumptions, we give adequate weight to the underlying data by the use of a non-parametric approach. Our rationale is that regardless what factors lead to non-uniformity, whether it is due to hexamer priming bias, local sequence bias, positional bias, RNA degradation, mapping bias or other unknown reasons, the probability that a fragment is sampled from a particular region will be reflected in the aligned data. This empirical approach thus maximally reflects the true underlying non-uniform read distribution. We evaluate the performance of PennSeq using both simulated data with known ground truth, and using two real Illumina RNA-Seq data sets including one with quantitative real time polymerase chain reaction measurements. Our results indicate superior performance of PennSeq over existing methods, particularly for isoforms demonstrating severe non-uniformity. PennSeq is freely available for download at http://sourceforge.net/projects/pennseq.

Testing genetic association with rare variants in admixed populations.

Recent studies suggest that rare variants play an important role in the etiology of many traits. Although a number of methods have been developed for genetic association analysis of rare variants, they all assume a relatively homogeneous population under study. Such an assumption may not be valid for samples collected from admixed populations such as African Americans and Hispanic Americans as there is a great extent of local variation in ancestry in these populations. To ensure valid and more powerful rare variant association tests performed in admixed populations, we have developed a local ancestry-based weighted dosage test, which is able to take into account local ancestry of rare alleles, uncertainties in rare variant imputation when imputed data are included, and the direction of effect that rare variants exert on phenotypic outcome. We used simulated sequence data to show that our proposed test has controlled type I error rates, whereas naïve application of existing rare variants tests and tests that adjust for global ancestry lead to inflated type I error rates. We showed that our test has higher power than tests without proper adjustment of ancestry. We also applied the proposed method to a candidate gene study on low-density lipoprotein cholesterol. Our results suggest that it is important to appropriately control for potential population stratification induced by local ancestry difference in the analysis of rare variants in admixed populations.

Identifying signatures of natural selection in Tibetan and Andean populations using dense genome scan data.

High-altitude hypoxia (reduced inspired oxygen tension due to decreased barometric pressure) exerts severe physiological stress on the human body. Two high-altitude regions where humans have lived for millennia are the Andean Altiplano and the Tibetan Plateau. Populations living in these regions exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. Although these responses have been well characterized physiologically, their underlying genetic basis remains unknown. We performed a genome scan to identify genes showing evidence of adaptation to hypoxia. We looked across each chromosome to identify genomic regions with previously unknown function with respect to altitude phenotypes. In addition, groups of genes functioning in oxygen metabolism and sensing were examined to test the hypothesis that particular pathways have been involved in genetic adaptation to altitude. Applying four population genetic statistics commonly used for detecting signatures of natural selection, we identified selection-nominated candidate genes and gene regions in these two populations (Andeans and Tibetans) separately. The Tibetan and Andean patterns of genetic adaptation are largely distinct from one another, with both populations showing evidence of positive natural selection in different genes or gene regions. Interestingly, one gene previously known to be important in cellular oxygen sensing, EGLN1 (also known as PHD2), shows evidence of positive selection in both Tibetans and Andeans. However, the pattern of variation for this gene differs between the two populations. Our results indicate that several key HIF-regulatory and targeted genes are responsible for adaptation to high altitude in Andeans and Tibetans, and several different chromosomal regions are implicated in the putative response to selection. These data suggest a genetic role in high-altitude adaption and provide a basis for future genotype/phenotype association studies necessary to confirm the role of selection-nominated candidate genes and gene regions in adaptation to altitude.

OPRM1 and EGFR contribute to skin pigmentation differences between Indigenous Americans and Europeans.

Contemporary variation in skin pigmentation is the result of hundreds of thousands years of human evolution in new and changing environments. Previous studies have identified several genes involved in skin pigmentation differences among African, Asian, and European populations. However, none have examined skin pigmentation variation among Indigenous American populations, creating a critical gap in our understanding of skin pigmentation variation. This study investigates signatures of selection at 76 pigmentation candidate genes that may contribute to skin pigmentation differences between Indigenous Americans and Europeans. Analysis was performed on two samples of Indigenous Americans genotyped on genome-wide SNP arrays. Using four tests for natural selection--locus-specific branch length (LSBL), ratio of heterozygosities (lnRH), Tajima's D difference, and extended haplotype homozygosity (EHH)--we identified 14 selection-nominated candidate genes (SNCGs). SNPs in each of the SNCGs were tested for association with skin pigmentation in 515 admixed Indigenous American and European individuals from regions of the Americas with high ground-level ultraviolet radiation. In addition to SLC24A5 and SLC45A2, genes previously associated with European/non-European differences in skin pigmentation, OPRM1 and EGFR were associated with variation in skin pigmentation in New World populations for the first time.

Evolution of locomotion in Anthropoidea: the semicircular canal evidence.

Our understanding of locomotor evolution in anthropoid primates has been limited to those taxa for which good postcranial fossil material and appropriate modern analogues are available. We report the results of an analysis of semicircular canal size variation in 16 fossil anthropoid species dating from the Late Eocene to the Late Miocene, and use these data to reconstruct evolutionary changes in locomotor adaptations in anthropoid primates over the last 35 Ma. Phylogenetically informed regression analyses of semicircular canal size reveal three important aspects of anthropoid locomotor evolution: (i) the earliest anthropoid primates engaged in relatively slow locomotor behaviours, suggesting that this was the basal anthropoid pattern; (ii) platyrrhines from the Miocene of South America were relatively agile compared with earlier anthropoids; and (iii) while the last common ancestor of cercopithecoids and hominoids likely was relatively slow like earlier stem catarrhines, the results suggest that the basal crown catarrhine may have been a relatively agile animal. The latter scenario would indicate that hominoids of the later Miocene secondarily derived their relatively slow locomotor repertoires.

Identifying positive selection candidate loci for high-altitude adaptation in Andean populations.

High-altitude environments (>2,500 m) provide scientists with a natural laboratory to study the physiological and genetic effects of low ambient oxygen tension on human populations. One approach to understanding how life at high altitude has affected human metabolism is to survey genome-wide datasets for signatures of natural selection. In this work, we report on a study to identify selection-nominated candidate genes involved in adaptation to hypoxia in one highland group, Andeans from the South American Altiplano. We analysed dense microarray genotype data using four test statistics that detect departures from neutrality. Using a candidate gene, single nucleotide polymorphism-based approach, we identified genes exhibiting preliminary evidence of recent genetic adaptation in this population. These included genes that are part of the hypoxia-inducible transcription factor ( HIF ) pathway, a biochemical pathway involved in oxygen homeostasis, as well as three other genomic regions previously not known to be associated with high-altitude phenotypes. In addition to identifying selection-nominated candidate genes, we also tested whether the HIF pathway shows evidence of natural selection. Our results indicate that the genes of this biochemical pathway as a group show no evidence of having evolved in response to hypoxia in Andeans. Results from particular HIF -targeted genes, however, suggest that genes in this pathway could play a role in Andean adaptation to high altitude, even if the pathway as a whole does not show higher relative rates of evolution. These data suggest a genetic role in high-altitude adaptation and provide a basis for genotype/phenotype association studies that are necessary to confirm the role of putative natural selection candidate genes and gene regions in adaptation to altitude.

Genetic evidence supports demic diffusion of Han culture.

The spread of culture and language in human populations is explained by two alternative models: the demic diffusion model, which involves mass movement of people; and the cultural diffusion model, which refers to cultural impact between populations and involves limited genetic exchange between them. The mechanism of the peopling of Europe has long been debated, a key issue being whether the diffusion of agriculture and language from the Near East was concomitant with a large movement of farmers. Here we show, by systematically analysing Y-chromosome and mitochondrial DNA variation in Han populations, that the pattern of the southward expansion of Han culture is consistent with the demic diffusion model, and that males played a larger role than females in this expansion. The Han people, who all share the same culture and language, exceed 1.16 billion (2000 census), and are by far the largest ethnic group in the world. The expansion process of Han culture is thus of great interest to researchers in many fields.

Alpha-adducin polymorphism associated with increased risk of adverse cardiovascular outcomes: results from GENEtic Substudy of the INternational VErapamil SR-trandolapril STudy (INVEST-GENES).

BACKGROUND: The alpha-adducin (ADD1) Gly460Trp polymorphism has been associated with hypertension and response to diuretic therapy, but controversy exists. METHODS: The present study was conducted to prospectively investigate the relationship among the ADD1 Gly460Trp polymorphism, diuretic use, and adverse cardiovascular outcomes among 5,979 patients with hypertensive coronary artery disease, who participated in the INVEST and provided genomic DNA. The primary outcome was defined as the first occurrence of nonfatal stroke, nonfatal myocardial infarction, or all-cause death. Secondary outcomes were the components of the primary outcome. Ancestry informative markers were used to control for population stratification. RESULTS: In blacks, ADD1 variant carriers were at higher risk for a primary outcome event than wild-type homozygotes (adjusted hazard ratio 2.62, 95% CI 1.23-5.58, P = .012), with a similar trend in whites and Hispanics, albeit a smaller magnitude of effect (adjusted hazard ratio 1.43, 0.86-2.39 in Hispanics; 1.24, 0.90-1.71 in whites). Secondary outcome analysis showed that the all-cause death was driving the differences in primary outcomes by genotype. There was no interaction between the ADD1 polymorphism and diuretic use for either primary outcome or secondary outcomes. CONCLUSIONS: In hypertensive patients with coronary artery disease, black ADD1 variant carriers showed a 2.6-fold excess risk for a primary outcome event and an 8-fold increase risk of death. White and Hispanic ADD1 variant carriers showed an increased but nonsignificant excess risk. However, the effect of diuretic use on risk of cardiovascular outcomes did not vary by ADD1 carrier status.

Finding the genes underlying adaptation to hypoxia using genomic scans for genetic adaptation and admixture mapping.

The complete sequencing the human genome and recent analytical advances have provided the opportunity to perform genome-wide studies of human variation. There is substantial potential for such population-genomic approaches to assist efforts to uncover the historical and demographic histories of human populations. Additionally, these genome-wide datasets allow for investigations of variability among genomic regions. Although all genomic regions in a population have experienced the same demographic events, they have not been affected by these events in precisely the same way. Much of the variability among genomic regions is simply the result of genetic drift (i.e., gene frequency changes resulting from the effects of small breeding-population size), but some is also the result of genetic adaptation, which will only affect the gene under selection and nearby regions. We have used a new DNA typing assay that allows for the genotyping of thousands of SNPs on hundreds of samples to identify regions most likely to have been affected by genetic adaptation. Populations that have inhabited different niches (e.g., high-altitude regions) can be used to identify genes underlying the physiological differences. We have used two methods (admixture mapping and genome scans for genetic adaptation) founded on the population-genomic paradigms to search for genes underlying population differences in response to chronic hypoxia. There is great promise that together these methods will facilitate the discovery of genes influencing hypoxic response.

[Genetic structure of Tujia as revealed by Y chromosomes].

Tujia people call themselves "Bizika", which means aboriginal. Genetic study of Tujia is virtually absent. To characterize the genetic structure of Tujia,the distribution of 14 Y haplogroups was studied in Tujia populations sampled from Enshi, Hubei (31 males) and Jishou, Hunan (68 males). A total of eight haplogroups were observed in the Enshi and Jishou populations. The haplogroup frequencies of Tujia were compared with the frequencies of other related ethnic groups, including Northern Han, Southern Han, Tibetan-Burman speaking populations,Daic and Hmong-Mien. The principal component (PC) analysis was conducted and the PCs were plotted to explore the historical migrations. In addition, partial correlation analysis was performed to study the relationship between the first three PCs and the haplogroups. The PC2 revealed a cluster of Tujia groups including Longshan,Yongshun, and Enshi with Lahu, suggesting possible interaction between Tujia and the Di-Qiang groups. However,a similarity between Han and Tujia populations, though differentiated, were also observed. We postulated, by incorporating the results of archaeological and historical evidences, that the Ba people, the ancestors of the Tujia,might be related with Di-Qiang groups and inhabited the Tujia area initially before a substantial interaction with Han and other ethnic groups.

Genetic evidence for the convergent evolution of light skin in Europeans and East Asians.

Human skin pigmentation shows a strong positive correlation with ultraviolet radiation intensity, suggesting that variation in skin color is, at least partially, due to adaptation via natural selection. We investigated the evolution of pigmentation variation by testing for the presence of positive directional selection in 6 pigmentation genes using an empirical F(ST) approach, through an examination of global diversity patterns of these genes in the Centre d'Etude du Polymorphisme Humain (CEPH)-Diversity Panel, and by exploring signatures of selection in data from the International HapMap project. Additionally, we demonstrated a role for MATP in determining normal skin pigmentation variation using admixture mapping methods. Taken together (with the results of previous admixture mapping studies), these results point to the importance of several genes in shaping the pigmentation phenotype and a complex evolutionary history involving strong selection. Polymorphisms in 2 genes, ASIP and OCA2, may play a shared role in shaping light and dark pigmentation across the globe, whereas SLC24A5, MATP, and TYR have a predominant role in the evolution of light skin in Europeans but not in East Asians. These findings support a case for the recent convergent evolution of a lighter pigmentation phenotype in Europeans and East Asians.

A genomewide admixture mapping panel for Hispanic/Latino populations.

Admixture mapping (AM) is a promising method for the identification of genetic risk factors for complex traits and diseases showing prevalence differences among populations. Efficient application of this method requires the use of a genomewide panel of ancestry-informative markers (AIMs) to infer the population of origin of chromosomal regions in admixed individuals. Genomewide AM panels with markers showing high frequency differences between West African and European populations are already available for disease-gene discovery in African Americans. However, no such a map is yet available for Hispanic/Latino populations, which are the result of two-way admixture between Native American and European populations or of three-way admixture of Native American, European, and West African populations. Here, we report a genomewide AM panel with 2,120 AIMs showing high frequency differences between Native American and European populations. The average intermarker genetic distance is ~1.7 cM. The panel was identified by genotyping, with the Affymetrix GeneChip Human Mapping 500K array, a population sample with European ancestry, a Mesoamerican sample comprising Maya and Nahua from Mexico, and a South American sample comprising Aymara/Quechua from Bolivia and Quechua from Peru. The main criteria for marker selection were both high information content for Native American/European ancestry (measured as the standardized variance of the allele frequencies, also known as "f value") and small frequency differences between the Mesoamerican and South American samples. This genomewide AM panel will make it possible to apply AM approaches in many admixed populations throughout the Americas.

Genetic structure of Hmong-Mien speaking populations in East Asia as revealed by mtDNA lineages.

Hmong-Mien (H-M) is a major language family in East Asia, and its speakers distribute primarily in southern China and Southeast Asia. To date, genetic studies on H-M speaking populations are virtually absent in the literature. In this report, we present the results of an analysis of genetic variations in the mitochondrial DNA (mtDNA) hypervariable segment 1 (HVS1) region and diagnostic variants in the coding regions in 537 individuals sampled from 17 H-M populations across East Asia. The analysis showed that the haplogroups that are predominant in southern East Asia, including B, R9, N9a, and M7, account for 63% (ranging from 45% to 90%) of mtDNAs in H-M populations. Furthermore, analysis of molecular variance (AMOVA), phylogenetic tree analysis, and principal component (PC) analysis demonstrate closer relatedness between H-M and other southern East Asians, suggesting a general southern origin of maternal lineages in the H-M populations. The estimated ages of the mtDNA lineages that are specific to H-M coincide with those based on archeological cultures that have been associated with H-M. Analysis of genetic distance and phylogenetic tree indicated some extent of difference between the Hmong and the Mien populations. Together with the higher frequency of north-dominating lineages observed in the Hmong people, our results indicate that the Hmong populations had experienced more contact with the northern East Asians, a finding consistent with historical evidence. Moreover, our data defined some new (sub-)haplogroups (A6, B4e, B4f, C5, F1a1, F1a1a, and R9c), which will direct further efforts to improve the phylogeny of East Asian mtDNAs.

SLC24A5, a putative cation exchanger, affects pigmentation in zebrafish and humans.

Lighter variations of pigmentation in humans are associated with diminished number, size, and density of melanosomes, the pigmented organelles of melanocytes. Here we show that zebrafish golden mutants share these melanosomal changes and that golden encodes a putative cation exchanger slc24a5 (nckx5) that localizes to an intracellular membrane, likely the melanosome or its precursor. The human ortholog is highly similar in sequence and functional in zebrafish. The evolutionarily conserved ancestral allele of a human coding polymorphism predominates in African and East Asian populations. In contrast, the variant allele is nearly fixed in European populations, is associated with a substantial reduction in regional heterozygosity, and correlates with lighter skin pigmentation in admixed populations, suggesting a key role for the SLC24A5 gene in human pigmentation.