Approximating identity-by-descent matrices using multiple haplotype configurations on pedigrees

Identity-by-descent (IBD) matrix calculation is an important step in quantitative trait loci (QTL) analysis using variance component models. To calculate IBD matrices efficiently for large pedigrees with large numbers of loci, an approximation method based on the reconstruction of haplotype configurations for the pedigrees is proposed. The method uses a subset of haplotype configurations with high likelihoods identified by a haplotyping method. The new method is compared with a Markov chain Monte Carlo (MCMC) method (Loki) in terms of QTL mapping performance on simulated pedigrees. Both methods yield almost identical results for the estimation of QTL positions and variance parameters, while the new method is much more computationally efficient than the MCMC approach for large pedigrees and large numbers of loci. The proposed method is also compared with an exact method (Merlin) in small simulated pedigrees, where both methods produce nearly identical estimates of position-specific kinship coefficients. The new method can be used for fine mapping with joint linkage disequilibrium and linkage analysis, which improves the power and accuracy of QTL mapping.     

Ancestral modal Y-STR haplotype shared among Romani and South Indian populations

One of the primary unanswered questions regarding the dispersal of Romani populations concerns the geographical region and/or the Indian caste/tribe that gave rise to the proto-Romani group. To shed light on this matter, 161 Y-chromosomes from Roma, residing in two different provinces of Serbia, were analyzed. Our results indicate that the paternal gene pool of both groups is shaped by several strata, the most prominent of which, H1-M52, comprises almost half of each collection's patrilineages. The high frequency of M52 chromosomes in the two Roma populations examined may suggest that they descend from a single founder that has its origins in the Indian subcontinent. Moreover, when the Y-STR profiles of haplogroup H derived individuals in our Roma populations were compared to those typed in the South Indian emigrants from Malaysia and groups from Madras, Karnataka (Lingayat and Vokkaliga castes) and tribal Soligas, sharing of the two most common haplotypes was observed. These similarities suggest that South India may have been one of the contributors to the proto-Romanis. European genetic signatures (i.e., haplogroups E1b1b1a1b-V13, G2a-P15, I-M258, J2-M172 and R1-M173), on the other hand, were also detected in both groups, but at varying frequencies. The divergent European genetic signals in each collection are likely the result of differential gene flow and/or admixture with the European host populations but may also be attributed to dissimilar endogamous practices following the initial founder effect. Our data also support the notion that a number of haplogroups including G2a-P15, J2a3b-M67(xM92), I-M258 and E1b1b1-M35 were incorporated into the proto-Romani paternal lineages as migrants moved from northern India through Southwestern Asia, the Middle East and/or Anatolia into the Balkans.     

Ancestral consanguinity and mortality among three endogamous populations of Chittoor District, Andhra Pradesh, India

Consanguineous marriages have been practiced around the globe by many societies from time immemorial, particularly in South India. Consanguineous marriages play a major role in the health of a population, and diseases leading to mortality of the progeny are a consequence of detrimental recessive genes. To evaluate the effects of ancestral consanguinity on mortality in relation to consanguineous marriage, we have ascertained data from 1,500 women belonging to three endogamous communities (Akuthota Reddy, Odde, and Madiga) of Chittoor District, Andhra Pradesh, India. There were 500 women from each community. For each marriage we drew a family pedigree, extended upward to two earlier generations on either side of the spouses, to determine the prevalence and pattern of consanguinity, with detailed information on fertility and mortality. We observed a significant difference in the mortality rates between consanguineous and nonconsanguineous marriages when all the marriages of the women, women's parents, and (women's) husband's parents were considered in all three communities. In inbreeding, the offspring of earlier generations might have passed on deleterious genes to later generations (under unfavorable conditions), which resulted in a negative aspect of reproduction (among the offspring of the present couple).     

Frequency-based haplotype reconstruction from deep sequencing data of bacterial populations

Clonal populations accumulate mutations over time, resulting in different haplotypes. Deep sequencing of such a population in principle provides information to reconstruct these haplotypes and the frequency at which the haplotypes occur. However, this reconstruction is technically not trivial, especially not in clonal systems with a relatively low mutation frequency. The low number of segregating sites in those systems adds ambiguity to the haplotype phasing and thus obviates the reconstruction of genome-wide haplotypes based on sequence overlap information.Therefore, we present EVORhA, a haplotype reconstruction method that complements phasing information in the non-empty read overlap with the frequency estimations of inferred local haplotypes. As was shown with simulated data, as soon as read lengths and/or mutation rates become restrictive for state-of-the-art methods, the use of this additional frequency information allows EVORhA to still reliably reconstruct genome-wide haplotypes. On real data, we show the applicability of the method in reconstructing the population composition of evolved bacterial populations and in decomposing mixed bacterial infections from clinical samples.     

Full-length haplotype reconstruction to infer the structure of heterogeneous virus populations

Next-generation sequencing (NGS) technologies enable new insights into the diversity of virus populations within their hosts. Diversity estimation is currently restricted to single-nucleotide variants or to local fragments of no more than a few hundred nucleotides defined by the length of sequence reads. To study complex heterogeneous virus populations comprehensively, novel methods are required that allow for complete reconstruction of the individual viral haplotypes. Here, we show that assembly of whole viral genomes of ∼8600 nucleotides length is feasible from mixtures of heterogeneous HIV-1 strains derived from defined combinations of cloned virus strains and from clinical samples of an HIV-1 superinfected individual. Haplotype reconstruction was achieved using optimized experimental protocols and computational methods for amplification, sequencing and assembly. We comparatively assessed the performance of the three NGS platforms 454 Life Sciences/Roche, Illumina and Pacific Biosciences for this task. Our results prove and delineate the feasibility of NGS-based full-length viral haplotype reconstruction and provide new tools for studying evolution and pathogenesis of viruses.     

C3 types in some endogamous populations of Andhra Pradesh

The distribution of C3 phenotypes was studied in some endogamous caste groups of Andhra Pradesh, South India. The C3F allele was found to be at a low frequency and comparable to those found in castes of other regions on the Indian subcontinent.     

Acid phosphatase polymorphism in some endogamous populations of Andhra Pradesh

The distribution of acid phosphatase (ACP) phenotypes in six endogamous Brahmin sub-sects and in the Vysya, Mala and Madiga castes of Andhra Pradesh, India, were investigated. The ACP A gene frequency ranged from 0.167 to 0.254. The ACP C allele was observed in three Brahmin sub-sects.     

Esterase D polymorphism in some endogamous populations of Andhra Pradesh

Esterase D polymorphism was investigated in six endogamous Brahmin sub-sects and in the Mala and Madiga castes of Andrah Pradesh. The ESD 2 gene frequency ranged between 0.184 and 0.405. The frequencies obtained did not show variation in the range of ESD 2 gene frequency from the other populations of India.     

Ancestral sequence reconstruction for protein engineers

In addition to its value in the study of molecular evolution, ancestral sequence reconstruction (ASR) has emerged as a useful methodology for engineering proteins with enhanced properties. Proteins generated by ASR often exhibit unique or improved activity, stability, and/or promiscuity, all of which are properties that are valued by protein engineers. Comparison between extant proteins and evolutionary intermediates generated by ASR also allows protein engineers to identify substitutions that have contributed to functional innovation or diversification within protein families. As ASR becomes more widely adopted as a protein engineering approach, it is important to understand the applications, limitations, and recent developments of this technique. This review highlights recent exemplifications of ASR, as well as technical aspects of the reconstruction process that are relevant to protein engineering.     

The impact of using related individuals for haplotype reconstruction in population studies

Recent studies have highlighted the dangers of using haplotypes reconstructed directly from population data for a fine-scale mapping analysis. Family data may help resolve ambiguity, yet can be costly to obtain. This study is concerned with the following question: How much family data (if any) should be used to facilitate haplotype reconstruction in a population study? We conduct a simulation study to evaluate how changes in family information can affect the accuracy of haplotype frequency estimates and phase reconstruction. To reconstruct haplotypes, we introduce an EM-based algorithm that can efficiently accommodate unrelated individuals, parent-child trios, and arbitrarily large half-sib pedigrees. Simulations are conducted for a diverse set of haplotype frequency distributions, all of which have been previously published in empirical studies. A wide variety of important results regarding the effectiveness of using pedigree data in a population study are presented in a coherent, unified framework. Insight into the different properties of the haplotype frequency distribution that can influence experimental design is provided. We show that a preliminary estimate of the haplotype frequency distribution can be valuable in large population studies with fixed resources.     

Glyoxylase I phenotypes in some endogamous populations of Andhra Pradesh, India

The distribution of glyoxylase (GLO) I phenotypes in six endogamous subgroups of Brahmins and in the Mala and Madiga castes of Andhra Pradesh was investigated. The GLO I gene frequencies ranged from 0.2444 to 0.3575. The frequency of 0.3565 found in the Mala is the highest recorded on the Indian subcontinent.     

Inference of relationships in population data using identity-by-descent and identity-by-state

It is an assumption of large, population-based datasets that samples are annotated accurately whether they correspond to known relationships or unrelated individuals. These annotations are key for a broad range of genetics applications. While many methods are available to assess relatedness that involve estimates of identity-by-descent (IBD) and/or identity-by-state (IBS) allele-sharing proportions, we developed a novel approach that estimates IBD0, 1, and 2 based on observed IBS within windows. When combined with genome-wide IBS information, it provides an intuitive and practical graphical approach with the capacity to analyze datasets with thousands of samples without prior information about relatedness between individuals or haplotypes. We applied the method to a commonly used Human Variation Panel consisting of 400 nominally unrelated individuals. Surprisingly, we identified identical, parent-child, and full-sibling relationships and reconstructed pedigrees. In two instances non-sibling pairs of individuals in these pedigrees had unexpected IBD2 levels, as well as multiple regions of homozygosity, implying inbreeding. This combined method allowed us to distinguish related individuals from those having atypical heterozygosity rates and determine which individuals were outliers with respect to their designated population. Additionally, it becomes increasingly difficult to identify distant relatedness using genome-wide IBS methods alone. However, our IBD method further identified distant relatedness between individuals within populations, supported by the presence of megabase-scale regions lacking IBS0 across individual chromosomes. We benchmarked our approach against the hidden Markov model of a leading software package (PLINK), showing improved calling of distantly related individuals, and we validated it using a known pedigree from a clinical study. The application of this approach could improve genome-wide association, linkage, heterozygosity, and other population genomics studies that rely on SNP genotype data.     

PGM1 subtype polymorphism in 14 endogamous Dravidian-speaking populations of South India

Red cell hemolysates from 1,004 persons belonging to 14 population groups drawn from four South Indian states, Andhra Pradesh, Tamil Nadu, Karnataka, and Kerala, were tested for PGM1 subtypes. The groups are characterized by a high frequency of phenotype 1+1+ (range 36.98-71.64%) and the allele 1+ (range 60-79%). The groups exhibit marked heterogeneity for PGM1 locus. The results show a clear demarcation between tribes and Brahmin groups.     

Genetic differentiation among some endogamous populations of Andhra Pradesh,India

Endogamous populations belonging to Brahmin, Vysya, Mala, and Madiga castes of Andhra Pradesh, South India, were investigated for certain red cell enzyme and serum protein genetic markers. Frequency values were statistically analyzed to assess genetic variation among the populations. Average heterozygosity of ten loci and genetic diversity within and between the populations were calculated by using the methods of Nei. Nei's index was used to calculate genetic distances between the pairs of populations. A dendrogram was drawn adapting the modified unweighted pair group method suggested by Li, which agreed with the history of the populations.     

Ancestral traits, parental populations, and hybrids

The distance traveled over the centuries by any human population is a rocky road; the fusion of the paths of two meandering populations often leads to an even rougher one. Tracing such a wayward child of history back to its long-lost parents presents complex problems for the physical anthropologist. For those parents may share unequally in their legacy to their offspring at the start of the journey. Many genes may disappear into the gutters of genetic drift; some byways may stop at the barricades of selection; fresh migrants may join by side roads to swell the stream; and inbreeding may narrow the roadbed. Searching for the genes that the hybrid population holds in common with alleged parents provides one clue to the relative role of the ancestors. Comparison of measurements of the living hybrid with the nearest known relatives of the presumed parents yields another clue. Resuscitation of ancestors from their bones and teeth supplies yet another link. Distinctive features of such skeletal remains are most valuable where their inheritance is understood. Bones may yield another clue from substances left within them of known genetic control. And finally all the evidence of history, every glimpse caught of the population along the route, gives knowledge of the contribution of the forebears to the living people. The ideal estimate will utilize a wide variety of information about the known groups.     

Ancestral origins of the Machado-Joseph disease mutation: a worldwide haplotype study

Machado-Joseph disease (MJD) is an autosomal dominant neurodegenerative disorder originally described in families of Portuguese-Azorean ancestry. The cloning of the MJD1 gene allowed identification of the disease in many other populations, and MJD is now known to be the most common cause of dominant spinocerebellar ataxia. The hypothesis that its present world distribution could result from the spread of an original founder mutation has been raised, both at historical and molecular levels. In the present study, we tested this hypothesis by linkage-disequilibrium analysis of tightly linked polymorphisms and by haplotype comparison, in 249 families from different countries. We typed five microsatellite markers surrounding the MJD1 locus (D14S1015, D14S995, D14S973, D14S1016, and D14S977), and three intragenic single-base-pair polymorphisms (A(669)TG/G(669)TG, C(987)GG/G(987)GG, and TAA(1118)/TAC(1118)). The results show two different haplotypes, specific to the island of origin, in families of Azorean extraction. In families from mainland Portugal, both Azorean haplotypes can be found. The majority of the non-Portuguese families also share the same intragenic haplotype seen in the families coming from the island of Flores, but at least three other haplotypes were seen. These findings suggest two introductions of the mutation into the Portuguese population. Worldwide, the sharing of one intragenic haplotype by the majority of the families studied implies a founder mutation in MJD.     

Comparisons of methods for linkage analysis and haplotype reconstruction using extended pedigree data

We compare and contrast the performance of SIMPLE, a Monte Carlo based software, with that of several other methods for linkage and haplotype analyses, focusing on the simulated data from the New York City population. First, a whole-genome scan study based on the microsatellite markers was performed using GENEHUNTER. Because GENEHUNTER had to drop individuals for many of the pedigrees, we performed a follow-up study focusing on several regions of interest using SIMPLE, which can handle all pedigrees in their entirety. Second, 3 haplotyping programs, including that in SIMPLE, were used to reconstruct haplotypic configurations in pedigrees. SIMPLE emerges clearly as a preferred tool, as it can handle large pedigrees and produces haplotypic configurations without double recombinant haplotypes. For this study, we had knowledge of the simulating models at the time we performed the analysis.     

Precision of methods for calculating identity-by-descent matrices using multiple markers

A rapid, deterministic method (DET) based on a recursive algorithm and a stochastic method based on Markov Chain Monte Carlo (MCMC) for calculating identity-by-descent (IBD) matrices conditional on multiple markers were compared using stochastic simulation. Precision was measured by the mean squared error (MSE) of the relationship coefficients in predicting the true IBD relationships, relative to MSE obtained from using pedigree only. Comparisons were made when varying marker density, allele numbers, allele frequencies, and the size of full-sib families. The precision of DET was 75–99% relative to MCMC, but was not simply related to the informativeness of individual loci. For situations mimicking microsatellite markers or dense SNP, the precision of DET was ≥ 95% relative to MCMC. Relative precision declined for the SNP, but not microsatellites as marker density decreased. Full-sib family size did not affect the precision. The methods were tested in interval mapping and marker assisted selection, and the performance was very largely determined by the MSE. A multi-locus information index considering the type, number, and position of markers was developed to assess precision. It showed a marked empirical relationship with the observed precision for DET and MCMC and explained the complex relationship between relative precision and the informativeness of individual loci.     

Ancestral state reconstruction analysis of hymenopteran sex determination mechanisms

We provide the first phylogenetic evidence supporting complementary sex determination (CSD) as the ancestral mechanism for haplodiploidy in the Hymenoptera. It is currently not possible, however, to distinguish the evolutionary polarity of single locus (sl) CSD and multiple‐locus (ml) CSD given the available data. In this light, we discuss the seemingly maladaptive hypothesis of ml‐CSD ancestry, suggesting that collapse from ml‐CSD to sl‐CSD should remain a viable evolutionary hypothesis based on (i) likely weakening of frequency‐dependent selection on sex alleles under ml‐CSD and (ii) recent findings with respect to the evolutionary novelty of the complementary sex determiner gene in honeybees. Our findings help provide a phylogenetically informed blueprint for future sampling of sex determination mechanisms in the Hymenoptera, as they yield hypotheses for many unsampled or ambiguous taxa and highlight taxa whose further sampling will influence reconstruction of the evolutionary polarity of sex determination mechanisms in major clades.