遗传算法和单体型组装加权最小字符翻转问题

单体型组装加权最小字符翻转（WMLF）问题指定个体联配的加权DNA片断数据,翻转权值和最小的SNP位点以推测出该个体的一对单体型。该问题是NP-难的,至今尚无实用的搜索寻优算法。根据DNA测序片段数据的特点提出了一种遗传算法。对于实际的生物实验数据,即使数据很大,该算法也可以在较短的时间得到WMLF问题的满意解,具有良好的可扩展性和较高的实用价值。     

A Long-Read Sequencing Approach for Direct Haplotype Phasing in Clinical Settings

The reconstruction of individual haplotypes can facilitate the interpretation of disease risks; however, high costs and technical challenges still hinder their assessment in clinical settings. Second-generation sequencing is the gold standard for variant discovery but, due to the production of short reads covering small genomic regions, allows only indirect haplotyping based on statistical methods. In contrast, third-generation methods such as the nanopore sequencing platform developed by Oxford Nanopore Technologies (ONT) generate long reads that can be used for direct haplotyping, with fewer drawbacks. However, robust standards for variant phasing in ONT-based target resequencing efforts are not yet available. In this study, we presented a streamlined proof-of-concept workflow for variant calling and phasing based on ONT data in a clinically relevant 12-kb region of the APOE locus, a hotspot for variants and haplotypes associated with aging-related diseases and longevity. Starting with sequencing data from simple amplicons of the target locus, we demonstrated that ONT data allow for reliable single-nucleotide variant (SNV) calling and phasing from as little as 60 reads, although the recognition of indels is less efficient. Even so, we identified the best combination of ONT read sets (600) and software (BWA/Minimap2 and HapCUT2) that enables full haplotype reconstruction when both SNVs and indels have been identified previously using a highly-accurate sequencing platform. In conclusion, we established a rapid and inexpensive workflow for variant phasing based on ONT long reads. This allowed for the analysis of multiple samples in parallel and can easily be implemented in routine clinical practice, including diagnostic testing.     

Implementation of Bayesian methods to identify SNP and haplotype regions with transmission ratio distortion across the whole genome: TRDscan v.1.0

Realized deviations from the expected Mendelian inheritance of alleles from heterozygous parents have been previously reported in a broad range of organisms (i.e., transmission ratio distortion; TRD). Various biological mechanisms affecting gametes, embryos, fetuses, or even postnatal offspring can produce patterns of TRD. However, knowledge about its prevalence and potential causes in livestock species is still scarce. Specific Bayesian models have been recently developed for the analyses of TRD for biallelic loci, which accommodated a wide range of population structures, enabling TRD investigation in livestock populations. The parameterization of these models is flexible and allows the study of overall (parent-unspecific) TRD and sire- and dam-specific TRD. This research aimed at deriving Bayesian models for fitting TRD on the basis of haplotypes, testing the models for both haplotype- and SNP-based methods in simulated data and actual Holstein genotypes, and developing a specific software for TRD analyses. Results obtained on simulated data sets showed that the statistical power of the analysis increased with sample size of trios (n), proportion of heterozygous parents, and the magnitude of the TRD. On the other hand, the statistical power to detect TRD decreased with the number of alleles at each loci. Bayesian analyses showed a strong Pearson correlation coefficient (≥0.97) between simulated and estimated TRD that reached the significance level of Bayes factor ≥10 for both single-marker and haplotype analyses when n ≥ 25. Moreover, the accuracy in terms of the mean absolute error decreased with the increase of the sample size and increased with the number of alleles at each loci. Using real data (55,732 genotypes of Holstein trios), SNP- and haplotype-based distortions were detected with overall TRD, sire-TRD, or dam-TRD, showing different magnitudes of TRD and statistical relevance. Additionally, the haplotype-based method showed more ability to capture TRD compared with individual SNP. To discard possible random TRD in real data, an approximate empirical null distribution of TRD was developed. The program TRDscan v.1.0 was written in Fortran 2008 language and provides a powerful statistical tool to scan for TRD regions across the whole genome. This developed program is freely available at http://www.casellas.info/files/TRDscan.zip.     

Short communication: new alleles of the bovine kappa-casein gene revealed by resequencing and haplotype inference analysis

We tested the hypothesis that extensive undiscovered genetic diversity exists in important functional genes from domestic and wild cattle species (Bos spp.). We resequenced 483 bp of a key exon (exon IV) from the kappa (κ)-casein gene (CSN3) for a panel of samples of domestic cattle from 8 countries and a close relative species, the gayal (Bos frontalis). Six single nucleotide polymorphisms were identified. Haplotype inference revealed 12 haplotypes, of which 8 were newly discovered. Among these 8 new haplotypes, 5 differed by one nonsynonymous mutation and 3 differed by one silent mutation from previously well-characterized CSN3 alleles. From those, one was shared by the gayal and Zebu, was different from CSN3*B at position Ile136Thr, and showed a close phylogenetic relationship with the banteng, gaur, and yak. The other 7 new haplotypes were detected in our panel of worldwide local cattle breeds but were absent from previously reported commercial breeds. These results support the hypothesis that genetic diversity at the coding region of CSN3 has been underestimated. This study also highlights how important it is to resequence functionally important genes in worldwide local cattle breeds, many of which are threatened by extinction or replacement by commercial breeds.     

Short communication: The combined use of linkage disequilibrium–based haploblocks and allele frequency–based haplotype selection methods enhances genomic evaluation accuracy in dairy cattle

The construction and use of haploblocks [adjacent single nucleotide polymorphisms (SNP) in strong linkage disequilibrium] for genomic evaluation is advantageous, because the number of effects to be estimated can be reduced without discarding relevant genomic information. Furthermore, haplotypes (the combination of 2 or more SNP) can increase the probability of capturing the quantitative trait loci effect compared with individual SNP markers. With regards to haplotypes, the allele frequency parameter is also of interest, because as a selection criterion, it allows the number of rare alleles to be reduced, and the effects of those alleles are usually difficult to estimate. We have proposed a simple pipeline that simultaneously incorporates linkage disequilibrium and allele frequency information in genomic evaluation, and here we present the first results obtained with this procedure. We used a population of 2,235 progeny-tested bulls from the Montbéliarde breed for the tests. Phenotype data were available in the form of daughter yield deviations on 5 production traits, and genotype data were available from the 50K SNP chip. We conducted a classical validation study by splitting the population into training (80% oldest animals) and validation (20% youngest animals) sets to emulate a real-life scenario in which the selection candidates had no available phenotype data. We measured all reported parameters for the validation set. Our results proved that the proposed method was indeed advantageous, and that the accuracy of genomic evaluation could be improved. Compared with results from a genomic BLUP analysis, correlations between daughter yield deviations (a proxy for true) and genomic estimated breeding values increased by an average of 2.7 percentage points for the 5 traits. Inflation of the genomic evaluation of the selection candidates was also significantly reduced. The proposed method outperformed the other SNP and haplotype-based tests we had evaluated in a previous study. The combination of linkage disequilibrium–based haploblocks and allele frequency–based haplotype selection methods is a promising way to improve the efficiency of genomic evaluation. Further work is needed to optimize each step in the proposed analysis pipeline.     

Antioxidant activity of yogurt made from milk characterized by different casein haplotypes and fortified with chestnut and sulla honeys

The aim of this work was to evaluate the antioxidant activity of yogurt made from milk characterized by different casein (CN) haplotypes (α_s1-, β-, κ-CN) and fortified with chestnut and sulla honeys. The CN haplotype was determined by isoelectric focusing, whereas antioxidant activity of yogurt was measured using 2,2’-azino-bis-3-ethylbenzothiazoline-6-sulfonic acid and ferric-reducing antioxidant power. The statistical analysis showed a significant effect of the studied factors. The results showed that chestnut honey presented the highest phenolic acid and flavonoid contents, which are closely associated with its high antioxidant activity. The antioxidant activity of fortified yogurt samples was affected both by different CN haplotypes and by type of honey added. Yogurts fortified with chestnut honey showed higher antioxidant activity than those fortified with sulla honey. The different behavior observed among the fortified yogurts led us to hypothesize that the effects of protein-polyphenol complex on antioxidant activity are interactive. The results suggest that milk proteins polymorphism and polyphenols play different roles in affecting the bioavailability and the antioxidant activity of yogurt.     

Variants of the Low Oxygen Sensors EGLN1 and HIF-1AN Associated with Acute Mountain Sickness

Two low oxygen sensors, Egl nine homolog 1 (EGLN1) and hypoxia-inducible factor 1-α inhibitor (HIF-1AN), play pivotal roles in the regulation of HIF-1α, and high altitude adaption may be involved in the pathology of acute mountain sickness (AMS). Here, we aimed to analyze single nucleotide polymorphisms (SNPs) in the untranslated regions of the EGLN1 and HIF-1AN genes and SNPs chosen from a genome-wide adaptation study of the Han Chinese population. To assess the association between EGLN1 and HIF-1AN SNPs and AMS in a Han Chinese population, a case–control study was performed including 190 patients and 190 controls. In total, thirteen SNPs were genotyped using the MassARRAY^® MALDI-TOF system. Multiple genetic models were tested; The Akaike’s information criterion (AIC) and Bayesian information criterion (BIC) values indicated that the dominant model may serve as the best-fit model for rs12406290 and rs2153364 of significant difference. However, these data were not significant after Bonferroni correction. No significant association was noted between AMS and rs12757362, rs1339894, rs1361384, rs2009873, rs2739513 or rs2486729 before and after Bonferroni correction. Further haplotype analyses indicated the presence of two blocks in EGLN1; one block consists of rs12406290-rs2153364, located upstream of the EGLN1 gene. Carriers of the “GG” haplotype of rs12406290-rs2153364 exhibited an increased risk of AMS after adjustments for age and smoking status. However, no significant association was observed among HIF-1AN 3''-untranslated region (3''-UTR) polymorphisms, haplotype and AMS. Our study indicates that variants in the EGLN1 5''-UTR influence the susceptibility to AMS in a Han Chinese population.     

Polymorphisms in the LPL and CETP Genes and Haplotype in the ESR1 Gene Are Associated with Metabolic Syndrome in Women from Southwestern Mexico

Metabolic syndrome (MetS) is a combination of metabolic disorders associated with an increased risk for cardiovascular disease (CVD). Studies in women reported associations between polymorphisms in ESR1, LPL and CETP genes and MetS. Our aim was to evaluate the association between variants in ESR1, LPL and CETP genes with MetS and its components. Four hundred and eighty women were analyzed, anthropometric features and biochemical profiles were evaluated, and genotyping was performed by real-time PCR. We found an association with elevated glucose levels (odds ratio (OR) = 2.9; p = 0.013) in carrying the AA genotype of rs1884051 in the ESR1 gene compared with the GG genotype, and the CC genotype of rs328 in the LPL gene was associated with MetS compared to the CG or GG genotype (OR = 2.8; p = 0.04). Moreover, the GA genotype of rs708272 in the CETP gene is associated with MetS compared to the GG or AA genotype (OR = 1.8; p = 0.006). In addition the ACTCCG haplotype in the ESR1 gene is associated with a decrease in the risk of MetS (OR = 0.02; p < 0.001). In conclusion, our results show the involvement of the variants of ESR1, LPL and CETP genes in metabolic events related to MetS or some of its features.     

三元家庭基因数据的单体分型和单体型频率估计

研究了在门德尔遗传定理和哈代-维恩伯格平衡假设下,三元家庭基因型数据的单体分型和单体型频率估计问题.过去的研究仅仅关注个体间没有联系或者含有一般家系信息的基因型数据,而对这种特殊的三元家庭关注得不够.考虑到HAPMAP数据库中有一部分数据就基于这种三元家庭,现在有越来越多的需求要求直接分析这种特殊的家系结构.提出一个两段式的三元家庭中单体型频率的估计方法:i) 分型阶段,找出每一个三元家庭零重组单体构型;ii) 频率估计阶段,在前一阶段得到的单体构型基础上,应用EM算法来估计单体型频率.在程序包TRIOHAP中用C语言实现了单体分型算法和EM算法,并且使用模拟和实际数据测试了TRIOHAP的有效性和效率.实验结果表明,TRIOHAP要比其他那些忽略了三元家庭信息的常见单体型频率估计软件运行快很多.进一步地,由于TRIOHAP利用了这些信息,其估计结果更加可靠.     

一种用带权SNP片断重建单体型的启发式算法

近年来,单体型检测问题已经得到了广泛的研究,成为计算生物学最热门的领域之一.本文对个体单体型重建问题进行研究,提出一种基于带权最少字符修改模型重建单体型的启发式算法HAW.HAW算法首先生成一对初始单体型,然后通过对初始单体型的不断扩充完成重建.实验结果表明HAW算法能有效求解模型,得到较以往算法更高的重建率,且算法运行速度较快,具有很高的实用价值.