Libin Deng

Natural selection on EPAS1 (HIF2α) associated with low hemoglobin concentration in Tibetan highlanders

By impairing both function and survival, the severe reduction in oxygen availability associated with high-altitude environments is likely to act as an agent of natural selection. We used genomic and candidate gene approaches to search for evidence of such genetic selection. First, a genome-wide allelic differentiation scan (GWADS) comparing indigenous highlanders of the Tibetan Plateau (3,200–3,500 m) with closely related lowland Han revealed a genome-wide significant divergence across eight SNPs located near EPAS1. This gene encodes the transcription factor HIF2α, which stimulates production of red blood cells and thus increases the concentration of hemoglobin in blood. Second, in a separate cohort of Tibetans residing at 4,200 m, we identified 31 EPAS1 SNPs in high linkage disequilibrium that correlated significantly with hemoglobin concentration. The sex-adjusted hemoglobin concentration was, on average, 0.8 g/dL lower in the major allele homozygotes compared with the heterozygotes. These findings were replicated in a third cohort of Tibetans residing at 4,300 m. The alleles associating with lower hemoglobin concentrations were correlated with the signal from the GWADS study and were observed at greatly elevated frequencies in the Tibetan cohorts compared with the Han. High hemoglobin concentrations are a cardinal feature of chronic mountain sickness offering one plausible mechanism for selection. Alternatively, as EPAS1 is pleiotropic in its effects, selection may have operated on some other aspect of the phenotype. Whichever of these explanations is correct, the evidence for genetic selection at the EPAS1 locus from the GWADS study is supported by the replicated studies associating function with the allelic variants.

SNP@Evolution: a hierarchical database of positive selection on the human genome.

BackgroundPositive selection is a driving force that has shaped the modern human. Recent developments in high throughput technologies and corresponding statistics tools have made it possible to conduct whole genome surveys at a population scale, and a variety of measurements, such as heterozygosity (HET), FST, and Tajimas D, have been applied to multiple datasets to identify signals of positive selection. However, great effort has been required to combine various types of data from individual sources, and incompatibility among datasets has been a common problem. SNP@Evolution, a new database which integrates multiple datasets, will greatly assist future work in this area.DescriptionAs part of our research scanning for evolutionary signals in HapMap Phase II and Phase III datasets, we built SNP@Evolution as a multi-aspect database focused on positive selection. Among its many features, SNP@Evolution provides computed FSTand HET of all HapMap SNPs, 5+ HapMap SNPs per qualified gene, and all autosome regions detected from whole genome window scanning. In an attempt to capture multiple selection signals across the genome, selection-signal enrichment strength (ES) values of HET, FST, and P-values of iHS of most annotated genes have been calculated and integrated within one frame for users to search for outliers. Genes with significant ES or P-values (with thresholds of 0.95 and 0.05, respectively) have been highlighted in color. Low diversity chromosome regions have been detected by sliding a 100 kb window in a 10 kb step. To allow this information to be easily disseminated, a graphical user interface (GBrowser) was constructed with the Generic Model Organism Database toolkit.ConclusionAvailable at http://bighapmap.big.ac.cn, SNP@Evolution is a hierarchical database focused on positive selection of the human genome. Based on HapMap Phase II and III data, SNP@Evolution includes 3,619,226/1,389,498 SNPs with their computed HET and FST, as well as qualified genes of 21,859/21,099 with ES values of HET and FST. In at least one HapMap population group, window scanning for selection signals has resulted in 1,606/10,138 large low HET regions. Among Phase II and III geographical groups, 660 and 464 regions show strong differentiation.

An unusual haplotype structure on human chromosome 8p23 derived from the inversion polymorphism

Chromosomal inversion is an important type of genomic variations involved in both evolution and disease pathogenesis. Here, we describe the refined genetic structure of a 3.8‐Mb inversion polymorphism at chromosome 8p23. Using HapMap data of 1,073 SNPs generated from 209 unrelated samples from CEPH—Utah residents with ancestry from northern and western Europe (CEU); Yoruba in Ibadan, Nigeria (YRI); and Asian (ASN) samples, which were comprised of Han Chinese from Beijing, China (CHB) and Japanese from Tokyo, Japan (JPT)—we successfully deduced the inversion orientations of all their 418 haplotypes. In particular, distinct haplotype subgroups were identified based on principal component analysis (PCA). Such genetic substructures were consistent with clustering patterns based on neighbor‐joining tree reconstruction, which revealed a total of four haplotype clades across all samples. Metaphase fluorescence in situ hybridization (FISH) in a subset of 10 HapMap samples verified their inversion orientations predicted by PCA or phylogenetic tree reconstruction. Positioning of the outgroup haplotype within one of YRI clades suggested that Human NCBI Build 36‐inverted order is most likely the ancestral orientation. Furthermore, the population differentiation test and the relative extended haplotype homozygosity (REHH) analysis in this region discovered multiple selection signals, also in a population‐specific manner. A positive selection signal was detected at XKR6 in the ASN population. These results revealed the correlation of inversion polymorphisms to population‐specific genetic structures, and various selection patterns as possible mechanisms for the maintenance of a large chromosomal rearrangement at 8p23 region during evolution. In addition, our study also showed that haplotype‐based clustering methods, such as PCA, can be applied in scanning for cryptic inversion polymorphisms at a genome‐wide scale. Hum Mutat 0, 1–8, 2008.

Whole genome HBV deletion profiles and the accumulation of preS deletion mutant during antiviral treatment

BackgroundHepatitis B virus (HBV), because of its error-prone viral polymerase, has a high mutation rate leading to widespread substitutions, deletions, and insertions in the HBV genome. Deletions may significantly change viral biological features complicating the progression of liver diseases. However, the clinical conditions correlating to the accumulation of deleted mutants remain unclear. In this study, we explored HBV deletion patterns and their association with disease status and antiviral treatment by performing whole genome sequencing on samples from 51 hepatitis B patients and by monitoring changes in deletion variants during treatment. Clone sequencing was used to analyze preS regions in another cohort of 52 patients.ResultsAmong the core, preS, and basic core promoter (BCP) deletion hotspots, we identified preS to have the highest frequency and the most complex deletion pattern using whole genome sequencing. Further clone sequencing analysis on preS identified 70 deletions which were classified into 4 types, the most common being preS2. Also, in contrast to the core and BCP regions, most preS deletions were in-frame. Most deletions interrupted viral surface epitopes, and are possibly involved in evading immuno-surveillance. Among various clinical factors examined, logistic regression showed that antiviral medication affected the accumulation of deletion mutants (OR = 6.81, 95% CI = 1.296 ~ 35.817, P = 0.023). In chronic carriers of the virus, and individuals with chronic hepatitis, the deletion rate was significantly higher in the antiviral treatment group (Fisher exact test, P = 0.007). Particularly, preS2 deletions were associated with the usage of nucleos(t)ide analog therapy (Fisher exact test, P = 0.023). Dynamic increases in preS1 or preS2 deletions were also observed in quasispecies from samples taken from patients before and after three months of ADV therapy. In vitro experiments demonstrated that preS2 deletions alone were not responsible for antiviral resistance, implying the coordination between wild type and mutant strains during viral survival and disease development.ConclusionsWe present the HBV deletion distribution patterns and preS deletion substructures in viral genomes that are prevalent in northern China. The accumulation of preS deletion mutants during nucleos(t)ide analog therapy may be due to viral escape from host immuno-surveillance.

A genome-wide association study combining pathway analysis for typical sporadic amyotrophic lateral sclerosis in Chinese Han populations

Sporadic amyotrophic lateral sclerosis (sALS) is a severe neurodegenerative disease that causes progressive motor neuron death. Although the etiology of sALS remains unknown, genetic variants are thought to predispose individuals to the disease. Several recent genome-wide association studies have identified a number of loci that increase sALS susceptibility, but these only explain a small proportion of the disease. To extend the current genetic evidence and to identify novel candidates of sALS, we performed a pooling genome-wide association study by 859,311 autosomal single-nucleotide polymorphisms of IlluminaHumanOmniZhongHua-8 combining pathway analysis in 250 typical sALS cases precluding age, clinical course, and phenotype interference and 250 control subjects from Chinese Han populations (CHP). The results revealed that 8 novel loci of 1p34.3, 3p21.1, 3p22.2, 10p15.2, 22q12.1, 3q13.11, 11q25, 12q24.33, and 5 previously reported loci of CNTN4 (kgp11325216), ATXN1 (kgp8327591), C9orf72 (kgp6016770), ITPR2 (kgp3041552), and SOD1 (kgp10760302) were associated with sALS from CHP. Furthermore, the pathway analysis based on the Gene Set Analysis Toolkit V2 showed that 10 top pathways were strongly associated with sALS from CHP, and among them, the 7 most potentially candidate pathways were phosphatidylinositol signaling system, Wnt signaling pathway, axon guidance, MAPK signaling pathway, neurotrophin signaling pathway, arachidonic acid metabolism, and T-cell receptor signaling pathway, a total of 39 significantly associate genes in 7 candidate pathways was suggested to involve in the pathogenesis of sALS from CHP. In conclusion, our results revealed several new loci and pathways related to sALS from CHP and extend the association evidence for partial loci, genes, and pathways, which were previously identified in other populations. Thus, our data provided new clues for exploring the pathogenesis of sALS.

Evolutionary selection associated with the multi-function of overlapping genes in the hepatitis B virus

It is a challenge to understand the discrete roles of each point mutation in viral evolution, but overlapping genes provide an excellent entrance for the investigation of this complicated process. We obtained 132 sequences from the largest overlapping region in the HBV genome. Based on the genetic divergence between genotypes B and C, we distinguished a set of related footprint mutations that are believed to be responsible for historical selection events. Examining the mutations in the functional domains, we found that the virus has adopted a coherent strategy in its evolutionary process that can be summarized as follows: (1) the distribution of mutations was non-random throughout the overlapping region, and more mutations were preserved in the sequence when one of the genes was under relaxed selection; (2) the viral domains were subject to different selective pressures; for instance, the PreS1 domain underwent a strict selection, whereas the overlapped Spacer domain was relatively relaxed with obvious tolerance of non-synonymous mutations with a high dN/dS ratio; (3) different selective pressures on two codon sites ultimately determined that every mutation persevered at a proper position. Taken together, the functional constraints of protein domains are believed to be primarily responsible for the different selection patterns exhibited by the distribution of mutations and amino acid changes in the region where overlapping genes reside.

Whole genome characterization of hepatitis B virus quasispecies with massively parallel pyrosequencing.

Viral quasispecies analysis is important for basic and clinical research. This study was designed to detect hepatitis B virus (HBV) genome-wide mutation profiling with detailed variant composition in individual patients, especially quasispecies evolution correlating with liver disease progression. We characterized viral populations by massively parallel pyrosequencing at whole HBV genome level in 17 patients with advanced liver disease (ALD) and 30 chronic carriers (CC). An average sequencing coverage of 2047× and 687× in ALD and CC groups, respectively, were achieved. Deep sequencing data resolved the landscapes of HBV substitutions and a more complicated quasispecies composition than previously observed. The values of substitution frequencies in quasispecies were clustered as either more than 80% or less than 20%, forming a unique U-shaped distribution pattern in both clinical groups. Furthermore, quantitative comparison of mutation frequencies of each site between two groups resulted in a spectrum of substitutions associated with liver disease progression, and among which, C2288A/T, C2304A, and A/G2525C/T were novel candidates. Moreover, distinct deletion patterns in preS, X, and C regions were shown between the two groups. In conclusion, pyrosequencing of the whole HBV genome revealed a panorama of viral quasispecies composition, characteristics of substitution distribution, and mutations correlating to severe liver disease.

Genetic Flux Between H1 and H2 Haplotypes of the 17q21.31 Inversion in European Population

The chromosome 17q21.31 inversion is a 900-kb common structural polymorphism found primarily in European population. Although the genetic flux within inversion region was assumed to be considerable suppressed, it is still unclear about the details of genetic exchange between the H1 (non-inverted sequence) and H2 (inverted sequence) haplotypes of this inversion. Here we describe a refined map of genetic exchanges between pairs of gene arrangements within the 17q21.31 region. Using HapMap phase II data of 1,546 single nucleotide polymorphisms, we successfully deduced 96 H1 and 24 H2 haplotypes in European samples by neighbor-joining tree reconstruction. Furthermore, we identified 15 and 26 candidate tracts with reciprocal and non-reciprocal genetic exchanges, respectively. In all 15 regions harboring reciprocal exchange, haplotypes reconstructed by clone sequencing did not support these exchange events, suggesting that such signals of exchange between two sister chromosomes in certain heterozygous individual were caused by phasing error regions. On the other hand, the finished clone sequencing across 4 of 26 tracts with non-reciprocal genetic flux confirmed that this kind of genetic exchange was caused by gene conversion. In summary, as crossover between pairs of gene arrangements had been considerably suppressed, gene conversion might be the most important mechanism for genetic exchange at 17q21.31.

Sterol regulatory element binding protein 1 and its target gene networks: Sterol regulatory element binding protein 1 and its target gene networks

Sterol regulatory element binding protein 1 (SREBP-1) is one of the important nuclear transcription factors. SREBP-1 can maintain lipids dynamic equilibrium by regulating the expression of enzymes required for synthesis of endogenous cholesterol, fatty acids, triglycerides and phospholipids. Anomalies of SREBP-1 and its target genes can cause a series of metabolic diseases such as insulin resistance, type Ⅱ diabetes, heart dysfunction, vascular complications and hepatic steatosis. In these years, the development of high-throughput technologies has greatly expanded our knowledge about SREBP-1 target genes and the pattern of transcriptional regulation. Here we reviewed recent research progress of SREBP-1, with a focus on the protein structure, activation process, DNA binding sites and target genes. Most importantly, we showed the transcriptional regulatory networks based on omics datasets, which will contribute to a better understanding of the role of SREBP-1 in lipid metabolism and provide new clues for the treatment of lipid metabolism disorders.

Erratum to: Polymorphism of rs3737597 in DISC1 Gene on Chromosome 1q42.2 in sALS Patients: a Chinese Han Population Case–Control Study

Although lots of genes have been revealed to relate to sporadic amyotrophic lateral sclerosis (sALS), its genetic mechanisms still need to be further explored. We aimed to search the novel genetic factors of sALS and assess their contribution. We constructed an integrative dataset based on the 3227 subsignificant genes (P value < 0.01) from two sALS-related genome-wide association studies (GWAS) (the US and Irish studies). A significant replication between both studies was confirmed by the gene set enrichment analysis in the integral level (P value < 10−4). Using the pathway overrepresentation analysis, we revealed the 34 shared Gene Ontology (GO) biological processes from the two independent studies (P value < 0.01). Among these pathways, the nervous system developmental pathway (NSD function, GO:0007399) was further supported by the previously reported genes related to sALS (P value = 3.28e−12). Importantly, four of 17 NSD-function-related target genes (disrupted-in-schizophrenia-1 (DISC1), CNTN4, NRXN3, and ERBB4) presented a considerable association with sALS in both studies. To further verify the association between the NSD function target genes and sALS, we preformed a two-stage case–control study based on 500 sALS patients and 500 controls of Chinese Han populations from mainland. A polymorphism of rs3737597 in DISC1 gene involved in the nervous system developmental pathway was closely associated with sALS. The nervous system developmental pathway is a potential pathogenesis of sALS, among them, the polymorphism of rs3737597 in DISC1 might play some roles.