Estimating a nucleotide substitution rate for maize from polymorphism at a major domestication locus

To estimate a rate for single nucleotide substitutions for maize (Zea mays ssp. mays), we have taken advantage of data from genetic and archaeological studies of the domestication of maize from its wild ancestor, teosinte (Z. mays ssp. parviglumis). Genetic studies have shown that the teosinte branched1 (tb1) gene was a major target of human selection during maize domestication, and sequence diversity in the intergenic region 5' to the tb1-coding sequence is extraordinarily low. We show that polymorphism in this region is consistent with new mutation following fixation for a small number of tb1 haplotypes during domestication. Archeological studies suggest that maize was domesticated approximately 6,250-10,000 years ago and subsequently the size of the maize population is thought to have expanded rapidly. Using the observed number of mutations within the region of selection at tb1, the approximate age of maize domestication, and approximations for the maize genealogy, we have derived estimates for the nucleotide substitution rate for the tb1 intergenic region. Using two approaches, one of which is a coalescent approach, we obtain rate estimates of approximately 2.9 x 10(-8) and 3.3 x 10(-8) substitutions per site per year. We also show that the pattern of polymorphism in the tb1 intergenic region appears to have been strongly affected by the mutagenic effect of DNA methylation. Excluding target sites of symmetric DNA methylation (CG and CNG sites) from analysis, the mutation rate estimates are reduced by approximately 50%-60%, while the rates for CG and CNG sites are nearly an order of magnitude higher. We use rate estimates from the tb1 region to estimate the timing of expansion of transposable elements in the maize genome and suggest that this expansion occurred primarily within the last million years.     

Population variation in linkage disequilibrium across the COMT gene considering promoter region and coding region variation

Catechol-O-methyl transferase (COMT) catalyzes the first step in one of the major pathways in the degradation of catecholamines. The COMT gene on chromosome 22 has been considered a candidate gene for many neuropsychiatric disorders, in part because an exon 4 single nucleotide polymorphism (SNP) in COMT causes an amino acid substitution associated with significantly altered enzyme activity. This functional variant, detected as an NlaIII restriction site polymorphism (RSP), is polymorphic in populations from around the world. A four-site haplotype spanning 28 kb effectively encompasses COMT. This haplotype is comprised of two novel polymorphisms [a tetranucleotide short tandem repeat polymorphism (STRP) in intron 1 and a HindIII RSP at the 5' end of COMT], the NlaIII site, and another previously published site - a BglI RSP at the 3' end of the gene. Overall linkage disequilibrium (LD) for this haplotype is strong and significant in 32 population samples from around the world. Conditional probabilities indicate that, in spite of moderate to strong disequilibrium in most non-African populations, the NlaIII site, although often used for prediction, would not always be a reliable predictor of allelic variation at the other sites. Because other functional variation might exist, especially regulatory variation, these findings indicate that haplotypes would be more effective indicators of possible involvement of COMT in disease etiology.     

The brown midrib3 (bm3) mutation in maize occurs in the gene encoding caffeic acid O-methyltransferase.

The brown midrib mutations are among the earliest described in maize. Plants containing a brown midrib mutation exhibit a reddish brown pigmentation of the leaf midrib starting when there are four to six leaves. These mutations are known to alter lignin composition and digestibility of plants and therefore constitute prime candidates in the breeding of silage maize. Here, we show that two independent brown midrib3 (bm3) mutations have resulted from structural changes in the COMT gene, which encodes the enzyme O-methyltransferase (COMT; EC 2.1.1.6), involved in lignin biosynthesis. Our results indicate that the bm3-1 allele (the reference mutant allele) has arisen from an insertional event producing a COMT mRNA altered in both size and amount. By sequencing a COMT cDNA clone obtained from bm3-1 maize, a retrotransposon with homology to the B5 element has been found to be inserted near the junction of the 3' coding region of the COMT gene intron. The second bm3 allele, bm3-2, has resulted from a deletion of part of the COMT gene. These alterations of the COMT gene were confirmed by DNA gel blot and polymerase chain reaction amplification analyses. These results clearly demonstrate that mutations at the COMT gene give a brown midrib3 phenotype. Thus, the gene genetically recognized as bm3 is the same as the one coding for COMT.     

Nucleotide Polymorphism in the Adh1 Locus of Pearl Millet (Pennisetum Glaucum) (Poaceae)

We investigated nucleotide polymorphism in the Adh1 locus of pearl millet (Pennisetum glaucum) (Poaceae) by determining the DNA sequence of 20 alleles from 10 individuals. The individuals were sampled from throughout pearl millet's indigenous range and represent both wild and cultivated accessions. Our results indicated that there is little nucleotide polymorphism in the Adh1 locus. Estimates of per site nucleotide polymorphism did not differ significantly between cultivated and wild millet accessions. We compared nucleotide polymorphism in pearl millet Adh1 with nucleotide polymorphism in maize (Zea mays) Adh1 and conclude that the maize Adh1 sample is more polymorphic. Increased polymorphism in maize Adh1 may be attributable, in part, to faster substitution rates in the maize lineage. Analysis suggests that substitution rates in the maize Adh1 lineage are ~1.7 times faster than substitution rates in the millet Adh1 lineage.     

Catechol-O-methyltransferase: effects of the val108met polymorphism on protein turnover in human cells

A single nucleotide polymorphism in the human COMT (catechol-O-methyltransferase) gene has been associated with increased risk for breast cancer and several CNS diseases and disorders. The G to A polymorphism causes a valine (val) to methionine (met) substitution at codon 108 soluble - (S)/158 membrane - (MB)-COMT, generating alleles encoding high and low-activity forms of the enzyme, COMT H and COMT L, respectively. Tissues and cells with a COMT LL genotype have decreased COMT activity compared to COMT HH cells. Previously, we reported that the decreased activity was due to decreased amounts of S-COMT L protein in human hepatocytes. In this study, we investigated the role of S-COMT protein synthesis and turnover as determinates of reduced COMT protein in COMT LL compared to COMT HH cells. No association between S-COMT protein synthesis and COMT genotype was detected. Using a pulse-chase protocol, the half-life of S-COMT H was determined to be 4.7 days, which was considerably longer than expected from the half-lives of other phase 2 enzyme proteins. The half-life of S-COMT L compared to S-COMT H protein was significantly shorter at 3.0 days, but the difference was affected by the medium used during the chase period. These results suggest that increased turnover may contribute to reduced COMT activity in cells and tissues from COMT LL individuals. Subtle differences appear to be able to affect the stability of the S-COMT L protein, and this may contribute to the differences observed in epidemiological studies on the association of this polymorphism with breast cancer risk.     

Distribution of the Val108/158Met polymorphism of the COMT gene in healthy Mexican population

Catechol-O-methyltransferase (COMT) inactivates the catecholamines adrenaline, noradrenaline and dopamine. On the other hand, some studies have reported that the enzymatic activity of COMT is partly genetically determined. With regard to the COMT gene, the most studied polymorphism is the functional variant Val108/158Met (rs4680), which results in substantial three- to four-fold variations in enzyme activity. To date, the rs4680 polymorphism of COMT has been associated with a number of disorders. In addition, this polymorphism has been found to have important differences in frequency according to the studied population. Therefore, the aim of the present study was to evaluate the frequency of a common single nucleotide polymorphism (SNP) Val108/158Met of the COMT gene in the Mexican population. Accordingly, we recruited 431 healthy volunteers. Our sample consisted of 111 healthy individuals from Mexico City and 320 individuals from the state of Tabasco, Mexico. We observed that Met was the most common allele, ranging from 57% (Tabasco) to 85% (Mexico City). In addition, we analyzed the frequency of Val108/158Met polymorphism of Caucasian (54% Met allele), Asian (29% Met allele) and African (34% Met allele) populations separately and also in comparison with Mexican (63% Met allele) population. In conclusion, the distribution of the Val108/158Met polymorphism distinguishes the Mexican population studied from other populations, but it is necessary to increase the size of the sample to get more conclusive results.     

The molecular evolution of terminal ear1, a regulatory gene in the genus Zea.

Nucleotide diversity in the terminal ear1 (te1) gene, a regulatory locus hypothesized to be involved in the morphological evolution of maize (Zea mays ssp. mays), was investigated for evidence of past selection. Nucleotide polymorphism in a 1.4-kb region of te1 was analyzed for a sample of 26 sequences isolated from 12 maize lines, five populations of the maize progenitor, Z. mays ssp. parviglumis, six other Zea populations, and two Tripsacum species. Although nucleotide diversity in te1 in maize is reduced relative to ssp. parviglumis, phylogenetic and statistical analyses of the pattern of polymorphism among these sequences provided no evidence of past selection, indicating that the region of the gene studied was probably not involved in maize evolution. The level of reduction in genetic diversity in te1 in maize relative to its progenitor is comparable to that found in previous reports for isozymes and other neutrally evolving maize genes and is consistent with a genome-wide reduction of genetic diversity resulting from a domestication bottleneck. An estimate of the age (1.2-1.4 million yr) of the maize gene pool based on te1 is roughly consistent with previous estimates based on other neutral genes, but may be biased by the apparently slow synonymous substitution rate at te1.     

Cloning and characterization of a putative GS3 ortholog involved in maize kernel development

The GS3 gene was the first identified gene controlling the grain size in rice. It has been proven to be involved in the evolution of grain size during domestication. We isolated the maize ortholog, ZmGS3 and investigated its role in the evolution of maize grain size. ZmGS3 has five exons encoding a protein with 198 amino acids, and has domains in common with the rice GS3 protein. Compared with teosinte, maize has reduced nucleotide diversity at ZmGS3, and the reduction is comparable to that found in neutrally evolving maize genes. No positive selection was detected along the length of the gene using either the Hudson–Kreitman–Aguadé or Tajima’s D tests. Phylogenetic analysis reveals a distribution of maize sequences among two different clades, with one clade including related teosinte sequences. The nucleotide polymorphism analysis, selection test and phylogenetic analysis reveal that ZmGS3 has not been subjected to selection, and appears to be a neutrally evolving gene. In maize, ZmGS3 is primarily expressed in immature ears and kernels, implying a role in maize kernel development. Association mapping analysis revealed one polymorphism in the fifth exon that is significantly associated with kernel length in two environments. Also one polymorphism in the promoter region was found to affect hundred kernel weight in both environments. Collectively, these results imply that ZmGS3 is involved in maize kernel development but with different functional polymorphisms and thus, possibly different mechanisms from that of the rice GS3 gene.     

Nucleotide diversity and linkage disequilibrium in wild avocado (Persea americana Mill.)

Resequencing studies provide the ultimate resolution of genetic diversity because they identify all mutations in a gene that are present within the sampled individuals. We report a resequencing study of Persea americana, a subtropical tree species native to Meso- and Central America and the progenitor of cultivated avocado. The sample includes 21 wild accessions from Mexico, Costa Rica, Ecuador, and the Dominican Republic. Estimated levels of nucleotide polymorphism and linkage disequilibrium (LD) are obtained from fully resolved haplotype data from 4 nuclear loci that span 5960 nucleotide sites. Results show that, although avocado is a subtropical tree crop and a predominantly outcrossing plant, the overall level of genetic variation is not exceptionally high (nucleotide diversity at silent sites, pi(sil) = 0.0102) compared with available estimates from temperate plant species. Intralocus LD decays rapidly to half the initial value within about 1 kb. Estimates of recombination rate (based on the sequence data) show that the rate is not exceptionally high when compared with annual plants such as wild barley or maize. Interlocus LD is significant owing to substantial population structure induced by mixing of the 3 botanical races of avocado.     

Pathways to Age of Onset of Heroin Use: A Structural Model Approach Exploring the Relationship of the COMT Gene,Impulsivity and Childhood Trauma

Background

The interaction of the association of dopamine genes, impulsivity and childhood trauma with substance abuse remains unclear.

Objectives

To clarify the impacts and the interactions of the Catechol -O-methyltransferase (COMT) gene, impulsivity and childhood trauma on the age of onset of heroin use among heroin dependent patients in China.

Methods

202 male and 248 female inpatients who meet DSM-IV criteria of heroin dependence were enrolled. Impulsivity and childhood trauma were measured using BIS-11 (Barratt Impulsiveness Scale-11) and ETISR-SF (Early Trauma Inventory Self Report-Short Form). The single nucleotide polymorphism (SNP) rs737866 on the COMT gene-which has previously been associated with heroin abuse, was genotyped using a DNA sequence detection system. Structural equations model was used to assess the interaction paths between these factors and the age of onset of heroin use.

Principal Findings

Chi-square test indicated the individuals with TT allele have earlier age of onset of heroin use than those with CT or CC allele. In the correlation analysis, the severity of childhood trauma was positively correlated to impulsive score, but both of them were negatively related to the age of onset of heroin use. In structure equation model, both the COMT gene and childhood trauma had impacts on the age of onset of heroin use directly or via impulsive personality.

Conclusions

Our findings indicated that the COMT gene, impulsive personality traits and childhood trauma experience were interacted to impact the age of onset of heroin use, which play a critical role in the development of heroin dependence. The impact of environmental factor was greater than the COMT gene in the development of heroin dependence.     

Cloning of cDNA encoding COMT from wheat which is differentially expressed in lodging-sensitive and -resistant cultivars

In the present study two cDNA fragments were cloned by nested-PCR using degenerate primers for COMT and found to be 93% identical at the nucleotide level. The deduced amino acid sequence of the two cDNAs showed a high degree of identity with COMT from other plants and were most similar to COMTs from monocots. RNA gel blot hybridization demonstrated that the wheat COMT gene W-cm5-1 was expressed in stem, root and leaf tissues. W-cm5-1 mRNA levels in the developing wheat stem were associated with the stem lodging character in two wheat cultivars. These results suggest that the action of the wheat COMT gene may be related to stem rigidity and lodging character in wheat.     

Isolation and sequence of a gene encoding a methionine-rich 10-kDa zein protein from maize

We have isolated the gene encoding a methionine-rich 10-kDa zein protein from a lambda EMBL3 maize genomic 'mini' library of the inbred line BSSS-53 and determined its nucleotide sequence. The sequence matches perfectly with a cDNA clone from the inbred line W22 (which has the same restriction fragment length polymorphism as many inbred lines tested) indicating that we have isolated a functional storage protein gene that is very conserved in maize. This comparison also excludes any splicing of any precursor mRNA and therefore any presence of introns. A number of potential regulatory sequences have been located in the flanking regions. The 10-kDa-zein gene represents the last size class in the zein multigene family to be characterized. Its structure allows us now to re-examine the relationship of all the zein proteins and also to compare the structure of a new class of storage proteins that are rich in methionine, an essential amino acid in livestock fodder.     

Structure and expression of the lignin O-methyltransferase gene from Zea mays L.

The isolation and characterization of cDNA and homologous genomic clones encoding the lignin O-methyltransferase (OMT) from maize is reported. The cDNA clone has been isolated by differential screening of maize root cDNA library. Southern analysis indicates that a single gene codes for this protein. The genomic sequence contains a single 916 bp intron. The deduced protein sequence from DNA shares significant homology with the recently reported lignin-bispecific caffeic acid/5-hydroxyferulic OMTs from alfalfa and aspen. It also shares homology with OMTs from bovine pineal glands and a purple non-sulfur photosynthetic bacterium. The mRNA of this gene is present at different levels in distinct organs of the plant with the highest accumulation detected in the elongation zone of roots. Bacterial extracts from clones containing the maize OMT cDNA show an activity in methylation of caffeic acid to ferulic acid comparable to that existing in the plant extracts. These results indicate that the described gene encodes the caffeic acid 3-O-methyltransferase (COMT) involved in the lignin biosynthesis of maize.     

Molecular characterization of a brown midrib3 deletion mutation in maize

The caffeic acid O-methyltransferase (COMT) gene plays an important role in the synthesis of lignin. We have used the polymerase chain reaction in conjuction with genomic analysis to characterize deletion mutations of this gene in maize. In addition, we have analyzed and compared regions of the COMT gene from three distinct heterotic groups. Both PCR and Southern analysis indicate that the active wild-type COMT gene can be polymorphic. We suggest that the intron domain of at least one heterotic inbred can contribute to the alteration of the wild-type gene. In addition, multiple deletion mutations have occurred at this locus. We have found a previously uncharacterized deletion mutation in which segments of both the intron and exon have been deleted and replaced by other sequences. Precise knowledge of its sequence has allowed us to develop an assay by which we can follow this mutation in a breeding program.     

DNA Sequence Variation within Maize and Melon: Observations from Polymerase Chain Reaction Amplification and Direct Sequencing

While compiling genetic linkage maps in several plant species based upon restriction fragment length polymorphisms (RFLPs), it was noted that the incidence of polymorphism differs among species. The basis of this disparity was investigated in this study by examining the nucleotide sequence at homologous loci among distinct cultivars within two species which exhibit considerably different levels of RFLPs. Using the polymerase chain reaction, homologous regions from different cultivars were first amplified and the nucleotide sequence of the products were determined. Four genomic regions of seven maize cultivars and three genomic regions of eight melon cultivars were examined to compare the respective levels of sequence variation between the two species. Levels of variation for both base substitutions and insertions/deletions varied widely among the maize sequences and between maize and melon for base substitutions. Estimates of theta (a measure of polymorphism) ranged from 0 to 0.002 in melon and from 0.006 to 0.040 for base substitutions and from 0.002 to 0.023 for insertions/deletions in maize. Critical value tests and chi-squared tests suggested that in maize the underlying processes generating and maintaining neutral mutations differ among the regions. The results not only suggest that several mechanisms are necessary to explain the variation seen in these two species, but also point to some basic dissimilarities in the organization and maintenance of the genomes of different plant species.     

Population structure and genetic diversity in a commercial maize breeding program assessed with SSR and SNP markers

Information about the genetic diversity and population structure in elite breeding material is of fundamental importance for the improvement of crops. The objectives of our study were to (a) examine the population structure and the genetic diversity in elite maize germplasm based on simple sequence repeat (SSR) markers, (b) compare these results with those obtained from single nucleotide polymorphism (SNP) markers, and (c) compare the coancestry coefficient calculated from pedigree records with genetic distance estimates calculated from SSR and SNP markers. Our study was based on 1,537 elite maize inbred lines genotyped with 359 SSR and 8,244 SNP markers. The average number of alleles per locus, of group specific alleles, and the gene diversity (D) were higher for SSRs than for SNPs. Modified Roger’s distance (MRD) estimates and membership probabilities of the STRUCTURE matrices were higher for SSR than for SNP markers but the germplasm organization in four heterotic pools was consistent with STRUCTURE results based on SSRs and SNPs. MRD estimates calculated for the two marker systems were highly correlated (0.87). Our results suggested that the same conclusions regarding the structure and the diversity of heterotic pools could be drawn from both markers types. Furthermore, although our results suggested that the ratio of the number of SSRs and SNPs required to obtain MRD or D estimates with similar precision is not constant across the various precision levels, we propose that between 7 and 11 times more SNPs than SSRs should be used for analyzing population structure and genetic diversity.     

Contrasting patterns of synonymous and nonsynonymous sequence evolution in asexual and sexual freshwater snail lineages

In asexual lineages, both synonymous and nonsynonymous sequence polymorphism may be reduced due to severe founder effects when asexual lineages originate. However, mildly deleterious (nonsynonymous) mutations may accumulate after asexual lineages are formed, because the efficiency of purifying selection is reduced even in the nonrecombining mitochondrial genome. Here we examine patterns of synonymous and nonsynonymous mitochondrial sequence polymorphism in asexual and sexual lineages of the freshwater snail Campeloma. Using clade-specific estimates, we found that synonymous sequence polymorphism was significantly reduced by 75% in asexuals relative to sexuals, whereas nonsynonymous sequence polymorphism did not differ significantly between sexuals and asexuals. Two asexual clades had high negative values for Tajima's D statistic. Coalescent simulations confirmed that various bottleneck scenarios can account for this result. We also used branch-specific estimates of the ratio of amino acid to silent substitutions, K(a)/K(s). Our study revealed that K(a)/K(s) ratios are six times higher in terminal branches of independent asexual lineages compared to sexuals. Coalescent-based reconstruction of gene networks for all sexual and asexual clades indicated that nonsynonymous mutations occurred at a higher frequency in recently derived asexual haplotypes. These findings suggest that patterns of synonymous and nonsynonymous nucleotide polymorphism in asexual snail lineages may be shaped by both severe founder effect and relaxed purifying selection.