The PHF21B gene is associated with major depression and modulates the stress response

Major depressive disorder (MDD) affects around 350 million people worldwide; however, the underlying genetic basis remains largely unknown. In this study, we took into account that MDD is a gene-environment disorder, in which stress is a critical component, and used whole-genome screening of functional variants to investigate the 'missing heritability' in MDD. Genome-wide association studies (GWAS) using single- and multi-locus linear mixed-effect models were performed in a Los Angeles Mexican-American cohort (196 controls, 203 MDD) and in a replication European-ancestry cohort (499 controls, 473 MDD). Our analyses took into consideration the stress levels in the control populations. The Mexican-American controls, comprised primarily of recent immigrants, had high levels of stress due to acculturation issues and the European-ancestry controls with high stress levels were given higher weights in our analysis. We identified 44 common and rare functional variants associated with mild to moderate MDD in the Mexican-American cohort (genome-wide false discovery rate, FDR, <0.05), and their pathway analysis revealed that the three top overrepresented Gene Ontology (GO) processes were innate immune response, glutamate receptor signaling and detection of chemical stimulus in smell sensory perception. Rare variant analysis replicated the association of the PHF21B gene in the ethnically unrelated European-ancestry cohort. The TRPM2 gene, previously implicated in mood disorders, may also be considered replicated by our analyses. Whole-genome sequencing analyses of a subset of the cohorts revealed that European-ancestry individuals have a significantly reduced (50%) number of single nucleotide variants compared with Mexican-American individuals, and for this reason the role of rare variants may vary across populations. PHF21b variants contribute significantly to differences in the levels of expression of this gene in several brain areas, including the hippocampus. Furthermore, using an animal model of stress, we found that Phf21b hippocampal gene expression is significantly decreased in animals resilient to chronic restraint stress when compared with non-chronically stressed animals. Together, our results reveal that including stress level data enables the identification of novel rare functional variants associated with MDD.

The human melanocortin-1 receptor locus: analysis of transcription unit, locus polymorphism and haplotype evolution

The complete sequence of the MC1R locus has been assembled, the coding region of the gene is intronless and placed within a 12 kb region flanked by the NULP1 and TUBB4 genes. The immediate promoter region has an E-box site with homology to the M-box consensus known to bind the microphthalmia transcription factor (MITF); however, promoter deletion analysis and transactivation studies have failed to show activation through this element by MITF. Polymorphism within the coding region, immediate 5' promoter region and a variable number tandem repeat (VNTR) minisatellite within the locus have been examined in a collection of Caucasian families and African individuals. Haplotype analysis shows linkage disequilibrium between the VNTR and MC1R coding region red hair variant alleles which can be used to estimate the age of these missense changes. Assuming a mean VNTR mutation rate of 1% and a star phylogeny, we estimate the Arg151Cys variant arose 7500 years before the present day, suggesting these variants may have arisen in the Caucasian population more recently than previously thought.

Mitochondrial DNA sequences in ancient Australians: Implications for modern human origins

DNA from ancient human remains provides perspectives on the origin of our species and the relationship between molecular and morphological variation. We report analysis of mtDNA from the remains of 10 ancient Australians. These include the morphologically gracile Lake Mungo 3 [ approximately 60 thousand years (ka) before present] and three other gracile individuals from Holocene deposits at Willandra Lakes (<10 ka), all within the skeletal range of living Australians, and six Pleistocene/early Holocene individuals (15 to <8 ka) from Kow Swamp with robust morphologies outside the skeletal range of contemporary indigenous Australians. Lake Mungo 3 is the oldest (Pleistocene) "anatomically modern" human from whom DNA has been recovered. His mtDNA belonged to a lineage that only survives as a segment inserted into chromosome 11 of the nuclear genome, which is now widespread among human populations. This lineage probably diverged before the most recent common ancestor of contemporary human mitochondrial genomes. This timing of divergence implies that the deepest known mtDNA lineage from an anatomically modern human occurred in Australia; analysis restricted to living humans places the deepest branches in East Africa. The other ancient Australian individuals we examined have mtDNA sequences descended from the most recent common ancestor of living humans. Our results indicate that anatomically modern humans were present in Australia before the complete fixation of the mtDNA lineage now found in all living people. Sequences from additional ancient humans may further challenge current concepts of modern human origins.

Testing for concordant equilibria between population samples

A substantial body of theory has been developed to assess the effect of evolutionary forces on the distribution of genotypes, both single and multilocus, within populations. One area where the potential for application of this theory has not been fully appreciated concerns the extent to which population samples differ. Within populations, the divergence of genotype or haplotype frequencies from that expected under Hardy-Weinberg (HW) or linkage equilibrium can be measured as disequilibria coefficients. To assess population samples for concordant equilibria, an analytical framework for comparing disequilibria coefficients between populations is necessary. Here we present log-linear models to evaluate such hypotheses. These models have broad utility ranging from conventional population genetics to genetic epidemiology. We demonstrate the use of these log-linear models (1) as a test for genetic association with disease and (2) as a test for different levels of linkage disequilibria between human populations.

Adaptive evolution of the tumour suppressor BRCA1 in humans and chimpanzees. Australian Breast Cancer Family Study

Mutations in BRCA1 (ref. 1) confer an increased risk of female breast cancer. In a genome-wide scan of linkage disequilibrium (LD), a high level of LD was detected among microsatellite markers flanking BRCA1 (ref. 3), raising the prospect that positive natural selection may have acted on this gene. We have used the predictions of evolutionary genetic theory to investigate this further. Using phylogeny-based maximum likelihood analysis of the BRCA1 sequences from primates and other mammals, we found that the ratios of replacement to silent nucleotide substitutions on the human and chimpanzee lineages were not different from one another (P=0.8), were different from those of other primate lineages (P=0.004) and were greater than 1 (P=0.04). This is consistent with the historic occurrence of positive darwinian selection pressure on the BRCA1 protein in the human and chimpanzee lineages. Analysis of genetic variation in a sample of female Australians of Northern European origin showed evidence for Hardy-Weinberg (HW) disequilibrium at polymorphic sites in BRCA1, consistent with the possibility that natural selection is affecting genotype frequencies in modern Europeans. The clustering of between-species variation in the region of the gene encoding the RAD51-interaction domain of BRCA1 suggests the maintenance of genomic integrity as a possible target of selection.

How important is DNA replication for mutagenesis?

Rates of mutation and substitution in mammals are generally greater in the germ lines of males. This is usually explained as resulting from the larger number of germ cell divisions during spermatogenesis compared with oogenesis, with the assumption made that mutations occur primarily during DNA replication. However, the rate of cell division is not the only difference between male and female germ lines, and mechanisms are known that can give rise to mutations independently of DNA replication. We investigate the possibility that there are other causes of male-biased mutation. First, we show that patterns of variation at approximately 5,200 short tandem repeat (STR) loci indicate a higher mutation rate in males. We estimate a ratio of male-to-female mutation rates of approximately 1.9. This is significantly greater than 1 and supports a greater rate of mutation in males, affecting the evolution of these loci. Second, we show that there are chromosome-specific patterns of nucleotide and dinucleotide composition in mammals that have been shaped by mutation at CpG dinucleotides. Comparable patterns occur in birds. In mammals, male germ lines are more methylated than female germ lines, and these patterns indicate that differential methylation has played a role in male-biased vertebrate evolution. However, estimates of male mutation bias obtained from both classes of mutation are substantially lower than estimates of cell division bias from anatomical data. This discrepancy, along with published data indicating slipped-strand mispairing arising at STR loci in nonreplicating DNA, suggests that a substantial percentage of mutation may occur in nonreplicating DNA.

A scan for linkage disequilibrium across the human genome

Linkage disequilibrium (LD), the tendency for alleles of linked loci to co-occur nonrandomly on chromosomal haplotypes, is an increasingly useful phenomenon for (1) revealing historic perturbation of populations including founder effects, admixture, or incomplete selective sweeps; (2) estimating elapsed time since such events based on time-dependent decay of LD; and (3) disease and phenotype mapping, particularly for traits not amenable to traditional pedigree analysis. Because few descriptions of LD for most regions of the human genome exist, we searched the human genome for the amount and extent of LD among 5048 autosomal short tandem repeat polymorphism (STRP) loci ascertained as specific haplotypes in the European CEPH mapping families. Evidence is presented indicating that approximately 4% of STRP loci separated by <4.0 cM are in LD. The fraction of locus pairs within these intervals that display small Fisher's exact test (FET) probabilities is directly proportional to the inverse of recombination distance between them (1/cM). The distribution of LD is nonuniform on a chromosomal scale and in a marker density-independent fashion, with chromosomes 2, 15, and 18 being significantly different from the genome average. Furthermore, a stepwise (locus-by-locus) 5-cM sliding-window analysis across 22 autosomes revealed nine genomic regions (2.2-6.4 cM), where the frequency of small FET probabilities among loci was greater than or equal to that presented by the HLA on chromosome 6, a region known to have extensive LD. Although the spatial heterogeneity of LD we detect in Europeans is consistent with the operation of natural selection, absence of a formal test for such genomic scale data prevents eliminating neutral processes as the evolutionary origin of the LD.

Dating the origin of the CCR5-Delta32 AIDS-resistance allele by the coalescence of haplotypes

The CCR5-Delta32 deletion obliterates the CCR5 chemokine and the human immunodeficiency virus (HIV)-1 coreceptor on lymphoid cells, leading to strong resistance against HIV-1 infection and AIDS. A genotype survey of 4,166 individuals revealed a cline of CCR5-Delta32 allele frequencies of 0%-14% across Eurasia, whereas the variant is absent among native African, American Indian, and East Asian ethnic groups. Haplotype analysis of 192 Caucasian chromosomes revealed strong linkage disequilibrium between CCR5 and two microsatellite loci. By use of coalescence theory to interpret modern haplotype genealogy, we estimate the origin of the CCR5-Delta32-containing ancestral haplotype to be approximately 700 years ago, with an estimated range of 275-1,875 years. The geographic cline of CCR5-Delta32 frequencies and its recent emergence are consistent with a historic strong selective event (e.g. , an epidemic of a pathogen that, like HIV-1, utilizes CCR5), driving its frequency upward in ancestral Caucasian populations.

Contrasting genetic influence of CCR2 and CCR5 variants on HIV-1 infection and disease progression. Hemophilia Growth and Development Study (HGDS), Multicenter AIDS Cohort Study (MACS), Multicenter Hemophilia Cohort Study (MHCS), San Francisco City Cohort (SFCC), ALIVE Study

The critical role of chemokine receptors (CCR5 and CXCR4) in human immunodeficiency virus-type 1 (HIV-1) infection and pathogenesis prompted a search for polymorphisms in other chemokine receptor genes that mediate HIV-1 disease progression. A mutation (CCR2-64I) within the first transmembrane region of the CCR2 chemokine and HIV-1 receptor gene is described that occurred at an allele frequency of 10 to 15 percent among Caucasians and African Americans. Genetic association analysis of five acquired immunodeficiency syndrome (AIDS) cohorts (3003 patients) revealed that although CCR2-64I exerts no influence on the incidence of HIV-1 infection, HIV-1-infected individuals carrying the CCR2-64I allele progressed to AIDS 2 to 4 years later than individuals homozygous for the common allele. Because CCR2-64I occurs invariably on a CCR5-+-bearing chromosomal haplotype, the independent effects of CCR5-Delta32 (which also delays AIDS onset) and CCR2-64I were determined. An estimated 38 to 45 percent of AIDS patients whose disease progresses rapidly (less than 3 years until onset of AIDS symptoms after HIV-1 exposure) can be attributed to their CCR2-+/+ or CCR5-+/+ genotype, whereas the survival of 28 to 29 percent of long-term survivors, who avoid AIDS for 16 years or more, can be explained by a mutant genotype for CCR2 or CCR5.

Genetic restriction of HIV-1 infection and progression to AIDS by a deletion allele of the CKR5 structural gene. Hemophilia Growth and Development Study, Multicenter AIDS Cohort Study, Multicenter Hemophilia Cohort Study, San Francisco City Cohort, ALIVE Study

The chemokine receptor 5 (CKR5) protein serves as a secondary receptor on CD4(+) T lymphocytes for certain strains of human immunodeficiency virus-type 1 (HIV-1). The CKR5 structural gene was mapped to human chromosome 3p21, and a 32-base pair deletion allele (CKR5Delta32) was identified that is present at a frequency of approximately0.10 in the Caucasian population of the United States. An examination of 1955 patients included among six well-characterized acquired immunodeficiency syndrome (AIDS) cohort studies revealed that 17 deletion homozygotes occurred exclusively among 612 exposed HIV-1 antibody-negative individuals (2.8 percent) and not at all in 1343 HIV-1-infected individuals. The frequency of CKR5 deletion heterozygotes was significantly elevated in groups of individuals that had survived HIV-1 infection for more than 10 years, and, in some risk groups, twice as frequent as their occurrence in rapid progressors to AIDS. Survival analysis clearly shows that disease progression is slower in CKR5 deletion heterozygotes than in individuals homozygous for the normal CKR5 gene. The CKR5Delta32 deletion may act as a recessive restriction gene against HIV-1 infection and may exert a dominant phenotype of delaying progression to AIDS among infected individuals.

Evolution of the chalcone synthase gene family in the genus Ipomoea

The evolution of the chalcone synthase [CHS; malonyl-CoA:4-coumaroyl-CoA malonyltransferase (cyclizing), EC 2.3.1.74] multigene family in the genus Ipomoea is explored. Thirteen CHS genes from seven Ipomoea species (family Convolvulaceae) were sequenced--three from genomic clones and the remainder from PCR amplification with primers designed from the 5' flanking region and the end of the 3' coding region of Ipomoea purpurea Roth. Analysis of the data indicates a duplication of CHS that predates the divergence of the Ipomoea species in this study. The Ipomoea CHS genes are among the most rapidly evolving of the CHS genes sequenced to date. The CHS genes in this study are most closely related to the Petunia CHS-B gene, which is also rapidly evolving and highly divergent from the rest of the Petunia CHS sequences.

Molecular evolution of the Ac/Ds transposable-element family in pearl millet and other grasses

We report an Ac-like sequence from pearl millet (Pennisetum glaucum) and deletion derivative Ac-like sequences from pearl millet and another grass species, Bambusa multiplex. Sequence relationships between the pearl millet and maize Ac elements suggest that Ac/Ds transposable-element family is ancient. Further, the sequence identity between the Bambusa Ac-like sequence and maize Ac implies that the Ac/Ds transposable-element family has been in the grass family since its inception. The Ac-like sequences reported from pearl millet and maize Ac are statistically heterogeneous in pair-wise distance comparisons to each other. Yet, we are unable to discriminate between differential selection or ectopic exchange (recombination and conversion) between nonidentical transposable element homologues, as the cause of the heterogeneity. However, the more extreme heterogeneity exhibited between the previously described pearl millet element and maize Ac seems likely to derive from ectopic exchange between elements with different levels of divergence.

Molecular evolutionary characterization of an Activator (Ac)-like transposable element sequence from pearl millet (Pennisetum glaucum) (Poaceae)

We present data on the evolution of the Ac/Ds family of transposable elements in select grasses (Poaceae). A defective Ac-like element was cloned from a DNA library of the grass Pennisetum glaucum (pearl millet) and its entire 4531 bp sequence has been determined. When the pearl millet Ac-like sequence is aligned with the maize Ac sequence, it is found that there is approximately 70% DNA similarity in the central region spanning most of maize Ac exon II and all of exon III. In addition, there are two smaller regions of similarity at the Ac terminii. Besides these three major structural similarities, Pennisetum Ac has two large regions, one 5' and one 3', that show little similarity to Zea Ac. Furthermore, most of the sequences corresponding to intron II in maize Ac are absent in pearl millet Ac. Kimura's evolutionary distance between the central region of maize and pearl millet Ac sequences is estimated to be 0.429 +/- 0.020 nucleotide substitutions per site. This value is not significantly different from the average number of synonymous substitutions for coding regions of the Adh1 gene between maize and pearl millet, which is 0.395 +/- 0.051 nucleotide substitutions per site. If we assume Ac and Adh1 divergence times are equivalent between maize and pearl millet, then the above calculations suggest Ac-like sequences have probably not been strongly constrained by natural selection. Conserved DNA and amino acid sequence motifs are also examined. The level of DNA sequence divergence between maize and pearl millet Ac sequences, the estimated date when maize and pearl millet diverged (25-40 million years ago), coupled with their reproductive isolation/lack of current genetic exchange, all support the theory that Ac-like sequences have not been recently introduced into pearl millet from maize. Instead, Ac-like sequences were probably present in the progenitor of maize and pearl millet and have thus existed in the grasses for at least 25 million years.