Modifier genes and variation in cystic fibrosis

The availability of molecular tools to carry out genotyping has led to a flurry of association studies between specific genes and clinical indices of disease or disease susceptibility. Human studies, for the most part, have a limited number of subjects available, precluding whole genome types of approaches. 'Candidate gene' strategies have consequently become widespread, probably in part due to the inherent similarity to clinical association studies. Such studies in cystic fibrosis have found tantalizing results in genes involved in infection and inflammation, but many other relevant pathways remain untapped. Genome scanning approaches may eventually uncover genes not currently recognized as important to cystic fibrosis. In the meantime, while thousands of polymorphisms are cataloged and other genomic resources become more available, the number of association studies with candidate genes will no doubt increase. To make sense of these studies, the choice of gene and phenotype must be carefully considered.     

Correlates of substitution rate variation in mammalian protein-coding sequences

Background  

Rates of molecular evolution in different lineages can vary widely, and some of this variation might be predictable from aspects of species' biology. Investigating such predictable rate variation can help us to understand the causes of molecular evolution, and could also help to improve molecular dating methods. Here we present a comprehensive study of the life history correlates of substitution rate variation across the mammals, comparing results for mitochondrial and nuclear loci, and for synonymous and non-synonymous sites. We use phylogenetic comparative methods, refined to take into account the special nature of substitution rate data. Particular attention is paid to the widespread correlations between the components of mammalian life history, which can complicate the interpretation of results.     

Tandem repeat copy-number variation in protein-coding regions of human genes

Background

Tandem repeat variation in protein-coding regions will alter protein length and may introduce frameshifts. Tandem repeat variants are associated with variation in pathogenicity in bacteria and with human disease. We characterized tandem repeat polymorphism in human proteins, using the UniGene database, and tested whether these were associated with host defense roles.

Results

Protein-coding tandem repeat copy-number polymorphisms were detected in 249 tandem repeats found in 218 UniGene clusters; observed length differences ranged from 2 to 144 nucleotides, with unit copy lengths ranging from 2 to 57. This corresponded to 1.59% (218/13,749) of proteins investigated carrying detectable polymorphisms in the copy-number of protein-coding tandem repeats. We found no evidence that tandem repeat copy-number polymorphism was significantly elevated in defense-response proteins (p = 0.882). An association with the Gene Ontology term 'protein-binding' remained significant after covariate adjustment and correction for multiple testing. Combining this analysis with previous experimental evaluations of tandem repeat polymorphism, we estimate the approximate mean frequency of tandem repeat polymorphisms in human proteins to be 6%. Because 13.9% of the polymorphisms were not a multiple of three nucleotides, up to 1% of proteins may contain frameshifting tandem repeat polymorphisms.

Conclusion

Around 1 in 20 human proteins are likely to contain tandem repeat copy-number polymorphisms within coding regions. Such polymorphisms are not more frequent among defense-response proteins; their prevalence among protein-binding proteins may reflect lower selective constraints on their structural modification. The impact of frameshifting and longer copy-number variants on protein function and disease merits further investigation.     

Tandem repeat copy-number variation in protein-coding regions of human genes

Background  

Tandem repeat variation in protein-coding regions will alter protein length and may introduce frameshifts. Tandem repeat variants are associated with variation in pathogenicity in bacteria and with human disease. We characterized tandem repeat polymorphism in human proteins, using the UniGene database, and tested whether these were associated with host defense roles.     

Epistatic selection between coding and regulatory variation in human evolution and disease

Interaction (nonadditive effects) between genetic variants has been highlighted as an important mechanism underlying phenotypic variation, but the discovery of genetic interactions in humans has proved difficult. In this study, we show that the spectrum of variation in the human genome has been shaped by modifier effects of cis-regulatory variation on the functional impact of putatively deleterious protein-coding variants. We analyzed 1000 Genomes population-scale resequencing data from Europe (CEU [Utah residents with Northern and Western European ancestry from the CEPH collection]) and Africa (YRI [Yoruba in Ibadan, Nigeria]) together with gene expression data from arrays and RNA sequencing for the same samples. We observed an underrepresentation of derived putatively functional coding variation on the more highly expressed regulatory haplotype, which suggests stronger purifying selection against deleterious coding variants that have increased penetrance because of their regulatory background. Furthermore, the frequency spectrum and impact size distribution of common regulatory polymorphisms (eQTLs) appear to be shaped in order to minimize the selective disadvantage of having deleterious coding mutations on the more highly expressed haplotype. Interestingly, eQTLs explaining common disease GWAS signals showed an enrichment of putative epistatic effects, suggesting that some disease associations might arise from interactions increasing the penetrance of rare coding variants. In conclusion, our results indicate that regulatory and coding variants often modify the functional impact of each other. This specific type of genetic interaction is detectable from sequencing data in a genome-wide manner, and characterizing these joint effects might help us understand functional mechanisms behind genetic associations to human phenotypes-including both Mendelian and common disease.     

Interaction between area effects and variation with habitat in Cepaea

Populations of Cepaea nemoralis in Warwickshire occupy habitats of considerable temporal stability, most being at least 250 years old, and some much older. As expected from earlier work, shell pattern polymorphisms in these populations show variation with habitat of a kind suggesting the operation of visual selection for crypsis. They also show patterns of microgeographical variation unrelated to habitat. Although of a less extreme character, this variation resembles the 'area effects' seen in downland populations of Cepaea , in the lack of coincidence of variation at different loci, and in the existence of stronger and larger scale patterns in banding than in colour morphs. A similar explanation is advanced for their occurrence: previous bottlenecks and colonization from small relicts with founder effects. The less marked character of the variation is expected from the greater habitat stability and continuity in Warwickshire compared with downland. A clear colonization effect is seen in the inverse relationship between age of habitat and frequency of yellow in woodland populations. The results suggest that population histories affect variation in Cepaea even in areas of relative habitat stability.     

Partitioning of genetic variation between regulatory and coding gene segments: the predominance of software variation in genes encoding introvert proteins

Summary  In considering genetic variation in eukaryotes, a fundamental distinction can be made between variation in regulatory (software) and coding (hardware) gene segments. For quantitative traits the bulk of variation, particularly that near the population mean, appears to reside in regulatory segments. The main exceptions to this rule concern proteins which handle extrinsic substances, here termed extrovert proteins. The immune system includes an unsually large proportion of this exceptional category, but even so its chief source of variation may well be polymorphism in regulatory gene segments. The main evidence for this view emerges from genome scanning for quantitative trait loci (QTL), which in the case of the immune system points to a major contribution of pro-inflammatory cytokine genes. Further support comes from sequencing of major histocompatibility complex (Mhc) class II promoters, where a high level of polymorphism has been detected. These Mhc promoters appear to act, in part at least, by gating the back-signal from T cells into antigen-presenting cells. Both these forms of polymorphism are likely to be sustained by the need for flexibility in the immune response. Future work on promoter polymorphism is likely to benefit from the input from genome informatics.     

Nucleotide correspondence between protein-coding sequences of Helicobacter pylori 26695 and J99 strains

Comparison of open-reading frames (ORFs) H. pylori 26695 and J99 strains has been revealed prevalence of nucleotide replacements as transitions (more than 3%) above transversions (less than 1%). Prevalence of nucleotide transitions is caused by high speed of C : G to T : A transitions in a coding strand of DNA (3.5-5.3%) and not coding strand (2.9-3.9%). The correspondence rate of transversion (A --> C, A --> T, C --> A, C --> G, G --> C, G --> T, T --> A and T --> G) did not exceed 0.84%. The highest correspondence frequency between C and T was detected in ACGT-ATGT (28.3%) - the site of methylation by active methyltransferase M.Hpy99XI in H. pylori 26695 and J99. Thus one can speculate that predominant transition taking place in H. pylori is mutation of C into T, which is realized through cytosine methylation-deamination mechanism.     

Evolutionary differences between alternative and constitutive protein-coding regions of alternatively spliced genes of Drosophila

A total of 790 Drosophila melanogaster genes that are alternatively spliced in a coding region and have orthologs in Drosophila pseudoobscura were studied. It proved that nucleotide substitutions are accumulated in alternative coding regions more rapidly than in constitutive coding regions. Moreover, the evolutionary patterns of alternative regions differing in insertion-deletion mechanisms (use of alternative promoters, splicing sites, or polyadenylation sites) differ significantly. The synonymous substitution rate in coding regions of genes varies more strongly than the nonsynonymous substitution rate. The patterns of substitutions in different classes of alternative regions of Drosophila melanogaster and mammals differ considerably.     

Rare and common regulatory variation in population-scale sequenced human genomes

Population-scale genome sequencing allows the characterization of functional effects of a broad spectrum of genetic variants underlying human phenotypic variation. Here, we investigate the influence of rare and common genetic variants on gene expression patterns, using variants identified from sequencing data from the 1000 genomes project in an African and European population sample and gene expression data from lymphoblastoid cell lines. We detect comparable numbers of expression quantitative trait loci (eQTLs) when compared to genotypes obtained from HapMap 3, but as many as 80% of the top expression quantitative trait variants (eQTVs) discovered from 1000 genomes data are novel. The properties of the newly discovered variants suggest that mapping common causal regulatory variants is challenging even with full resequencing data; however, we observe significant enrichment of regulatory effects in splice-site and nonsense variants. Using RNA sequencing data, we show that 46.2% of nonsynonymous variants are differentially expressed in at least one individual in our sample, creating widespread potential for interactions between functional protein-coding and regulatory variants. We also use allele-specific expression to identify putative rare causal regulatory variants. Furthermore, we demonstrate that outlier expression values can be due to rare variant effects, and we approximate the number of such effects harboured in an individual by effect size. Our results demonstrate that integration of genomic and RNA sequencing analyses allows for the joint assessment of genome sequence and genome function.     

Accounting for evolutionary rate variation among sequence sites consistently changes universal phylogenies deduced from rRNA and protein-coding genes.

Phylogenetic analyses of gene and protein sequences have led to two major competing views of the universal phylogeny, the evolutionary tree relating the three kinds of living organisms, Bacteria, Archaea, and Eukarya. In the first scheme, called "the archaebacterial tree, " organisms of the same type are clustered together. In the second scenario, called "the eocyte tree," the archaeal phylum of Crenarchaeota is more closely related to eukaryotes than are other Archaea. A major property of the evolution of functional ribosomal and protein-encoding genes is that the rate of nucleotide and amino acid substitution varies across sequence sites. Here, using distance-based and maximum-likelihood methods, we show that universal phylogenies of ribosomal RNAs and RNA polymerases built by ignoring this variation are biased toward the archaebacterial tree because of attraction between long branches. In contrast, taking among-site rate variability into account gives support for the eocyte tree.     

Genetic diversity of Asian water buffalo (Bubalus bubalis): microsatellite variation and a comparison with protein-coding loci

Twenty-one microsatellite loci in 11 populations of Asian water buffalo (eight swamp, three river type) were analysed and, within and among populations, genetic variability was compared with results from 25 polymorphic protein-coding loci. Within-population mean heterozygosity ranged from 0·380–0·615, approximately twice that estimated from the protein-coding loci (0·184– 0·346). Only eight significant departures from Hardy–Weinberg equilibrium (involving four loci) were detected; global tests showed significant heterozygote deficiencies for these four loci. Non-amplifying alleles are likely to be segregating in some or all populations for one of these loci, and probably for the other three. There was significant differentiation between the swamp and river types of water buffalo, and among populations within each buffalo type. Estimates of θ (measure of population differentiation) for each locus for the eight swamp populations were all highly significant (mean θ = 0·168 ± 0·018). Mean θ for protein-coding loci was not significantly different (0·182 ± 0·041). The variance among protein-coding loci was significantly higher than among microsatellite loci, suggesting balancing selection affecting allele frequencies at some protein-coding loci. Genetic distances show clear separation of the swamp and river types, which were estimated to have diverged at least 10 000–15 000 years ago. The topology of the swamp populations’ microsatellite tree is consistent with their geographical distribution and their presumed spread through south-east Asia. By contrast, the tree based on the protein-coding loci distances is quite different, being clearly distorted by a bottleneck effect in one population, and possibly in at least two others. As many domestic livestock breeds are possibly descended from small numbers of founders, microsatellite-based trees are to be preferred in assessing breed genetic relationships.     

A variation in the structure of the protein-coding region of the human p53 gene

An extensive analysis of genomic DNA preparations from a number of normal and malignant tissues revealed BglII site polymorphism of the human p53 gene. Approximately 10% of p53 gene alleles were found to contain an additional BglII site localized in a region of intron I. This allelic form of p53 gene was also responsible for p53 protein having altered electrophoretic mobility. Molecular cloning and sequencing of both the alleles of p53 gene revealed a base-pair change in codon 72 causing arginine → proline substitution in the allele with the additional BglII site. Both variants of the p53 gene may occur in homozygous state and are therefore functional.