中国人群参考基因组及基因组变异图谱资源库

随着人类基因组计划和国际千人基因组计划的实施,已公开数百个中国人个体的全基因组数据。建立高精度的中国人群参考基因组序列,发现并解析中国人群特有的序列变异,是我国未来精准医学研究的基础。为满足未来精准医学研究中国人基因组数据持续增长的科学管理和深入研究的需求,中国科学院北京基因组研究所发展并建立了基于中国人群全基因组测序数据的虚拟中国人基因组数据库(Virtual Chinese Genome Database, VCGDB)和中国人群基因组变异数据库(Genome Variation Map, GVM),面向国内外用户提供数据检索、共享、下载和在线分析服务。本文重点介绍了这两个数据库的特点和功能,以及未来发展与应用前景,以期为中国人群参考基因组及基因组变异图谱资源库的推广使用、发展完善提供有益信息。     

Sharing reference data and including cows in the reference population improve genomic predictions in Danish Jersey

Small reference populations limit the accuracy of genomic prediction in numerically small breeds, such like Danish Jersey. The objective of this study was to investigate two approaches to improve genomic prediction by increasing size of reference population in Danish Jersey. The first approach was to include North American Jersey bulls in Danish Jersey reference population. The second was to genotype cows and use them as reference animals. The validation of genomic prediction was carried out on bulls and cows, respectively. In validation on bulls, about 300 Danish bulls (depending on traits) born in 2005 and later were used as validation data, and the reference populations were: (1) about 1050 Danish bulls, (2) about 1050 Danish bulls and about 1150 US bulls. In validation on cows, about 3000 Danish cows from 87 young half-sib families were used as validation data, and the reference populations were: (1) about 1250 Danish bulls, (2) about 1250 Danish bulls and about 1150 US bulls, (3) about 1250 Danish bulls and about 4800 cows, (4) about 1250 Danish bulls, 1150 US bulls and 4800 Danish cows. Genomic best linear unbiased prediction model was used to predict breeding values. De-regressed proofs were used as response variables. In the validation on bulls for eight traits, the joint DK-US bull reference population led to higher reliability of genomic prediction than the DK bull reference population for six traits, but not for fertility and longevity. Averaged over the eight traits, the gain was 3 percentage points. In the validation on cows for six traits (fertility and longevity were not available), the gain from inclusion of US bull in reference population was 6.6 percentage points in average over the six traits, and the gain from inclusion of cows was 8.2 percentage points. However, the gains from cows and US bulls were not accumulative. The total gain of including both US bulls and Danish cows was 10.5 percentage points. The results indicate that sharing reference data and including cows in reference population are efficient approaches to increase reliability of genomic prediction. Therefore, genomic selection is promising for numerically small population.     

Derivation and characterization of human embryonic stem cell lines from the Chinese population

Human embryonic stem cells (hESCs) can self-renew indefinitely and differentiate into all cell types in the human body. Therefore, they are valuable in regenerative medicine, human developmental biology and drug discovery. A number of hESC lines have been derived from the Chinese population,but limited of them are available for research purposes. Here we report the derivation and characterization of two hESC lines derived from human blastocysts of Chinese origin. These hESCs express alkaline phosphatase and hESC-specific markers, including Oct4, Nanog, SSEA-3, SSEA-4,TRA-1-60 and TRA-1-81. They also have high levels of telomerase activity and normal karyotypes. These cells can form embryoid body in vitro and can be differentiated into all three germ layers in vivo by teratoma formation. The newly established hESCs will be distributed for research purposes.The availability of hESC lines from the Chinese population will facilitate studies on the differences in hESCs from different ethnic groups.     

Derivation and characterization of human embryonic stem cell lines from the Chinese population

Human embryonic stem cells(hESCs) can self-renew indefinitely and differentiate into all cell types in the human body.Therefore,they are valuable in regenerative medicine,human developmental biology and drug discovery.A number of hESC lines have been derived from the Chinese population, but limited of them are available for research purposes.Here we report the derivation and characterization of two hESC lines derived from human blastocysts of Chinese origin.These hESCs express alkaline phosphatase and hE...     

widgetcon: A website and program for quick conversion among common population genetic data formats

One of the most tedious steps in genetic data analyses is the reformatting data generated with one program for use with other applications. This conversion is necessary because comprehensive evaluation of the data may be based on different algorithms included in diverse software, each requiring a distinct input format. A platform‐independent and freely available program or a web‐based tool dedicated to such reformatting can save time and efforts in data processing. Here, we report widgetcon , a website and a program which has been developed to quickly and easily convert among various molecular data formats commonly used in phylogenetic analysis, population genetics, and other fields. The web‐based service is available at https://www.widgetcon.net . The program and the website convert the major data formats in four basic steps in less than a minute. The resource will be a useful tool for the research community and can be updated to include more formats and features in the future.     

Symmetric tests and confidence intervals for survival probabilities and quantiles of censored survival data

We describe existing tests and introduce two new tests concerning the value of a survival function. These tests may be used to construct a confidence interval for the survival probability at a given time or for a quantile of the survival distribution. Simulation studies show that error rates can differ substantially from their nominal values, particularly at survival probabilities close to zero or one. We recommend our new constrained bootstrap test for its good overall performance.     

Use of field data in pig genomic selection schemes: a simulation study

The aim of this study was to test how genetic gain for a trait not measured on the nucleus animals could be obtained within a genomic selection pig breeding scheme. Stochastic simulation of a pig breeding program including a breeding nucleus, a multiplier to produce and disseminate semen and a production tier where phenotypes were recorded was performed to test (1) the effect of obtaining phenotypic records from offspring of nucleus animals, (2) the effect of genotyping production animals with records for the purpose of including them in a genomic selection reference population or (3) to combine the two approaches. None of the tested strategies affected genetic gain if the trait under investigation had a low economic value of only 10% of the total breeding goal. When the relative economic weight was increased to 30%, a combination of the methods was most effective. Obtaining records from offspring of already genotyped nucleus animals had more impact on genetic gain than to genotype more distant relatives with phenotypes to update the reference population. When records cannot be obtained from offspring of nucleus animals, genotyping of production animals could be considered for traits with high economic importance.     

大中型兽类种群数量估算的研究进展

大中型兽类种群数量的估算是动物生态学中重要的基本问题, 受到研究者、管理者和公众的共同关注。国际上从20世纪中期开始研究该问题, 已出现了多种研究方法和相应案例, 且还在快速发展, 但世界各地仍有很多物种的种群数量尚未知晓。在我国, 从20世纪80年代开始调查大中型兽类种群数量, 取得了重要进展, 也还有很多物种的种群数量尚不清楚。因此, 有必要归纳国际上种群数量估算的研究进展, 同时, 总结国内研究的现状、优势和趋势, 供研究者参考。本文首先选择估算大中型兽类种群数量的原理、数据来源和模型这3个要素归纳出简明的研究框架, 将现有的多种方法置于其中予以阐述。在该框架下, 根据估算原理分为4大类方法, 为距离取样法、标志重捕法、基于遇见率法和遥感影像直接计数法。针对每一大类方法, 论述其基本原理模型和模型假设, 说明能实现该原理的相应数据来源(视觉观测、红外相机拍摄、DNA微卫星识别、卫星定位跟踪、声音监测或遥感影像)的特点及如何实现该原理, 评价其适用性及优缺点, 并选择其中具有可比性的方法予以比较评价。其次, 参照该研究框架, 总结我国的研究现状, 分析未来发展的优势和趋势: 我国的红外相机数据积累充分, 可以发展以此为数据源的距离取样法、标志重捕法和基于遇见率法; 发展以粪便样品为数据来源的距离取样法和粪便DNA标志重捕法; 相比地面调查数据, 获取高分辨率遥感影像数据更容易, 尽量以此估算符合适用条件的大中型兽类的种群数量。最后, 本文提出了适用于我国大中型兽类种群数量的估算方法的选择流程, 供研究者参考。     

The macroecology of population dynamics: taxonomic and biogeographic patterns in population cycles

Regular cycles in population abundance are fascinating phenomena, but are they common in natural populations? How are they distributed among taxa? Are there differences between different regions of the world, or along latitudinal gradients? Using the new Global Population Dynamics Database we analysed nearly 700 long (25 + years) time series of animal field populations, looking for large-scale patterns in cycles. Nearly 30% of the time series were cyclic. Cycle incidence varied among taxonomic classes, being most common in mammal and fish populations, but only in fish did cycle incidence vary among orders. Cycles were equally common in European and North American populations, but were more common in Atlantic fish than Pacific fish. The incidence of cycles increased with latitude in mammals only. There was no latitudinal gradient in cycle period, but cycle amplitude declined with latitude in some groups of fish. Even after considering the biases in the data source and expected type I error, population cycles seem common enough to warrant ecological attention.     

Accounting for missing data in the estimation of contemporary genetic effective population size (Ne)

Theoretical models are often applied to population genetic data sets without fully considering the effect of missing data. Researchers can deal with missing data by removing individuals that have failed to yield genotypes and/or by removing loci that have failed to yield allelic determinations, but despite their best efforts, most data sets still contain some missing data. As a consequence, realized sample size differs among loci, and this poses a problem for unbiased methods that must explicitly account for random sampling error. One commonly used solution for the calculation of contemporary effective population size (N_e) is to calculate the effective sample size as an unweighted mean or harmonic mean across loci. This is not ideal because it fails to account for the fact that loci with different numbers of alleles have different information content. Here we consider this problem for genetic estimators of contemporary effective population size (N_e). To evaluate bias and precision of several statistical approaches for dealing with missing data, we simulated populations with known N_e and various degrees of missing data. Across all scenarios, one method of correcting for missing data (fixed‐inverse variance‐weighted harmonic mean) consistently performed the best for both single‐sample and two‐sample (temporal) methods of estimating N_e and outperformed some methods currently in widespread use. The approach adopted here may be a starting point to adjust other population genetics methods that include per‐locus sample size components.     

Analysis of genetic diversity and population structure of Chinese yak breeds （Bos grunniens） using microsatellite markers

Nine Chinese yak breeds (Maiwa,Tianzhu White,Qinghai Plateau,Sibu,Zhongdian,Pall,Tibetan High Mountain,Jiulong,and Xin-jiang) and Gayal were analyzed by means of 16 microsatellite markers to determine the level of genetic variation within populations,genetic relationship between populations,and population structure for each breed.A total of 206 microsatellite alleles were observed.Mean F-statistics (0.056) for 9 yak breeds indicated that 94.4% of the genetic variation was observed within yak breeds and 5.6% of the genetic variation existed amongst breeds.The Neighbor-Joining phylogenetic free was constructed based on Nei's standard genetic dis-tances and two clusters were obtained.The Gayal separated from the yaks far away and formed one cluster and 9 yak breeds were grouped together.The analysis of population structure for 9 yak breeds and the Gayal showed that they resulted in four clusters; one clus-ter includes yaks from Tibet Autonomous Region and Qinghai Province,one cluster combines Zhongdian,Maiwa,and Tianzhu White,and Jiulong and Xinjiang come into the third cluster.Pali was mainly in the first cluster (90%),Jiulong was mainly in the second cluster (87.1%),Zhongdian was primarily in the third cluster (83%),and the other yak breeds were distributed in two to three clusters.The Gayal was positively left in the fourth cluster (99.3%).     

Cell type specific DNA methylation in cord blood: A 450K-reference data set and cell count-based validation of estimated cell type composition

Epigenome-wide association studies of prenatal exposure to different environmental factors are becoming increasingly common. These studies are usually performed in umbilical cord blood. Since blood comprises multiple cell types with specific DNA methylation patterns, confounding caused by cellular heterogeneity is a major concern. This can be adjusted for using reference data consisting of DNA methylation signatures in cell types isolated from blood. However, the most commonly used reference data set is based on blood samples from adult males and is not representative of the cell type composition in neonatal cord blood. The aim of this study was to generate a reference data set from cord blood to enable correct adjustment of the cell type composition in samples collected at birth. The purity of the isolated cell types was very high for all samples (>97.1%), and clustering analyses showed distinct grouping of the cell types according to hematopoietic lineage. We explored whether this cord blood and the adult peripheral blood reference data sets impact the estimation of cell type composition in cord blood samples from an independent birth cohort (MoBa, n = 1092). This revealed significant differences for all cell types. Importantly, comparison of the cell type estimates against matched cell counts both in the cord blood reference samples (n = 11) and in another independent birth cohort (Generation R, n = 195), demonstrated moderate to high correlation of the data. This is the first cord blood reference data set with a comprehensive examination of the downstream application of the data through validation of estimated cell types against matched cell counts.     

Resequencing studies of nonmodel organisms using closely related reference genomes: optimal experimental designs and bioinformatics approaches for population genomics

Decreasing costs of next‐generation sequencing (NGS) experiments have made a wide range of genomic questions open for study with nonmodel organisms. However, experimental designs and analysis of NGS data from less well‐known species are challenging because of the lack of genomic resources. In this work, we investigate the performance of alternative experimental designs and bioinformatics approaches in estimating variability and neutrality tests based on the site‐frequency‐spectrum (SFS) from individual resequencing data. We pay particular attention to challenges faced in the study of nonmodel organisms, in particular the absence of a species‐specific reference genome, although phylogenetically close genomes are assumed to be available. We compare the performance of three alternative bioinformatics approaches – genotype calling, genotype–haplotype calling and direct estimation without calling genotypes. We find that relying on genotype calls provides biased estimates of population genetic statistics at low to moderate read depth (2–8×). Genotype–haplotype calling returns more accurate estimates irrespective of the divergence to the reference genome, but requires moderate depth (8–20×). Direct estimation without calling genotypes returns the most accurate estimates of variability and of most SFS tests investigated, including at low read depth (2–4×). Studies without species‐specific reference genome should thus aim for low read depth and avoid genotype calling whenever individual genotypes are not essential. Otherwise, aiming for moderate to high depth at the expense of number of individuals, and using genotype–haplotype calling, is recommended.     

proseq: A software for preparation and evolutionary analysis of DNA sequence data sets

proseq is an integrated user‐friendly windows based program for convenient sequence editing and evolutionary analysis. It is designed to simplify preparation and analysis of DNA sequence data sets in population genetic, phylogenetic and molecular ecology studies. Sequence editor features include editing of chromatogram files, contig assembly, sequence alignment, translation and other utilities. Analysis features include calculation of genetic diversity, divergence, population subdivision and gene flow with permutation‐based significance testing and various tests of neutrality. A tool for coalescent simulations implements models with intragenic recombination, population subdivision and population growth.     

Genetic analysis of 15 STR loci in Chinese Han population from West China

Allele frequencies for 15 short tandem repeat （STR） loci （D8S1179, D21S11, D7S820, CSFIPO, D3S1358, TH01, D13S317, D16S539, D2S1338, D19S433, vWA, TPOX, D18S51, D5S818, and FGA） were obtained from 7,636 unrelated individuals of Chinese Han population living in Qinghai and Chongqing, China. Totally 206 alleles were observed, with the corresponding allele frequencies ranging from 0.0001-0.4982. Chi-square test showed that all of the STR loci agreed with the Hardy-Weinberg equilibrium. We also compared our data with previously published population data of other ethnics or areas. The results are valuable for human identification and paternity testing in Chinese Han population.     

安徽省扬子鳄饲养种群繁殖力的参数分析及种群增长预测

通过对宣城扬子鳄繁殖研究中心三代繁殖鳄1982—2003年繁殖力有关参数的分析,结果表明近三年来繁殖研究中心扬子鳄产卵量明显减少;亲代受精率为(83.80±2.37)%、孵化率为(89.53±0.86)%;发现F1受精率为(79.38±2.74)%、孵化率为(83.78±1.95)%;F2代个体受精率为(68.7±1.84)%、孵化率为(88.16±1.68)%;通过单因素方差分析F(受精率)=4.33(P<0.05)、F(孵化率)=4.56(P<0.05)差异显著。分析影响扬子鳄繁殖的因素发现:产卵量明显减少可能与亲代鳄逐渐退出繁殖期、F1代鳄整体产卵能力有下降趋势、F2代鳄没有大规模进入繁殖期有关;分析历史资料发现扬子鳄产卵窝数与4月份的温度相关性高(R=0.979,P<0.01);作者认为F1、F2代个体与亲代之间的受精率、孵化率有显著差异,近交衰退、遗传多样性的丧失有密切关系。利用时间序列分析建立理想状态下饲养扬子鳄种群数量估计模型,预测10年内扬子鳄饲养将可能达到19000尾;根据该种群增长模式及种群年龄结构特点,应扩大饲养规模,特别是育成鳄饲养池及幼鳄饲养池的建设。     

Evaluation of polymorphisms in microRNA‐binding sites and pancreatic cancer risk in Chinese population

As promising biomarkers and therapy targets, microRNAs (miRNAs) are involved in various physiological and tumorigenic processes. Genetic variants in miRNA‐binding sites can lead to dysfunction of miRNAs and contribute to disease. However, systematic investigation of the miRNA‐related single nucleotide polymorphisms (SNPs) for pancreatic cancer (PC) risk remains elusive. We performed integrative bioinformatics analyses to select 31 SNPs located in miRNA‐target binding sites using the miRNASNP v2.0, a solid database providing miRNA‐related SNPs for genetic research, and investigated their associations with risk of PC in two large case‐control studies totally including 1847 cases and 5713 controls. We observed that the SNP rs3802266 is significantly associated with increased risk of PC (odds ratio (OR) = 1.21, 95% confidence intervals (CI) = 1.11‐1.31, P = 1.29E‐05). Following luciferase reporter gene assays show that rs3802266‐G creates a stronger binding site for miR‐181a‐2‐3p in 3′ untranslated region (3′UTR) of the gene ZHX2. Expression quantitative trait loci (eQTL) analysis suggests that ZHX2 expression is lower in individuals carrying rs3802266‐G with increased PC risk. In conclusion, our findings highlight the involvement of miRNA‐binding SNPs in PC susceptibility and provide new clues for PC carcinogenesis.     

Testing hypotheses of population structuring in the Northeast Atlantic Ocean and Mediterranean Sea using the common cuttlefish Sepia officinalis

Population structuring in species inhabiting marine environments such as the Northeast Atlantic Ocean (NEA) and Mediterranean Sea (MS) has usually been explained based on past and present physical barriers to gene flow and isolation by distance (IBD). Here, we examined the relative importance of these factors on population structuring of the common cuttlefish Sepia officinalis by using methods of phylogenetic inference and hypothesis testing coupled with coalescent and classical population genetic parameter estimation. Individuals from 10 Atlantic and 15 Mediterranean sites were sequenced for 659 bp of the mitochondrial COI gene (259 sequences). IBD seems to be the main factor driving present and past genetic structuring of Sepia populations across the NEA-MS, both at large and small geographical scales. Such an evolutionary process agrees well with some of the biological features characterizing this cuttlefish species (short migrations, nektobenthic habit, benthic eggs hatching directly to benthic juveniles). Despite the many barriers to migration/gene flow suggested in the NEA-MS region, genetic population fragmentation due to past isolation of water masses (Pleistocene; 0.56 million years ago) and/or present-day oceanographic currents was only detected between the Aegean-Ionian and western Mediterranean Seas. Restricted gene flow associated with the Almería-Oran hydrographic front was also suggested between southern and eastern Spanish populations. Distinct population boundaries could not be clearly determined, except for the Aegean-Ionian stock. Two Atlantic and five Mediterranean samples showed evidence of current decline in genetic diversity, which may indicate over-exploitation of Sepia in both marine regions.