DnaSP, DNA polymorphism analyses by the coalescent and other methods

SUMMARY

DnaSP is a software package for the analysis of DNA polymorphism data. Present version introduces several new modules and features which, among other options allow: (1) handling big data sets (approximately 5 Mb per sequence); (2) conducting a large number of coalescent-based tests by Monte Carlo computer simulations; (3) extensive analyses of the genetic differentiation and gene flow among populations; (4) analysing the evolutionary pattern of preferred and unpreferred codons; (5) generating graphical outputs for an easy visualization of results.

AVAILABILITY

The software package, including complete documentation and examples, is freely available to academic users from: http://www.ub.es/dnasp

The map-based sequence of the rice genome

Rice, one of the world's most important food plants, has important syntenic relationships with the other cereal species and is a model plant for the grasses. Here we present a map-based, finished quality sequence that covers 95% of the 389 Mb genome, including virtually all of the euchromatin and two complete centromeres. A total of 37,544 non-transposable-element-related protein-coding genes were identified, of which 71% had a putative homologue in Arabidopsis. In a reciprocal analysis, 90% of the Arabidopsis proteins had a putative homologue in the predicted rice proteome. Twenty-nine per cent of the 37,544 predicted genes appear in clustered gene families. The number and classes of transposable elements found in the rice genome are consistent with the expansion of syntenic regions in the maize and sorghum genomes. We find evidence for widespread and recurrent gene transfer from the organelles to the nuclear chromosomes. The map-based sequence has proven useful for the identification of genes underlying agronomic traits. The additional single-nucleotide polymorphisms and simple sequence repeats identified in our study should accelerate improvements in rice production.

Role of genotype in the cycle of violence in maltreated children

We studied a large sample of male children from birth to adulthood to determine why some children who are maltreated grow up to develop antisocial behavior, whereas others do not. A functional polymorphism in the gene encoding the neurotransmitter-metabolizing enzyme monoamine oxidase A (MAOA) was found to moderate the effect of maltreatment. Maltreated children with a genotype conferring high levels of MAOA expression were less likely to develop antisocial problems. These findings may partly explain why not all victims of maltreatment grow up to victimize others, and they provide epidemiological evidence that genotypes can moderate children's sensitivity to environmental insults.

To err (meiotically) is human: the genesis of human aneuploidy

Aneuploidy (trisomy or monosomy) is the most commonly identified chromosome abnormality in humans, occurring in at least 5% of all clinically recognized pregnancies. Most aneuploid conceptuses perish in utero, which makes this the leading genetic cause of pregnancy loss. However, some aneuploid fetuses survive to term and, as a class, aneuploidy is the most common known cause of mental retardation. Despite the devastating clinical consequences of aneuploidy, relatively little is known of how trisomy and monosomy originate in humans. However, recent molecular and cytogenetic approaches are now beginning to shed light on the non-disjunctional processes that lead to aneuploidy.

Initial sequence of the chimpanzee genome and comparison with the human genome

Here we present a draft genome sequence of the common chimpanzee (Pan troglodytes). Through comparison with the human genome, we have generated a largely complete catalogue of the genetic differences that have accumulated since the human and chimpanzee species diverged from our common ancestor, constituting approximately thirty-five million single-nucleotide changes, five million insertion/deletion events, and various chromosomal rearrangements. We use this catalogue to explore the magnitude and regional variation of mutational forces shaping these two genomes, and the strength of positive and negative selection acting on their genes. In particular, we find that the patterns of evolution in human and chimpanzee protein-coding genes are highly correlated and dominated by the fixation of neutral and slightly deleterious alleles. We also use the chimpanzee genome as an outgroup to investigate human population genetics and identify signatures of selective sweeps in recent human evolution.

5-HTTLPR polymorphism impacts human cingulate-amygdala interactions: a genetic susceptibility mechanism for depression

Carriers of the short allele of a functional 5' promoter polymorphism of the serotonin transporter gene have increased anxiety-related temperamental traits, increased amygdala reactivity and elevated risk of depression. Here, we used multimodal neuroimaging in a large sample of healthy human subjects to elucidate neural mechanisms underlying this complex genetic association. Morphometrical analyses showed reduced gray matter volume in short-allele carriers in limbic regions critical for processing of negative emotion, particularly perigenual cingulate and amygdala. Functional analysis of those regions during perceptual processing of fearful stimuli demonstrated tight coupling as a feedback circuit implicated in the extinction of negative affect. Short-allele carriers showed relative uncoupling of this circuit. Furthermore, the magnitude of coupling inversely predicted almost 30% of variation in temperamental anxiety. These genotype-related alterations in anatomy and function of an amygdala-cingulate feedback circuit critical for emotion regulation implicate a developmental, systems-level mechanism underlying normal emotional reactivity and genetic susceptibility for depression.

Human Gene Mutation Database (HGMD): 2003 update

The Human Gene Mutation Database (HGMD) constitutes a comprehensive core collection of data on germ-line mutations in nuclear genes underlying or associated with human inherited disease (www.hgmd.org). Data catalogued includes: single base-pair substitutions in coding, regulatory and splicing-relevant regions; micro-deletions and micro-insertions; indels; triplet repeat expansions as well as gross deletions; insertions; duplications; and complex rearrangements. Each mutation is entered into HGMD only once in order to avoid confusion between recurrent and identical-by-descent lesions. By March 2003, the database contained in excess of 39,415 different lesions detected in 1,516 different nuclear genes, with new entries currently accumulating at a rate exceeding 5,000 per annum. Since its inception, HGMD has been expanded to include cDNA reference sequences for more than 87% of listed genes, splice junction sequences, disease-associated and functional polymorphisms, as well as links to data present in publicly available online locus-specific mutation databases. Although HGMD has recently entered into a licensing agreement with Celera Genomics (Rockville, MD), mutation data will continue to be made freely available via the Internet.

Multifactor-dimensionality reduction reveals high-order interactions among estrogen-metabolism genes in sporadic breast cancer

One of the greatest challenges facing human geneticists is the identification and characterization of susceptibility genes for common complex multifactorial human diseases. This challenge is partly due to the limitations of parametric-statistical methods for detection of gene effects that are dependent solely or partially on interactions with other genes and with environmental exposures. We introduce multifactor-dimensionality reduction (MDR) as a method for reducing the dimensionality of multilocus information, to improve the identification of polymorphism combinations associated with disease risk. The MDR method is nonparametric (i.e., no hypothesis about the value of a statistical parameter is made), is model-free (i.e., it assumes no particular inheritance model), and is directly applicable to case-control and discordant-sib-pair studies. Using simulated case-control data, we demonstrate that MDR has reasonable power to identify interactions among two or more loci in relatively small samples. When it was applied to a sporadic breast cancer case-control data set, in the absence of any statistically significant independent main effects, MDR identified a statistically significant high-order interaction among four polymorphisms from three different estrogen-metabolism genes. To our knowledge, this is the first report of a four-locus interaction associated with a common complex multifactorial disease.

Replication validity of genetic association studies

The rapid growth of human genetics creates countless opportunities for studies of disease association. Given the number of potentially identifiable genetic markers and the multitude of clinical outcomes to which these may be linked, the testing and validation of statistical hypotheses in genetic epidemiology is a task of unprecedented scale. Meta-analysis provides a quantitative approach for combining the results of various studies on the same topic, and for estimating and explaining their diversity. Here, we have evaluated by meta-analysis 370 studies addressing 36 genetic associations for various outcomes of disease. We show that significant between-study heterogeneity (diversity) is frequent, and that the results of the first study correlate only modestly with subsequent research on the same association. The first study often suggests a stronger genetic effect than is found by subsequent studies. Both bias and genuine population diversity might explain why early association studies tend to overestimate the disease protection or predisposition conferred by a genetic polymorphism. We conclude that a systematic meta-analytic approach may assist in estimating population-wide effects of genetic risk factors in human disease.

Molecular basis of human hypertension: role of angiotensinogen

Essential hypertension is a common human disease believed to result from the interplay of multiple genetic and environmental determinants. In genetic studies of two large panels of hypertensive sibships from widely separated geographical areas, we obtained evidence of genetic linkage between the angiotensinogen gene (AGT) and hypertension, demonstrated association of AGT molecular variants with the disease, and found significant differences in plasma concentrations of angiotensinogen among hypertensive subjects with different AGT genotypes. The corroboration and replication afforded by these results support the interpretation that molecular variants of AGT constitute inherited predispositions to essential hypertension in humans.

Hitchhiking under positive Darwinian selection

Positive selection can be inferred from its effect on linked neutral variation. In the restrictive case when there is no recombination, all linked variation is removed. If recombination is present but rare, both deterministic and stochastic models of positive selection show that linked variation hitchhikes to either low or high frequencies. While the frequency distribution of variation can be influenced by a number of evolutionary processes, an excess of derived variants at high frequency is a unique pattern produced by hitchhiking (derived refers to the nonancestral state as determined from an outgroup). We adopt a statistic, H, to measure an excess of high compared to intermediate frequency variants. Only a few high-frequency variants are needed to detect hitchhiking since not many are expected under neutrality. This is of particular utility in regions of low recombination where there is not much variation and in regions of normal or high recombination, where the hitchhiking effect can be limited to a small (<1 kb) region. Application of the H test to published surveys of Drosophila variation reveals an excess of high frequency variants that are likely to have been influenced by positive selection.

The amphioxus genome and the evolution of the chordate karyotype

Lancelets ('amphioxus') are the modern survivors of an ancient chordate lineage, with a fossil record dating back to the Cambrian period. Here we describe the structure and gene content of the highly polymorphic approximately 520-megabase genome of the Florida lancelet Branchiostoma floridae, and analyse it in the context of chordate evolution. Whole-genome comparisons illuminate the murky relationships among the three chordate groups (tunicates, lancelets and vertebrates), and allow not only reconstruction of the gene complement of the last common chordate ancestor but also partial reconstruction of its genomic organization, as well as a description of two genome-wide duplications and subsequent reorganizations in the vertebrate lineage. These genome-scale events shaped the vertebrate genome and provided additional genetic variation for exploitation during vertebrate evolution.

Genetics of gene expression surveyed in maize, mouse and man

Treating messenger RNA transcript abundances as quantitative traits and mapping gene expression quantitative trait loci for these traits has been pursued in gene-specific ways. Transcript abundances often serve as a surrogate for classical quantitative traits in that the levels of expression are significantly correlated with the classical traits across members of a segregating population. The correlation structure between transcript abundances and classical traits has been used to identify susceptibility loci for complex diseases such as diabetes and allergic asthma. One study recently completed the first comprehensive dissection of transcriptional regulation in budding yeast, giving a detailed glimpse of a genome-wide survey of the genetics of gene expression. Unlike classical quantitative traits, which often represent gross clinical measurements that may be far removed from the biological processes giving rise to them, the genetic linkages associated with transcript abundance affords a closer look at cellular biochemical processes. Here we describe comprehensive genetic screens of mouse, plant and human transcriptomes by considering gene expression values as quantitative traits. We identify a gene expression pattern strongly associated with obesity in a murine cross, and observe two distinct obesity subtypes. Furthermore, we find that these obesity subtypes are under the control of different loci.

An amino-acid taste receptor

The sense of taste provides animals with valuable information about the nature and quality of food. Mammals can recognize and respond to a diverse repertoire of chemical entities, including sugars, salts, acids and a wide range of toxic substances. Several amino acids taste sweet or delicious (umami) to humans, and are attractive to rodents and other animals. This is noteworthy because L-amino acids function as the building blocks of proteins, as biosynthetic precursors of many biologically relevant small molecules, and as metabolic fuel. Thus, having a taste pathway dedicated to their detection probably had significant evolutionary implications. Here we identify and characterize a mammalian amino-acid taste receptor. This receptor, T1R1+3, is a heteromer of the taste-specific T1R1 and T1R3 G-protein-coupled receptors. We demonstrate that T1R1 and T1R3 combine to function as a broadly tuned L-amino-acid sensor responding to most of the 20 standard amino acids, but not to their D-enantiomers or other compounds. We also show that sequence differences in T1R receptors within and between species (human and mouse) can significantly influence the selectivity and specificity of taste responses.

A second genetic polymorphism in methylenetetrahydrofolate reductase (MTHFR) associated with decreased enzyme activity

A common mutation in methylenetetrahydrofolate reductase (MTHFR), C677T, results in a thermolabile variant with reduced activity. Homozygous mutant individuals (approximately 10% of North Americans) are predisposed to mild hyperhomocysteinemia, when their folate status is low. This genetic-nutrient interactive effect is believed to increase the risk for neural tube defects and vascular disease. In this communication, we characterize a second common variant in MTHFR (A1298C), an E to A substitution. Homozygosity was observed in approximately 10% of Canadian individuals. This polymorphism was associated with decreased enzyme activity; homozygotes had approximately 60% of control activity in lymphocytes. Heterozygotes for both the C677T and the A1298C mutation, approximately 15% of individuals, had 50-60% of control activity, a value that was lower than that seen in single heterozygotes for the C677T variant. No individuals were homozygous for both mutations. Additional studies of the A1298C mutation, in the absence and presence of the C677T mutation, are warranted, to adequately address the role of this new genetic variant in complex traits. A silent genetic variant, T1317C, was identified in the same exon. It was relatively infrequent (allele frequency 5%) in our study group, but was quite common in a small sample of African individuals (allele frequency 39%).

Multifactor dimensionality reduction software for detecting gene-gene and gene-environment interactions

MOTIVATION

Polymorphisms in human genes are being described in remarkable numbers. Determining which polymorphisms and which environmental factors are associated with common, complex diseases has become a daunting task. This is partly because the effect of any single genetic variation will likely be dependent on other genetic variations (gene-gene interaction or epistasis) and environmental factors (gene-environment interaction). Detecting and characterizing interactions among multiple factors is both a statistical and a computational challenge. To address this problem, we have developed a multifactor dimensionality reduction (MDR) method for collapsing high-dimensional genetic data into a single dimension thus permitting interactions to be detected in relatively small sample sizes. In this paper, we describe the MDR approach and an MDR software package.

RESULTS

We developed a program that integrates MDR with a cross-validation strategy for estimating the classification and prediction error of multifactor models. The software can be used to analyze interactions among 2-15 genetic and/or environmental factors. The dataset may contain up to 500 total variables and a maximum of 4000 study subjects.

AVAILABILITY

Information on obtaining the executable code, example data, example analysis, and documentation is available upon request.

SUPPLEMENTARY INFORMATION

All supplementary information can be found at http://phg.mc.vanderbilt.edu/Software/MDR.

Massive migration from the steppe was a source for Indo-European languages in Europe

We generated genome-wide data from 69 Europeans who lived between 8,000-3,000 years ago by enriching ancient DNA libraries for a target set of almost 400,000 polymorphisms. Enrichment of these positions decreases the sequencing required for genome-wide ancient DNA analysis by a median of around 250-fold, allowing us to study an order of magnitude more individuals than previous studies and to obtain new insights about the past. We show that the populations of Western and Far Eastern Europe followed opposite trajectories between 8,000-5,000 years ago. At the beginning of the Neolithic period in Europe, ∼8,000-7,000 years ago, closely related groups of early farmers appeared in Germany, Hungary and Spain, different from indigenous hunter-gatherers, whereas Russia was inhabited by a distinctive population of hunter-gatherers with high affinity to a ∼24,000-year-old Siberian. By ∼6,000-5,000 years ago, farmers throughout much of Europe had more hunter-gatherer ancestry than their predecessors, but in Russia, the Yamnaya steppe herders of this time were descended not only from the preceding eastern European hunter-gatherers, but also from a population of Near Eastern ancestry. Western and Eastern Europe came into contact ∼4,500 years ago, as the Late Neolithic Corded Ware people from Germany traced ∼75% of their ancestry to the Yamnaya, documenting a massive migration into the heartland of Europe from its eastern periphery. This steppe ancestry persisted in all sampled central Europeans until at least ∼3,000 years ago, and is ubiquitous in present-day Europeans. These results provide support for a steppe origin of at least some of the Indo-European languages of Europe.

Linkage disequilibrium in humans: models and data

In this review, we describe recent empirical and theoretical work on the extent of linkage disequilibrium (LD) in the human genome, comparing the predictions of simple population-genetic models to available data. Several studies report significant LD over distances longer than those predicted by standard models, whereas some data from short, intergenic regions show less LD than would be expected. The apparent discrepancies between theory and data present a challenge-both to modelers and to human geneticists-to identify which important features are missing from our understanding of the biological processes that give rise to LD. Salient features may include demographic complications such as recent admixture, as well as genetic factors such as local variation in recombination rates, gene conversion, and the potential segregation of inversions. We also outline some implications that the emerging patterns of LD have for association-mapping strategies. In particular, we discuss what marker densities might be necessary for genomewide association scans.

Are rare variants responsible for susceptibility to complex diseases?

Little is known about the nature of genetic variation underlying complex diseases in humans. One popular view proposes that mapping efforts should focus on identification of susceptibility mutations that are relatively old and at high frequency. It is generally assumed-at least for modeling purposes-that selection against complex disease mutations is so weak that it can be ignored. In this article, I propose an explicit model for the evolution of complex disease loci, incorporating mutation, random genetic drift, and the possibility of purifying selection against susceptibility mutations. I show that, for the most plausible range of mutation rates, neutral susceptibility alleles are unlikely to be at intermediate frequencies and contribute little to the overall genetic variance for the disease. Instead, it seems likely that the bulk of genetic variance underlying diseases is due to loci where susceptibility mutations are mildly deleterious and where there is a high overall mutation rate to the susceptible class. At such loci, the total frequency of susceptibility mutations may be quite high, but there is likely to be extensive allelic heterogeneity at many of these loci. I discuss some practical implications of these results for gene mapping efforts.

Interleukin-10: new perspectives on an old cytokine

Interleukin-10 (IL-10) has long been recognized to have potent and broad-spectrum anti-inflammatory activity, which has been unequivocally established in various models of infection, inflammation, and even in cancer. However, because of the marginal successes of the initial clinical trials using recombinant IL-10, some of the interest in this cytokine as an anti-inflammatory therapeutic has diminished. New work showing IL-10 production from regulatory T cells and even T-helper 1 T cells has reinvigorated the field and revealed the power of this cytokine to influence immune responses. Furthermore, new preclinical studies suggest that combination therapies, using antibodies to IL-10 along with chemotherapy, can be effective in treating bacterial, viral, or neoplastic diseases. Studies to understand IL-10 gene expression in the various cell types may lead to new therapeutics to enhance or inhibit IL-10 production. In this review, we summarize what is known about the regulation of IL-10 gene expression by various immune cells. We speculate on the promise that this cytokine holds to influence immune responses and mitigate immune pathologies.

Genetic polymorphisms of cytochrome P450 2D6 (CYP2D6): clinical consequences, evolutionary aspects and functional diversity

CYP2D6 is of great importance for the metabolism of clinically used drugs and about 20-25% of those are metabolised by this enzyme. In addition, the enzyme utilises hydroxytryptamines as endogenous substrates. The polymorphism of the enzyme results in poor, intermediate, efficient or ultrarapid metabolisers (UMs) of CYP2D6 drugs. It is plausible that the UM genotype, where more than one active gene on one allele occurs, is the outcome of selective dietary selection in certain populations in North East Africa. The UM phenotype affects 5.5% of the population in Western Europe. A hypothesis for the evolutionary basis behind selection for CYP2D6 gene duplications is presented in relation to selection for Cyp6 variants in insecticide resistant Drosophila strains. The polymorphism of CYP2D6 significantly affects the pharmacokinetics of about 50% of the drugs in clinical use, which are CYP2D6 substrates. The consequences of the polymorphism at ordinary drug doses can be either adverse drug reactions or no drug response. Examples are presented where CYP2D6 polymorphism affects the efficacy and costs of drug treatment. Predictive CYP2D6 genotyping is estimated by the author to be beneficial for treatment of about 30-40% of CYP2D6 drug substrates, that is, for about 7-10% of all drugs clinically used, although prospective clinical studies are necessary to evaluate the exact benefit of drug selection and dosage based on the CYP2D6 genotype.

Moderation of the effect of adolescent-onset cannabis use on adult psychosis by a functional polymorphism in the catechol-O-methyltransferase gene: longitudinal evidence of a gene X environment interaction

BACKGROUND

Recent evidence documents that cannabis use by young people is a modest statistical risk factor for psychotic symptoms in adulthood, such as hallucinations and delusions, as well as clinically significant schizophrenia. The vast majority of cannabis users do not develop psychosis, however, prompting us to hypothesize that some people are genetically vulnerable to the deleterious effects of cannabis.

METHODS

In a longitudinal study of a representative birth cohort followed to adulthood, we tested why cannabis use is associated with the emergence of psychosis in a minority of users, but not in others.

RESULTS

A functional polymorphism in the catechol-O-methyltransferase (COMT) gene moderated the influence of adolescent cannabis use on developing adult psychosis. Carriers of the COMT valine158 allele were most likely to exhibit psychotic symptoms and to develop schizophreniform disorder if they used cannabis. Cannabis use had no such adverse influence on individuals with two copies of the methionine allele.

CONCLUSIONS

These findings provide evidence of a gene x environment interaction and suggest that a role of some susceptibility genes is to influence vulnerability to environmental pathogens.

A functional polymorphism in the monoamine oxidase A gene promoter

We describe a new polymorphism upstream of the gene for monoamine oxidase A (MAOA), an important enzyme in human physiology and behavior. The polymorphism, which is located 1.2 kb upstream of the MAOA coding sequences, consists of a 30-bp repeated sequence present in 3, 3.5, 4, or 5 copies. The polymorphism is in linkage disequilibrium with other MAOA and MAOB gene markers and displays significant variations in allele frequencies across ethnic groups. The polymorphism has been shown to affect the transcriptional activity of the MAOA gene promoter by gene fusion and transfection experiments involving three different cell types. Alleles with 3.5 or 4 copies of the repeat sequence are transcribed 2-10 times more efficiently than those with 3 or 5 copies of the repeat, suggesting an optimal length for the regulatory region. This promoter region polymorphism may be useful as both a functional and an anonymous genetic marker for MAOA.