Association of Genetic Variants with Self-Assessed Color Categories in Brazilians

The Brazilian population was formed by extensive admixture of three different ancestral roots: Amerindians, Europeans and Africans. Our previous work has shown that at an individual level, ancestry, as estimated using molecular markers, was a poor predictor of color in Brazilians. We now investigate if SNPs known to be associated with human skin pigmentation can be used to predict color in Brazilians. For that, we studied the association of fifteen SNPs, previously known to be linked with skin color, in 243 unrelated Brazilian individuals self-identified as White, Browns or Blacks from Rio de Janeiro and 212 unrelated Brazilian individuals self-identified as White or Blacks from São Paulo. The significance of association of SNP genotypes with self-assessed color was evaluated using partial regression analysis. After controlling for ancestry estimates as covariates, only four SNPs remained significantly associated with skin pigmentation: rs1426654 and rs2555364 within SLC24A5, rs16891982 at SLC45A2 and rs1042602 at TYR. These loci are known to be involved in melanin synthesis or transport of melanosomes. We found that neither genotypes of these SNPs, nor their combination with biogeographical ancestry in principal component analysis, could predict self-assessed color in Brazilians at an individual level. However, significant correlations did emerge at group level, demonstrating that even though elements other than skin, eye and hair pigmentation do influence self-assessed color in Brazilians, the sociological act of self-classification is still substantially dependent of genotype at these four SNPs.     

Total ancestry measure: quantifying the similarity in tree-like classification, with genomic applications

MOTIVATION: Many classifications of protein function such as Gene Ontology (GO) are organized in directed acyclic graph (DAG) structures. In these classifications, the proteins are terminal leaf nodes; the categories 'above' them are functional annotations at various levels of specialization and the computation of a numerical measure of relatedness between two arbitrary proteins is an important proteomics problem. Moreover, analogous problems are important in other contexts in large-scale information organization--e.g. the Wikipedia online encyclopedia and the Yahoo and DMOZ web page classification schemes. RESULTS: Here we develop a simple probabilistic approach for computing this relatedness quantity, which we call the total ancestry method. Our measure is based on counting the number of leaf nodes that share exactly the same set of 'higher up' category nodes in comparison to the total number of classified pairs (i.e. the chance for the same total ancestry). We show such a measure is associated with a power-law distribution, allowing for the quick assessment of the statistical significance of shared functional annotations. We formally compare it with other quantitative functional similarity measures (such as, shortest path within a DAG, lowest common ancestor shared and Azuaje's information-theoretic similarity) and provide concrete metrics to assess differences. Finally, we provide a practical implementation for our total ancestry measure for GO and the MIPS functional catalog and give two applications of it in specific functional genomics contexts. AVAILABILITY: The implementations and results are available through our supplementary website at: http://gersteinlab.org/proj/funcsim. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Demographic history shapes genomic ancestry in hybrid zones

Demographic factors such as migration rate and population size can impede or facilitate speciation. In hybrid zones, reproductive boundaries between species are tested and demography mediates the opportunity for admixture between lineages that are partially isolated. Genomic ancestry is a powerful tool for revealing the history of admixed populations, but models and methods based on local ancestry are rarely applied to structured hybrid zones. To understand the effects of demography on ancestry in hybrids zones, we performed individual‐based simulations under a stepping‐stone model, treating migration rate, deme size, and hybrid zone age as parameters. We find that the number of ancestry junctions (the transition points between genomic regions with different ancestries) and heterogenicity (the genomic proportion heterozygous for ancestry) are often closely connected to demographic history. Reducing deme size reduces junction number and heterogenicity. Elevating migration rate increases heterogenicity, but migration affects junction number in more complex ways. We highlight the junction frequency spectrum as a novel and informative summary of ancestry that responds to demographic history. A substantial proportion of junctions are expected to fix when migration is limited or deme size is small, changing the shape of the spectrum. Our findings suggest that genomic patterns of ancestry could be used to infer demographic history in hybrid zones.     

Signatures of hybridization and speciation in genomic patterns of ancestry*

Genomes sampled from hybrid zones between nascent species provide important clues into the speciation process. With advances in genome sequencing and single nucleotide polymorphism (SNP) genotyping, it is now feasible to measure variation in gene flow with high genomic resolution. This progress motivates the development of conceptual and analytical frameworks for hybrid zones that complement well‐established cline approaches. We extend the perspective that genomic distributions of ancestry are sensitive indicators of hybridization history. We use simulations to examine the behavior of the number of ancestry junctions—a simple summary of genomic patterns—in hybrid zones under increasingly realistic scenarios. Neutral simulations revealed that ancestry junction number is shaped by population structure, migration rate, and population size. Modeling multiple genetic architectures of hybrid dysfunction, with an emphasis on epistatic hybrid incompatibilities, showed that selection reduces junction number near loci that confer reproductive barriers. The magnitude of this signature was affected by the form of selection, dominance, and genomic location (autosome vs. sex chromosome) of incompatible loci. Our results suggest that researchers can identify loci involved in reproductive isolation by scanning hybrid genomes for local reductions in junction number. We outline necessary directions for future theory and method development to realize this goal.     

Genetic mapping and annotation of genomic microsatellites isolated from globe artichoke

Cynara cardunculus includes three taxa, the globe artichoke (subsp. scolymus L. Hegi), the cultivated cardoon (var. altilis) and their progenitor, the wild cardoon (var. sylvestris). Globe artichoke is an important component of the Mediterranean rural economy, but its improvement through breeding has been rather limited and its genome organization remains largely unexplored. Here, we report the isolation of 61 new microsatellite loci which amplified a total of 208 alleles in a panel of 22 C. cardunculus genotypes. Of these, 51 were informative for linkage analysis and 39 were used to increase marker density in the available globe artichoke genetic maps. Sequence analysis of the 22 loci associated with genes showed that 9 are located within coding sequence, with the repetitive domain probably being involved in DNA binding or in protein–protein interactions. The expression of the genes associated with 9 of the 22 microsatellite loci was demonstrated by RT-PCR.     

The breakdown of genomic ancestry blocks in hybrid lineages given a finite number of recombination sites

When a lineage originates from hybridization genomic blocks of contiguous ancestry from different ancestors are fragmented through genetic recombination. The resulting blocks are delineated by so called junctions, which accumulate with every generation that passes. Modeling the accumulation of ancestry block junctions can elucidate processes and timeframes of genomic admixture. Previous models have not addressed ancestry block dynamics for chromosomes that consist of a finite number of recombination sites. However, genomic data typically consist of informative markers that are interspersed with fragments for which no ancestry information is available. Hence, repeated recombination events may occur between markers, effectively removing existing junctions. Here, we present an analytical treatment of the dynamics of the mean number of junctions over time, taking into account the number of recombination sites per chromosome, population size, genetic map length, and the frequency of the ancestral species in the founding hybrid swarm. We describe the expected number of junctions using equidistant molecular markers and estimate the number of junctions using random markers. This extended theory of junctions thus reflects properties of empirical data and can serve to study the genomic patterns following admixture.     

Glottal fry and voice disguise: a case study in forensic phonetics

In recent legal proceedings, forensic phoneticians were called upon to analyse a tape-recorded message intended for the blackmail of a bank manager following the kidnap of his wife. The brief was to establish the likelihood that the tape recording may have been made by any one of three suspects, samples of whose speech were also made available. The comparison was greatly complicated by voice disguise employed by the speaker who recorded the kidnap tape. This disguise comprised a form of phonation described phonetically as ‘glottal fry’ or vocal ‘creak’. This form of phonation occurs normally in normal speech, but it has received most attention in relation to voice pathologies. On the other hand there are few references to its use as a form of voice disguise. This paper discusses the nature of the creak, and examines its effectiveness as voice disguise. In addition, a method is described for speaker identification regardless of the disguise. Results indicate that trained listeners without repeated presentations or instrumentation are able to match speakers with 65% accuracy when one voice is creaky, compared with 90% accuracy for undisguised voices. Using a Euclidean metric to compare the power spectra of the [s] sound, we find that creaky disguised voices may be correctly matched with the undisguised voice of the same speaker (9 distracters) in 5 cases out of 10. However, when the computer's task is made more similar to the perceptual task, selecting one speaker out of two, it achieves an accuracy of 81%. Implications for forensic phonetics are discussed.     

A comparative study of microsatellites among crocodiles and development of genomic resources for the critically endangered Indian gharial

Genetica - Next-generation sequencing has allowed us to explore new methods, where comparative and population genomics can be used simultaneously. Keeping this in mind, we surveyed and analyzed the...     

Drawbacks with the use of microsatellites in phylogeography: the song sparrow Melospiza melodia as a case study

Microsatellites are commonly used molecular markers in phylogeography, and many view them as superior to mitochondrial DNA (mtDNA) gene trees. Being based on frequencies of alleles, and not gene trees, microsatellites exhibit the same analytical drawbacks that resulted in the abandonment of allozymes in genetic studies of population history. I illustrate some these familiar drawbacks by reanalyzing microsatellite data on the song sparrow. Subspecies were previously evaluated with hierarchical analyses of molecular variance, suggesting that subspecies explain 8% of the total variance in microsatellite frequencies. However, this useful heuristic technique only evaluates a priori groupings, and the objective of the study ought to be to discover such groupings, not assume them. In fact, other arbitrary groupings of samples explained the same or greater amounts of variance, and I suggest that for testing subspecies limits, a gene tree is preferable. Grouping population samples by subspecies in the San Francisco Bay area accounts for 1.2% of the microsatellite variation, and despite claims that this informs conservation planning, the data do not support any particular population or subspecies as being genetically or evolutionarily significant. A distance phenogram was used to infer a sequential colonization of the Aleutian Islands, but because individuals were pooled into a priori groups and the phenogram was arbitrarily rooted, this conclusion is tenuous. A plot of heterozygosity vs number of alleles per sample showed that an equally parsimonious interpretation is that current genetic diversity tracks effective population size. Microsatellites should be replaced in nuclear-gene phylogeography by analyses of sequences, which will benefit the study of phylogeography, comparison of nuclear and mtDNA results, and aid in interpreting the results in a conservation context.     

Non-clonability correlates with genomic instability: a case study of a unique DNA region

Instability of eukaryotic DNA in constructs propagated in prokaryotic hosts is a frequently observed phenomenon. With the exception of a very high A+T-content and the presence of multiple repetitions, no general rule at the basis of this phenomenon is actually known. The intergenic spacer located between the pi and alpha(D) chicken alpha-type globin genes is frequently deleted from recombinant phages and plasmids. Here we have cloned this DNA fragment using a specially designed bacterial strain (SURE competent cells, Stratogene). Comparative analysis of DNA of recombinant clones bearing deletions and clones containing the intact genomic DNA fragment has revealed two important DNA sequence motifs that contribute to the unclonability of eukaryotic DNA in prokaryotic cells. First, the similarity to bacterial transposons (i.e. the presence of repeats flanking a several kilobase DNA fragment) may cause the loss of the fragment during propagation of the recombinant DNA in E. coli. Second, a high content of rotationally correlated kinkable elements (TG*CA steps) may result in non-clonability of the DNA sequence. Interestingly, the latter type of "unclonable" DNA sequence motifs identified in the globin gene domain is unstable (frequently rearranged) also in the eukaryotic chromosome resulting in a local polymorphism. In the chicken domain of alpha globin genes this unstable DNA sequence seems to be partially protected by interaction with nuclear matrix proteins.     

New polymorphic microsatellite loci for Theobroma cacao: isolation and characterization of microsatellites from enriched genomic libraries

Seventeen polymorphic microsatellite markers were isolated from enriched genomic libraries for Theobroma cacao, providing additional tools for studying the genetic diversity and map saturation of this species. These markers were characterized in 32 accessions of the T. cacao germplasm collection from the Centro de Pesquisas do Cacau. The number of alleles at each locus varied from 2 to 8, with an average of 4.41 alleles per locus. The polymorphism information content varied from 0.060 to 0.695, with an average of 0.333. The markers characterized in this study will be employed in map saturation studies and diversity assessments of cacao genotypes.     

AFLPs mark different genomic regions compared with RFLPs: a case study in tomato.

A set of three tomato chromosome 7 introgression lines (ILs) containing overlapping segments of Lycopersicon pennellii DNA was screened with a set of 10 EcoRI-MseI and 10 PstI-MseI AFLP primer combinations. A large number of markers were identified that mapped to one of the four regions of chromosome 7 defined by the set of three ILs. Because many of the identified markers have known map positions in three tomato reference maps, their location on the tomato genome could be verified. It was demonstrated that the three chromosome 7 ILs carried a chromosome 10 region harbouring a cluster of six AFLP markers that had not been detected before using RFLPs. The causes and implications of this observation are discussed.     

Population growth of human Y chromosomes: a study of Y chromosome microsatellites

We use variation at a set of eight human Y chromosome microsatellite loci to investigate the demographic history of the Y chromosome. Instead of assuming a population of constant size, as in most of the previous work on the Y chromosome, we consider a model which permits a period of recent population growth. We show that for most of the populations in our sample this model fits the data far better than a model with no growth. We estimate the demographic parameters of this model for each population and also the time to the most recent common ancestor. Since there is some uncertainty about the details of the microsatellite mutation process, we consider several plausible mutation schemes and estimate the variance in mutation size simultaneously with the demographic parameters of interest. Our finding of a recent common ancestor (probably in the last 120,000 years), coupled with a strong signal of demographic expansion in all populations, suggests either a recent human expansion from a small ancestral population, or natural selection acting on the Y chromosome.     

Parasite-specific immune response in adult Drosophila melanogaster: a genomic study

Insects of the order Diptera are vectors for parasitic diseases such as malaria, sleeping sickness and leishmania. In the search for genes encoding proteins involved in the antiparasitic response, we have used the protozoan parasite Octosporea muscaedomesticae for oral infections of adult Drosophila melanogaster. To identify parasite-specific response molecules, other flies were exposed to virus, bacteria or fungi in parallel. Analysis of gene expression patterns after 24 h of microbial challenge, using Affymetrix oligonucleotide microarrays, revealed a high degree of microbe specificity. Many serine proteases, key intermediates in the induction of insect immune responses, were uniquely expressed following infection of the different organisms. Several lysozyme genes were induced in response to Octosporea infection, while in other treatments they were not induced or downregulated. This suggests that lysozymes are important in antiparasitic defence.     

Cattle ancestry in bison: explanations for higher mtDNA than autosomal ancestry

Understanding and documenting the process of hybridization and introgression between related species is a major focus of recent evolutionary research using molecular techniques. Many North American bison herds have cattle ancestry introduced by crossbreeding over a century ago. Molecular estimates of this ancestry have shown much higher levels for cattle mtDNA than for autosomal cattle genes. A large part of this difference appears to be the result of partial reproductive isolation between the two species where only bison bull × domestic cow crosses are successful, and all the surviving progeny are females. In addition, selection against autosomal cattle genes in bison may have contributed to differential levels of cattle ancestry. The impact of selection against cattle mtDNA and gene flow of bison mtDNA are examined to explain particular combinations of mtDNA and autosomal cattle ancestry. A bottleneck, after the level of cattle ancestry in bison was reduced to a low level, is consistent with the high variance over autosomal loci observed for cattle ancestry, and differential selection among cattle loci in bison does not need to be invoked. Further examination of the cattle genome in bison may shed light on whether these markers, or their associated regions, are indeed neutral.     

Characterization and identification of (CT)n microsatellites in soybean using sheared genomic libraries.

Three small insert (300 to approximately 600 bp) sheared genomic libraries were constructed by pipetting and DNase I treatment of soybean DNA. About 15,000 clones from each library were screened for CT- simple sequence repeats (CT-SSRs). The CT-SSRs were abundant in the soybean genome at an estimated frequency of approximately one SSR per 110 kb of genomic DNA. Following the sequencing of 129 positive clones, the repeat types and frequency of CT repeats among the positive clones were characterized. Forty-nine primer pairs were designed and preliminarily evaluated for their ability to amplify genomic DNA from a set of six varieties, including parents of a mapping family. Amplified products were analyzed by 10% PAGE. Eighty-eight percent of the designed primers were able to amplify all these genomic DNAs using a single PCR profile of 53 degrees C annealing temperature, of which 22 (45%) were polymorphic in the six varieties, and 14 of them were polymorphic in the parents of the mapping family. The polymorphic primer sets were further assessed for allelic information using DNA from 16 soybean cultivars. The average number of alleles was 4, ranging from 2 to 7 with the highest polymorphism information content value 0.84. Fourteen of these SSRs were mapped, using an existing soybean RFLP map. The findings presented here will advance our understanding of the soybean genome, and assist in the mapping genome and discrimination of closely related varieties of this species.     

A simple and sensitive method for quantifying human genomic DNA in forensic specimen extracts

The analysis of DNA restriction fragment length polymorphisms by Southern blot hybridization requires that sufficient quantities of high molecular weight genomic DNA be extracted from biological specimens. Prior to analysis, it is necessary to determine the quantity and quality of the extracted DNA. For many applications, it is also desirable to determine the amount of DNA which is of human origin. In this report, we describe a simple and highly sensitive procedure for the specific quantification of human genomic DNA in forensic extracts or any biological sample. A small fraction of the extract is immobilized onto a nylon membrane and subsequently hybridized to p17H8 (D17Z1), a cloned probe which detects highly repetitive, primate-specific alpha satellite DNA. The procedure requires less than four hours to complete and can be used to quantify subnanogram amounts of hybridizable human genomic DNA.     

Allozymes,mtDNA and microsatellites study introgression in a stocked trout population in France

A comparative study of disturbed (downstream) and undisturbed (upstream) populations in a French Mediterranean river (the Sorgue) was conducted to assess the fate of introgressing hatchery-genes in the aboriginal gene pool. Variation was assessed at 27 protein-coding loci, three microsatellite loci and two mtDNA restriction sites. The results were compared to those from another Mediterranean drainage basin (that of the river Orb). The genetic pattern observed among markers and drainage basins was found to be noticeably different. A reduction in population size is hypothesized to explain the lower variability observed in the upper part of the river Sorgue. Estimations of hatchery gene flow were different between loci but in the same range of values between markers.     

Isolation, characterization and cross-species testing of microsatellites obtained from a sand smelt (Atherina boyeri) genomic library

This study reports the isolation and characterization of 11 polymorphic microsatellites from a sand smelt (Atherina boyeri) genomic library. Enrichment was performed with di-, tri- and tetranucleotide motifs following the FIASCO procedure (fast isolation by AFLP of sequences containing repeats). All loci were found to be in linkage and in Hardy-Weinberg equilibrium. This represents the first microsatellite isolation for the family Atherinidae and the isolated loci were accordingly tested on four additional species of the family: two recognized (A. presbyter and A. hepsetus) and two proposed ('punctata' and 'non-punctata' forms). Moreover their cross-species suitability on Menidia menidia, belonging to the same order but to the family Atherinopsidae, was also tested.