Reconstructing pedigrees: a stochastic perspective

A pedigree is a directed graph that describes how individuals are related through ancestry in a sexually-reproducing population. In this paper we explore the question of whether one can reconstruct a pedigree by just observing sequence data for present day individuals. This is motivated by the increasing availability of genomic sequences, but in this paper we take a more theoretical approach and consider what models of sequence evolution might allow pedigree reconstruction (given sufficiently long sequences). Our results complement recent work that showed that pedigree reconstruction may be fundamentally impossible if one uses just the degrees of relatedness between different extant individuals. We find that for certain stochastic processes, pedigrees can be recovered up to isomorphism from sufficiently long sequences.     

Reconstructing pedigrees: some identifiability questions for a recombination-mutation model

Pedigrees are directed acyclic graphs that represent ancestral relationships between individuals in a population. Based on a schematic recombination process, we describe two simple Markov models for sequences evolving on pedigrees—Model R (recombinations without mutations) and Model RM (recombinations with mutations). For these models, we ask an identifiability question: is it possible to construct a pedigree from the joint probability distribution of extant sequences? We present partial identifiability results for general pedigrees: we show that when the crossover probabilities are sufficiently small, certain spanning subgraph sequences can be counted from the joint distribution of extant sequences. We demonstrate how pedigrees that earlier seemed difficult to distinguish are distinguished by counting their spanning subgraph sequences.     

Reconstruction of pedigrees in clonal plant populations

We present a Bayesian method for the reconstruction of pedigrees in clonal populations using co-dominant genomic markers such as microsatellites and single nucleotide polymorphisms (SNPs). The accuracy of the algorithm is demonstrated for simulated data. We show that the joint estimation of parameters of interest such as the rate of self-fertilization is possible with high accuracy even with marker panels of moderate power. Classical methods can only assign a very limited number of statistically significant parentages in this case and would therefore fail. Statistical confidence is estimated by Markov Chain Monte Carlo (MCMC) sampling. The method is implemented in a fast and easy to use open source software that scales to large datasets with many thousand individuals.     

The power comparison of the haplotype-based collapsing tests and the variant-based collapsing tests for detecting rare variants in pedigrees

Background

Both common and rare genetic variants have been shown to contribute to the etiology of complex diseases. Recent genome-wide association studies (GWAS) have successfully investigated how common variants contribute to the genetic factors associated with common human diseases. However, understanding the impact of rare variants, which are abundant in the human population (one in every 17 bases), remains challenging. A number of statistical tests have been developed to analyze collapsed rare variants identified by association tests. Here, we propose a haplotype-based approach. This work inspired by an existing statistical framework of the pedigree disequilibrium test (PDT), which uses genetic data to assess the effects of variants in general pedigrees. We aim to compare the performance between the haplotype-based approach and the rare variant-based approach for detecting rare causal variants in pedigrees.

Results

Extensive simulations in the sequencing setting were carried out to evaluate and compare the haplotype-based approach with the rare variant methods that drew on a more conventional collapsing strategy. As assessed through a variety of scenarios, the haplotype-based pedigree tests had enhanced statistical power compared with the rare variants based pedigree tests when the disease of interest was mainly caused by rare haplotypes (with multiple rare alleles), and vice versa when disease was caused by rare variants acting independently. For most of other situations when disease was caused both by haplotypes with multiple rare alleles and by rare variants with similar effects, these two approaches provided similar power in testing for association.

Conclusions

The haplotype-based approach was designed to assess the role of rare and potentially causal haplotypes. The proposed rare variants-based pedigree tests were designed to assess the role of rare and potentially causal variants. This study clearly documented the situations under which either method performs better than the other. All tests have been implemented in a software, which was submitted to the Comprehensive R Archive Network (CRAN) for general use as a computer program named rvHPDT.     

Reconstructing the Indian origin and dispersal of the European Roma: a maternal genetic perspective

Previous genetic, anthropological and linguistic studies have shown that Roma (Gypsies) constitute a founder population dispersed throughout Europe whose origins might be traced to the Indian subcontinent. Linguistic and anthropological evidence point to Indo-Aryan ethnic groups from North-western India as the ancestral parental population of Roma. Recently, a strong genetic hint supporting this theory came from a study of a private mutation causing primary congenital glaucoma. In the present study, complete mitochondrial control sequences of Iberian Roma and previously published maternal lineages of other European Roma were analyzed in order to establish the genetic affinities among Roma groups, determine the degree of admixture with neighbouring populations, infer the migration routes followed since the first arrival to Europe, and survey the origin of Roma within the Indian subcontinent. Our results show that the maternal lineage composition in the Roma groups follows a pattern of different migration routes, with several founder effects, and low effective population sizes along their dispersal. Our data allowed the confirmation of a North/West migration route shared by Polish, Lithuanian and Iberian Roma. Additionally, eleven Roma founder lineages were identified and degrees of admixture with host populations were estimated. Finally, the comparison with an extensive database of Indian sequences allowed us to identify the Punjab state, in North-western India, as the putative ancestral homeland of the European Roma, in agreement with previous linguistic and anthropological studies.     

Synthetic biology,genome mining,and combinatorial biosynthesis of NRPS-derived antibiotics: a perspective

Combinatorial biosynthesis of novel secondary metabolites derived from nonribosomal peptide synthetases (NRPSs) has been in slow development for about a quarter of a century. Progress has been hampered by the complexity of the giant multimodular multienzymes. More recently, advances have been made on understanding the chemical and structural biology of these complex megaenzymes, and on learning the design rules for engineering functional hybrid enzymes. In this perspective, I address what has been learned about successful engineering of complex lipopeptides related to daptomycin, and discuss how synthetic biology and microbial genome mining can converge to broaden the scope and enhance the speed and robustness of combinatorial biosynthesis of NRPS-derived natural products for drug discovery.     

Likelihood approach for detecting imprinting and in utero maternal effects using general pedigrees from prospective family-based association studies

Genetic imprinting and in utero maternal effects are causes of parent-of-origin effect but they are confounded with each other. Tests attempting to detect only one of these effects would have a severely inflated type I error rate if the assumption of the absence of the other effect is violated. Some existing methods avoid the potential confounding by modeling imprinting and in utero maternal effect simultaneously. However, these methods are not amendable to extended families, which are commonly recruited in family-based studies. In this article, we propose a likelihood approach for detecting imprinting and maternal effects (LIME) using general pedigrees from prospective family-based association studies. LIME formulates the probability of familial genotypes without the Hardy-Weinberg equilibrium assumption by introducing a novel concept called conditional mating type between marry-in founders and their nonfounder spouses. Further, a logit link is used to model the penetrance. To deal with the issue of incomplete pedigree genotypic data, LIME imputes the unobserved genotypes implicitly by considering all compatible ones conditional on the observed genotypes. We carried out a simulation study to evaluate the relative power and type I error of LIME and two existing methods. The results show that the use of extended pedigree data, even with incomplete information, can achieve much greater power than using nuclear families for detecting imprinting and in utero maternal effects without leading to inflated type I error rates.     

The death factors: a combinatorial analysis

This paper presents theoretical and computational tools to understand how a small group of proteins, the death factors, are involved in widely different behavior of the cell. Experiments were done using a virtual laboratory that can simulate cellular response to different external stimuli. WARNING: It is not certain which of the theoretical protein clusters described here really occur in nature. In addition, the rules of cluster assembly are combinatorial, and thus an oversimplification to describe the real situation.     

HaploPainter: a tool for drawing pedigrees with complex haplotypes

SUMMARY: HaploPainter is a user-friendly pedigree-drawing application with special features for easy visualization of complex haplotype information. It has been developed to facilitate gene mapping in Mendelian diseases in terms of fast and reliable definition of the smallest critical interval harbouring the underlying gene defect. HaploPainter is written in Perl and may be used for visualization of haplotypes calculated by any of the common linkage programs. With special features like haplotype compression or the ability of marker section cut-out it particularly addresses the requirements for viewing large haplotypes as obtained by using for genome scans high-density marker panels of many thousands of single nucleotide polymorphisms (SNPs). AVAILABILITY: http://haplopainter.sourceforge.net/ CONTACT: holger.thiele@uni-koeln.de.     

Reconstructing ancestral character states: a critical reappraisal

Using parsimony to reconstruct ancestral character states on a phylogenetic tree has become a popular method for testing ecological and evolutionary hypotheses. Despite its popularity, the assumptions and uncertainties of reconstructing the ancestral states of a single character have received less attention than the much less challenging endeavor of reconstructing phylogenetic trees from many characters. Recent research suggests that parsimony reconstructions are often sensitive to violations of the almost universal assumption of equal probabilities of gains and losses. In addition, maximum likelihood has been developed as an alternative to parsimony reconstruction, and has also revealed a surprising amount of uncertainty in ancestral reconstructions.     

Wild pedigrees: the way forward

Metrics derived from pedigrees are key to investigating several major issues in evolutionary biology, including the quantitative genetic architecture of traits, inbreeding depression, and the evolution of cooperation and inbreeding avoidance. There is merit in studying these issues in natural populations experiencing spatially and temporally variable environmental conditions, since these analyses may yield different results from laboratory studies and allow us to understand population responses to rapid environmental change. Partial pedigrees are now available for several natural populations which are the subject of long-term individual-based studies, and analyses using these pedigrees are leading to important insights. Accurate pedigree construction supported by molecular genetic data is now feasible across a wide range of taxa, and even where only imprecise pedigrees are available it is possible to estimate the consequences of imprecision for the questions of interest. In outbred diploid populations, the pedigree approach is superior to analyses based on marker-based pairwise estimators of coancestry.     

Digestion-after-shuffling: a method for combinatorial library construction

We have developed a convenient method for family shuffling of amino acid sequences, termed digestion-after-shuffling. After DNA shuffling of homologous genes, plasmid mixture is extracted from a library and used for several double digestions with restriction enzymes. For each double digestion, two restriction enzymes are selected, corresponding to the single restriction sites of different parental genes. After digestions, fragments with expected sizes are obtained by gel purification and religated to construct recombinant plasmids. Thus, the obtained genes should be chimeras and have at least two restriction sites originating from different parental sequences.     

Haplotyping for disease association: a combinatorial approach

We consider a combinatorial problem derived from haplotyping a population with respect to a genetic disease, either recessive or dominant. Given a set of individuals, partitioned into healthy and diseased, and the corresponding sets of genotypes, we want to infer "bad' and "good' haplotypes to account for these genotypes and for the disease. Assume e.g. the disease is recessive. Then, the resolving haplotypes must consist of bad and good haplotypes, so that (i) each genotype belonging to a diseased individual is explained by a pair of bad haplotypes and (ii) each genotype belonging to a healthy individual is explained by a pair of haplotypes of which at least one is good. We prove that the associated decision problem is NP-complete. However, we also prove that there is a simple solution, provided the data satisfy a very weak requirement.     

A Markov chain Monte Carlo strategy for sampling from the joint posterior distribution of pedigrees and population parameters under a Fisher-Wright model with partial selfing

A simple population genetic model is presented for a hermaphrodite annual species, allowing both selfing and outcrossing. Those male gametes (pollen) responsible for outcrossing are assumed to disperse much further than seeds. Under this model, the pedigree of a sample from a single locality is loop-free. A novel Markov chain Monte Carlo strategy is presented for sampling from the joint posterior distribution of the pedigree of such a sample and the parameters of the population genetic model (including the selfing rate) given the genotypes of the sampled individuals at unlinked marker loci. The computational costs of this Markov chain Monte Carlo strategy scale well with the number of individuals in the sample, and the number of marker loci, but increase exponentially with the age (time since colonisation from the source population) of the local population. Consequently, this strategy is particularly suited to situations where the sample has been collected from a population which is the result of a recent colonisation process.     

海洋病毒荧光显微计数法的优化与应用

【目的】建立一种快速、稳定、可靠的海洋病毒计数方法。【方法】海水水样经福尔马林固定后,滤过孔径为0.02μm的Anodisc Al2O3膜。滤膜经SYBR Green I染色后,在相应波长的激发光下进行观察。借助荧光显微镜目镜网格尺,计数视野中的病毒颗粒,换算后获得样品中病毒的浓度。【结果】对具体实验方法进行了优化,可快速、稳定地对海水中的病毒计数。【结论】建立了一种适用于国内实验条件的、可靠的海洋病毒计数方法。     

Estimating genealogies from unlinked marker data: a Bayesian approach

An issue often encountered in statistical genetics is whether, or to what extent, it is possible to estimate the degree to which individuals sampled from a background population are related to each other, on the basis of the available genotype data and some information on the demography of the population. In this article, we consider this question using explicit modelling of the pedigrees and gene flows at unlinked marker loci, but then restricting ourselves to a relatively recent history of the population, that is, considering the genealogy at most some tens of generations backwards in time. As a computational tool we use a Markov chain Monte Carlo numerical integration on the state space of genealogies of the sampled individuals. As illustrations of the method, we consider the question of relatedness at the level of genes/genomes (IBD estimation), using both simulated and real data.     

Hepatoimmunology: a perspective

Premises for the subspecialty of hepatoimmunology include the recognition that the liver is a lymphoid organ with unique immunological properties. These properties ensure efficient innate defence against intestinal microbes and toxins, confer a particular capacity for induction of tolerance, and provide for apoptotic disposal of redundant lymphocytes. Pathological responses within the liver are elicited when: (i) hepatotropic viruses (hepatitis virus B and C) escape immune elimination and reside in hepatocytes; (ii) the liver becomes the site of autoimmune responses directed against either hepatocytes (autoimmune hepatitis) or biliary ductules (primary biliary cirrhosis); or (iii) the liver in the course of disposal of drugs generates neoantigens that provoke adverse allergic responses. Recent advances in the understanding of the immunopathogenesis of these entities are reviewed.     

Orthogonal combinatorial mutagenesis: a codon-level combinatorial mutagenesis method useful for low multiplicity and amino acid-scanning protocols

We describe here a method to generate combinatorial libraries of oligonucleotides mutated at the codon-level, with control of the mutagenesis rate so as to create predictable binomial distributions of mutants. The method allows enrichment of the libraries with single, double or larger multiplicity of amino acid replacements by appropriate choice of the mutagenesis rate, depending on the concentration of synthetic precursors. The method makes use of two sets of deoxynucleoside-phosphoramidites bearing orthogonal protecting groups [4,4′-dimethoxytrityl (DMT) and 9-fluorenylmethoxycarbonyl (Fmoc)] in the 5′ hydroxyl. These phosphoramidites are divergently combined during automated synthesis in such a way that wild-type codons are assembled with commercial DMT-deoxynucleoside-methyl-phosphoramidites while mutant codons are assembled with Fmoc-deoxynucleoside-methyl-phosphoramidites in an NNG/C fashion in a single synthesis column. This method is easily automated and suitable for low mutagenesis rates and large windows, such as those required for directed evolution and alanine scanning. Through the assembly of three oligonucleotide libraries at different mutagenesis rates, followed by cloning at the polylinker region of plasmid pUC18 and sequencing of 129 clones, we concluded that the method performs essentially as intended.