Ancestral genomes, sex, and the population structure of Trypanosoma cruzi

Acquisition of detailed knowledge of the structure and evolution of Trypanosoma cruzi populations is essential for control of Chagas disease. We profiled 75 strains of the parasite with five nuclear microsatellite loci, 24Salpha RNA genes, and sequence polymorphisms in the mitochondrial cytochrome oxidase subunit II gene. We also used sequences available in GenBank for the mitochondrial genes cytochrome B and NADH dehydrogenase subunit 1. A multidimensional scaling plot (MDS) based in microsatellite data divided the parasites into four clusters corresponding to T. cruzi I (MDS-cluster A), T. cruzi II (MDS-cluster C), a third group of T. cruzi strains (MDS-cluster B), and hybrid strains (MDS-cluster BH). The first two clusters matched respectively mitochondrial clades A and C, while the other two belonged to mitochondrial clade B. The 24Salpha rDNA and microsatellite profiling data were combined into multilocus genotypes that were analyzed by the haplotype reconstruction program PHASE. We identified 141 haplotypes that were clearly distributed into three haplogroups (X, Y, and Z). All strains belonging to T. cruzi I (MDS-cluster A) were Z/Z, the T. cruzi II strains (MDS-cluster C) were Y/Y, and those belonging to MDS-cluster B (unclassified T. cruzi) had X/X haplogroup genotypes. The strains grouped in the MDS-cluster BH were X/Y, confirming their hybrid character. Based on these results we propose the following minimal scenario for T. cruzi evolution. In a distant past there were at a minimum three ancestral lineages that we may call, respectively, T. cruzi I, T. cruzi II, and T. cruzi III. At least two hybridization events involving T. cruzi II and T. cruzi III produced evolutionarily viable progeny. In both events, the mitochondrial recipient (as identified by the mitochondrial clade of the hybrid strains) was T. cruzi II and the mitochondrial donor was T. cruzi III.     

Unique features of HLA-mediated HIV evolution in a Mexican cohort: a comparative study

Background

Mounting evidence indicates that HLA-mediated HIV evolution follows highly stereotypic pathways that result in HLA-associated footprints in HIV at the population level. However, it is not known whether characteristic HLA frequency distributions in different populations have resulted in additional unique footprints.

Methods

The phylogenetic dependency network model was applied to assess HLA-mediated evolution in datasets of HIV pol sequences from free plasma viruses and peripheral blood mononuclear cell (PBMC)-integrated proviruses in an immunogenetically unique cohort of Mexican individuals. Our data were compared with data from the IHAC cohort, a large multi-center cohort of individuals from Canada, Australia and the USA.

Results

Forty three different HLA-HIV codon associations representing 30 HLA-HIV codon pairs were observed in the Mexican cohort (q < 0.2). Strikingly, 23 (53%) of these associations differed from those observed in the well-powered IHAC cohort, strongly suggesting the existence of unique characteristics in HLA-mediated HIV evolution in the Mexican cohort. Furthermore, 17 of the 23 novel associations involved HLA alleles whose frequencies were not significantly different from those in IHAC, suggesting that their detection was not due to increased statistical power but to differences in patterns of epitope targeting. Interestingly, the consensus differed in four positions between the two cohorts and three of these positions could be explained by HLA-associated selection. Additionally, different HLA-HIV codon associations were seen when comparing HLA-mediated selection in plasma viruses and PBMC archived proviruses at the population level, with a significantly lower number of associations in the proviral dataset.

Conclusion

Our data support universal HLA-mediated HIV evolution at the population level, resulting in detectable HLA-associated footprints in the circulating virus. However, it also strongly suggests that unique genetic backgrounds in different HIV-infected populations may influence HIV evolution in a particular direction as particular HLA-HIV codon associations are determined by specific HLA frequency distributions. Our analysis also suggests a dynamic HLA-associated evolution in HIV with fewer HLA-HIV codon associations observed in the proviral compartment, which is likely enriched in early archived HIV sequences, compared to the plasma virus compartment. These results highlight the importance of comparative HIV evolutionary studies in immunologically different populations worldwide.     

Immuno-informatics: Mining genomes for vaccine components

The complete genome sequences of more than 60 microbes have been completed in the past decade. Concurrently, a series of new informatics tools, designed to harness this new wealth of information, have been developed. Some of these new tools allow researchers to select regions of microbial genomes that trigger immune responses. These regions, termed epitopes, are ideal components of vaccines. When the new tools are used to search for epitopes, this search is usually coupled with in vitro screening methods; an approach that has been termed computational immunology or immuno-informatics.Researchers are now implementing these combined methods to scan genomic sequences for vaccine components. They are thereby expanding the number of different proteins that can be screened for vaccine development, while narrowing this search to those regions of the proteins that are extremely likely to induce an immune response.As the tools improve, it may soon be feasible to skip over many of the in vitro screening steps, moving directly from genome sequence to vaccine design. The present article reviews the work of several groups engaged in the development of immuno-informatics tools and illustrates the application of these tools to the process of vaccine discovery.     

Simple sequences are ubiquitous repetitive components of eukaryotic genomes

Simple sequences are stretches of DNA which consist of only one, or a few tandemly repeated nucleotides, for example poly (dA) X poly (dT) or poly (dG-dT) X poly (dC-dA). These two types of simple sequence have been shown to be repetitive and interspersed in many eukaryotic genomes. Several other types have been found by sequencing eukaryotic DNA. In this report we have undertaken a systematical survey for simple sequences. We hybridized synthetical simple sequence DNA to genome blots of phylogenetically different organisms. We found that many, probably even all possible types of simple sequence are repetitive components of eukaryotic genomes. We propose therefore that they arise by common mechanisms namely slippage replication and unequal crossover and that they might have no general function with regards to gene expression. This latter inference is supported by the fact that we have detected simple sequences only in the metabolically inactive micronucleus of the protozoan Stylonychia, but not in the metabolically active macronucleus which is derived from the micronucleus by chromosome diminution.     

Association of serum lipid components and obesity with genetic ancestry in an admixed population of elderly women

The prevalence of metabolic disorders varies among ethnic populations and these disorders represent a critical health care issue for elderly women. This study investigated the correlation between genetic ancestry and body composition, metabolic traits and clinical status in a sample of elderly women. Clinical, nutritional and anthropometric data were collected from 176 volunteers. Genetic ancestry was estimated using 23 ancestry-informative markers. Pearsons correlation test was used to examine the relationship between continuous variables and an independent samples t-test was used to compare the means of continuous traits within categorical variables. Overall ancestry was a combination of European (57.49%), Native American (25.78%) and African (16.73%). Significant correlations were found for European ancestry with body mass index (r = 0.165; p = 0.037) and obesity (mean difference (MD) = 5.3%; p = 0.042). African ancestry showed a significant correlation with LDL (r = 0.159, p = 0.035), VLDL (r = −0.185; p = 0.014), hypertriglyceridemia (MD = 6.4%; p = 0.003) and hyperlipidemia (MD = 4.8%; p = 0.026). Amerindian ancestry showed a significant correlation with triglyceride levels (r = 0.150; p = 0.047) and hypertriglyceridemia (MD = 4.5%; p = 0.039). These findings suggest that genetic admixture may influence the etiology of lipid metabolism-related diseases and obesity in elderly women.     

Estimating kinship in admixed populations

Genome-wide association studies (GWASs) are commonly used for the mapping of genetic loci that influence complex traits. A problem that is often encountered in both population-based and family-based GWASs is that of identifying cryptic relatedness and population stratification because it is well known that failure to appropriately account for both pedigree and population structure can lead to spurious association. A number of methods have been proposed for identifying relatives in samples from homogeneous populations. A strong assumption of population homogeneity, however, is often untenable, and many GWASs include samples from structured populations. Here, we consider the problem of estimating relatedness in structured populations with admixed ancestry. We propose a method, REAP (relatedness estimation in admixed populations), for robust estimation of identity by descent (IBD)-sharing probabilities and kinship coefficients in admixed populations. REAP appropriately accounts for population structure and ancestry-related assortative mating by using individual-specific allele frequencies at SNPs that are calculated on the basis of ancestry derived from whole-genome analysis. In simulation studies with related individuals and admixture from highly divergent populations, we demonstrate that REAP gives accurate IBD-sharing probabilities and kinship coefficients. We apply REAP to the Mexican Americans in Los Angeles, California (MXL) population sample of release 3 of phase III of the International Haplotype Map Project; in this sample, we identify third- and fourth-degree relatives who have not previously been reported. We also apply REAP to the African American and Hispanic samples from the Women's Health Initiative SNP Health Association Resource (WHI-SHARe) study, in which hundreds of pairs of cryptically related individuals have been identified.     

Non-LTR retrotransposons (LINEs) as ubiquitous components of plant genomes

During the course of work aimed at isolating a rice gene from Oryza australiensis by PCR, the oligonucleotide primers used were found to generate a fragment that showed sequence homology to the endonuclease (EN) region of the maize non-LTR retrotransposon (LINE) Cin4. We carried out further PCRs using oligonucleotide primers that hybridized to these sequences, and found that they amplified several fragments, each with homology to the EN regions, from Oryza sativa cv. Nipponbare as well as O. australiensis. We mapped the approximate locations of two rice LINE homologues by screening clones in a YAC library made from a rice (O. sativa) genome, and found that each homologue was present in a low copy number apparently at nonspecific regions on rice chromosomes. We then carried out PCR using degenerate oligonucleotide primers which hybridized to the rice LINE homologues and Cin4 to ascertain whether LINE homologues are present in a variety of members of the plant kingdom, including angiosperms, gymnosperms, bracken, horsetail and liverwort. Cloning and nucleotide sequencing revealed that 53 clones obtained from 27 out of 33 plant species contained LINE homologues. In addition to these homologues, we identified four homologues with EN regions in the Arabidopsis thaliana genome by a computer search of databases. The nucleotide sequences of almost all the LINE homologues were greatly diverged, but the derived amino acid sequences were well conserved, and all contained glutamic acid and tyrosine residues at almost the same relative positions as in the the active site regions of AP (apurinic/apyrimidinic)-endonucleases. The EN regions in the LINE homologues from closely related plant species show a closer phylogenetic relationship, indicating that sequence divergence during vertical transmission has been a major influence upon the evolution of plant LINEs. Received: 13 July 1998 / Accepted: 13 October 1998     

Mapping of disease-associated variants in admixed populations

Recent developments in high-throughput genotyping and whole-genome sequencing will enhance the identification of disease loci in admixed populations. We discuss how a more refined estimation of ancestry benefits both admixture mapping and association mapping, making disease loci identification in admixed populations more powerful.     

Estimating local ancestry in admixed populations

Large-scale genotyping of SNPs has shown a great promise in identifying markers that could be linked to diseases. One of the major obstacles involved in performing these studies is that the underlying population substructure could produce spurious associations. Population substructure can be caused by the presence of two distinct subpopulations or a single pool of admixed individuals. In this work, we focus on the latter, which is significantly harder to detect in practice. New advances in this research direction are expected to play a key role in identifying loci that are different among different populations and are still associated with a disease. We evaluated current methods for inference of population substructure in such cases and show that they might be quite inaccurate even in relatively simple scenarios. We therefore introduce a new method, LAMP (Local Ancestry in adMixed Populations), which infers the ancestry of each individual at every single-nucleotide polymorphism (SNP). LAMP computes the ancestry structure for overlapping windows of contiguous SNPs and combines the results with a majority vote. Our empirical results show that LAMP is significantly more accurate and more efficient than existing methods for inferring locus-specific ancestries, enabling it to handle large-scale datasets. We further show that LAMP can be used to estimate the individual admixture of each individual. Our experimental evaluation indicates that this extension yields a considerably more accurate estimate of individual admixture than state-of-the-art methods such as STRUCTURE or EIGENSTRAT, which are frequently used for the correction of population stratification in association studies.     

Ancestral TSH mechanism signals summer in a photoperiodic mammal

In mammals, day-length-sensitive (photoperiodic) seasonal breeding cycles depend on the pineal hormone melatonin, which modulates secretion of reproductive hormones by the anterior pituitary gland [1]. It is thought that melatonin acts in the hypothalamus to control reproduction through the release of neurosecretory signals into the pituitary portal blood supply, where they act on pituitary endocrine cells [2]. Contrastingly, we show here that during the reproductive response of Soay sheep exposed to summer day lengths, the reverse applies: Melatonin acts directly on anterior-pituitary cells, and these then relay the photoperiodic message back into the hypothalamus to control neuroendocrine output. The switch to long days causes melatonin-responsive cells in the pars tuberalis (PT) of the anterior pituitary to increase production of thyrotrophin (TSH). This acts locally on TSH-receptor-expressing cells in the adjacent mediobasal hypothalamus, leading to increased expression of type II thyroid hormone deiodinase (DIO2). DIO2 initiates the summer response by increasing hypothalamic tri-iodothyronine (T3) levels. These data and recent findings in quail [3] indicate that the TSH-expressing cells of the PT play an ancestral role in seasonal reproductive control in vertebrates. In mammals this provides the missing link between the pineal melatonin signal and thyroid-dependent seasonal biology.     

Reconstructing genetic ancestry blocks in admixed individuals

A chromosome in an individual of recently admixed ancestry resembles a mosaic of chromosomal segments, or ancestry blocks, each derived from a particular ancestral population. We consider the problem of inferring ancestry along the chromosomes in an admixed individual and thereby delineating the ancestry blocks. Using a simple population model, we infer gene-flow history in each individual. Compared with existing methods, which are based on a hidden Markov model, the Markov-hidden Markov model (MHMM) we propose has the advantage of accounting for the background linkage disequilibrium (LD) that exists in ancestral populations. When there are more than two ancestral groups, we allow each ancestral population to admix at a different time in history. We use simulations to illustrate the accuracy of the inferred ancestry as well as the importance of modeling the background LD; not accounting for background LD between markers may mislead us to false inferences about mixed ancestry in an indigenous population. The MHMM makes it possible to identify genomic blocks of a particular ancestry by use of any high-density single-nucleotide-polymorphism panel. One application of our method is to perform admixture mapping without genotyping special ancestry-informative-marker panels.     

Linkage analysis of a complex disease through use of admixed populations

Linkage disequilibrium arising from the recent admixture of genetically distinct populations can be potentially useful in mapping genes for complex diseases. McKeigue has proposed a method that conditions on parental admixture to detect linkage. We show that this method tests for linkage only under specific assumptions, such as equal admixture in the parental generation and admixture that occurs in a single generation. In practice, these assumptions are unlikely to hold for natural populations, resulting in an inflation of the type I error rate when testing for linkage by this method. In this article, we generalize McKeigue's approach of testing for linkage to allow two different admixture models: (1) intermixture admixture and (2) continuous gene flow. We calculate the sample size required for a genomewide search by this method under different disease models: multiplicative, additive, recessive, and dominant. Our results show that the sample size required to obtain 90% power to detect a putative mutant allele at a genomewide significance level of 5% can usually be achieved in practice if informative markers are available at a density of 2 cM.     

Match probabilities in racially admixed populations.

The calculation of match probabilities is the most contentious issue dividing prosecution and defense experts in the forensic applications of DNA fingerprinting. In particular, defense experts question the applicability of the population genetic laws of Hardy-Weinberg and linkage equilibrium to racially admixed American populations. Linkage equilibrium justifies the product rule for computing match probabilities across loci. The present paper suggests a method of bounding match probabilities that depends on modeling gene descent from ancestral populations to contemporary populations under the assumptions of Hardy-Weinberg and linkage equilibrium only in the ancestral populations. Although these bounds are conservative from the defendant's perspective, they should be small enough in practice to satisfy prosecutors.     

Identification and functional characterization of gene components of Type VI Secretion system in bacterial genomes

A new secretion system, called the Type VI Secretion system (T6SS), was recently reported in Vibrio cholerae, Pseudomonas aeruginosa and Burkholderia mallei. A total of 18 genes have been identified to be belonging to this secretion system in V. cholerae. Here we attempt to identify presence of T6SS in other bacterial genomes. This includes identification of orthologous sequences, conserved motifs, domains, families, 3D folds, genomic islands containing T6SS components, phylogenetic profiles and protein-protein association of these components. Our analysis indicates presence of T6SS in 42 bacteria and its absence in most of their non-pathogenic species, suggesting the role of T6SS in imparting pathogenicity to an organism. Analysis of genomic regions containing T6SS components, phylogenetic profiles and protein-protein association of T6SS components indicate few additional genes which could be involved in this secretion system. Based on our studies, functional annotations were assigned to most of the components. Except one of the genes, we could group all the other genes of T6SS into those belonging to the puncturing device, and those located in the outer membrane, transmembrane and inner membrane. Based on our analysis, we have proposed a model of T6SS and have compared the same with the other bacterial secretion systems.     

Ancestral Mechanisms in Modern Environments

It is commonly assumed that the desire for a thin female physique and its pathological expression in eating disorders result from a social pressure for thinness. However, such widespread behavior may be better understood not merely as the result of arbitrary social pressure, but as an exaggerated expression of behavior that may have once been adaptive. The reproductive suppression hypothesis suggests that natural selection shaped a mechanism for adjusting female reproduction to socioecological conditions by altering the amount of body fat. In modern Western culture, social and ecological cues, which would have signaled the need for temporary postponement of reproduction in ancestral environments, may now be experienced to an unprecedented intensity and duration.     

Combinations of Ancestral Modules in Proteins

Twenty-seven protein sequence elements, six to nine amino acids long, were extracted from 15 phylogenetically diverse complete prokaryotic proteomes. The elements are present in all of these proteomes, with at least one copy each (omnipresent elements), and have presumably been conserved since the last universal common ancestor (LUCA). All these omnipresent elements are identified in crystallized protein structures as parts of highly conserved closed loops, 25–30 residues long, thus representing the closed-loop modules discovered in 2000 by Berezovsky et al. The omnipresent peptides make up seven distinct groups, of which the largest groups, Aleph and Beth, contain 18 and four elements, respectively, which are related but different, while five other groups are represented by only one element each. The LUCA modules appear with one or several copies per protein molecule in a variety of combinations depending on the functional identity of the corresponding protein. The functional involvement of individual LUCA modules is outlined on the basis of known protein annotations. Analyses of all the related sequences in a large, formatted protein sequence space suggest that many, if not all, of the 27 omnipresent elements have a common sequence origin. This sequence space network analysis may lead to elucidation of the earliest stages of protein evolution.     

Case-control studies of association in structured or admixed populations

Case-control tests for association are an important tool for mapping complex-trait genes. But population structure can invalidate this approach, leading to apparent associations at markers that are unlinked to disease loci. Family-based tests of association can avoid this problem, but such studies are often more expensive and in some cases--particularly for late-onset diseases--are impractical. In this review article we describe a series of approaches published over the past 2 years which use multilocus genotype data to enable valid case-control tests of association, even in the presence of population structure. These tests can be classified into two categories. "Genomic control" methods use the independent marker loci to adjust the distribution of a standard test statistic, while "structured association" methods infer the details of population structure en route to testing for association. We discuss the statistical issues involved in the different approaches and present results from simulations comparing the relative performance of the methods under a range of models.