Candidate gene identification following linkage: Search for hypertension susceptibility genes

Following up genetic linkage studies to identify the underlying susceptibility gene(s) for complex disease traits is an arduous yet biologically and clinically important task. Complex traits, such as hypertension, are considered polygenic with many genes influencing risk, each with small effects. Chromosome 2 has been consistently identified as a genomic region with genetic linkage evidence suggesting that one or more loci contribute to blood pressure levels and hypertension status. Using combined positional candidate gene methods, the Family Blood Pressure Program has concentrated efforts in investigating this region of chromosome 2 in an effort to identify underlying candidate hypertension susceptibility gene(s). Initial informatics efforts identified the boundaries of the region and the known genes within it. A total of 82 polymorphic sites in eight positional candidate genes were genotyped in a large hypothesis-generating sample consisting of 1640 African Americans, 1339 whites, and 1616 Mexican Americans. To adjust for multiple comparisons, resampling-based false discovery adjustment was applied, extending traditional resampling methods to sibship samples. Following this adjustment for multiple comparisons, SLC4A5, a sodium bicarbonate transporter, was identified as a primary candidate gene for hypertension. Polymorphisms in SLC4A5 were subsequently genotyped and analyzed for validation in two populations of African Americans (N = 461; N = 778) and two of whites (N = 550; N = 967). Again, SNPs within SLC4A5 were significantly associated with blood pressure levels and hypertension status. While not identifying a single causal DNA sequence variation that is significantly associated with blood pressure levels and hypertension status across all samples, the results further implicate SLC4A5 as a candidate hypertension susceptibility gene, validating previous evidence for one or more genes on chromosome 2 that influence hypertension related phenotypes in the population-at-large. The methodology and results reported provide a case study of one approach for following up the results of genetic linkage analyses to identify genes influencing complex traits. ^

Genome-wide gene-gene interaction analysis for cardiovascular disease

Numerous studies have been carried out to try to better understand the genetic predisposition for cardiovascular disease. Although it is widely believed that multifactorial diseases such as cardiovascular disease is the result from effects of many genes which working alone or interact with other genes, most genetic studies have been focused on identifying of cardiovascular disease susceptibility genes and usually ignore the effects of gene-gene interactions in the analysis. The current study applies a novel linkage disequilibrium based statistic for testing interactions between two linked loci using data from a genome-wide study of cardiovascular disease. A total of 53,394 single nucleotide polymorphisms (SNPs) are tested for pair-wise interactions, and 8,644 interactions are found to be significant with p-values less than 3.5×10-11. Results indicate that known cardiovascular disease susceptibility genes tend not to have many significantly interactions. One SNP in the CACNG1 (calcium channel, voltage-dependent, gamma subunit 1) gene and one SNP in the IL3RA (interleukin 3 receptor, alpha) gene are found to have the most significant pair-wise interactions. Findings from the current study should be replicated in other independent cohort to eliminate potential false positive results.^

Single nucleotide polymorphisms (SNPs) associated with TGF-beta pathway and their significance in systemic sclerosis - A multilevel analysis

Systemic sclerosis (SSc) or Scleroderma is a complex disease and its etiopathogenesis remains unelucidated. Fibrosis in multiple organs is a key feature of SSc and studies have shown that transforming growth factor-β (TGF-β) pathway has a crucial role in fibrotic responses. For a complex disease such as SSc, expression quantitative trait loci (eQTL) analysis is a powerful tool for identifying genetic variations that affect expression of genes involved in this disease. In this study, a multilevel model is described to perform a multivariate eQTL for identifying genetic variation (SNPs) specifically associated with the expression of three members of TGF-β pathway, CTGF, SPARC and COL3A1. The uniqueness of this model is that all three genes were included in one model, rather than one gene being examined at a time. A protein might contribute to multiple pathways and this approach allows the identification of important genetic variations linked to multiple genes belonging to the same pathway. In this study, 29 SNPs were identified and 16 of them located in known genes. Exploring the roles of these genes in TGF-β regulation will help elucidate the etiology of SSc, which will in turn help to better manage this complex disease. ^

Human carbonic anhydrase XII: cDNA cloning, expression, and chromosomal localization of a carbonic anhydrase gene that is overexpressed in some renal cell cancers

We report the cloning and characterization of a tumor-associated carbonic anhydrase (CA) that was identified in a human renal cell carcinoma (RCC) by serological expression screening with autologous antibodies. The cDNA sequence predicts a 354-amino acid polypeptide with a molecular mass of 39,448 Da that has features of a type I membrane protein. The predicted sequence includes a 29-amino acid signal sequence, a 261-amino acid CA domain, an additional short extracellular segment, a 26-amino acid hydrophobic transmembrane domain, and a hydrophilic C-terminal cytoplasmic tail of 29 amino acids that contains two potential phosphorylation sites. The extracellular CA domain shows 30–42% homology with known human CAs, contains all three Zn-binding histidine residues found in active CAs, and contains two potential sites for asparagine glycosylation. When expressed in COS cells, the cDNA produced a 43- to 44-kDa protein in membranes that had around one-sixth the CA activity of membranes from COS cells transfected with the same vector expressing bovine CA IV. We have designated this human protein CA XII. Northern blot analysis of normal tissues demonstrated a 4.5-kb transcript only in kidney and intestine. However, in 10% of patients with RCC, the CA XII transcript was expressed at much higher levels in the RCC than in surrounding normal kidney tissue. The CA XII gene was mapped by using fluorescence in situ hybridization to 15q22. CA XII is the second catalytically active membrane CA reported to be overexpressed in certain cancers. Its relationship to oncogenesis and its potential as a clinically useful tumor marker clearly merit further investigation.

Direct sequencing of bacterial and P1 artificial chromosome-nested deletions for identifying position-specific single-nucleotide polymorphisms

A loxP-transposon retrofitting strategy for generating large nested deletions from one end of the insert DNA in bacterial artificial chromosomes and P1 artificial chromosomes was described recently [Chatterjee, P. K. & Coren, J. S. (1997) Nucleic Acids Res. 25, 2205–2212]. In this report, we combine this procedure with direct sequencing of nested-deletion templates by using primers located in the transposon end to illustrate its value for position-specific single-nucleotide polymorphism (SNP) discovery from chosen regions of large insert clones. A simple ampicillin sensitivity screen was developed to facilitate identification and recovery of deletion clones free of transduced transposon plasmid. This directed approach requires minimal DNA sequencing, and no in vitro subclone library generation; positionally oriented SNPs are a consequence of the method. The procedure is used to discover new SNPs as well as physically map those identified from random subcloned libraries or sequence databases. The deletion templates, positioned SNPs, and markers are also used to orient large insert clones into a contig. The deletion clone can serve as a ready resource for future functional genomic studies because each carries a mammalian cell-specific antibiotic resistance gene from the transposon. Furthermore, the technique should be especially applicable to the analysis of genomes for which a full genome sequence or radiation hybrid cell lines are unavailable.

The E5 gene product of rhesus papillomavirus is an activator of endogenous Ras and phosphatidylinositol-3′-kinase in NIH 3T3 cells

We examined the effect of two rhesus papillomavirus 1 (RhPV) oncogenes on cytokine-induced signal transduction pathways leading to the possible activation of Ras protein (p21ras) and phosphatidylinositol kinase. p21ras in both the activated (GTP-bound) and inactivated (GDP-bound) states were quantitated. NIH 3T3 cell lines expressing the RhPV 1 E5 gene or epidermal growth factor receptor cDNA had about a sixfold higher ratio of p21ras-bound GTP to p21ras-bound GDP as compared with parental NIH 3T3 cells or a cell line expressing the RhPV 1 E7 gene under normal culture conditions, yet expressed similar levels of p21ras. Quiescent cells had dramatically reduced levels of activated p21ras, except those containing RhPV 1 E7. Levels were restored by stimulation with epidermal growth factor or platelet-derived growth factor. Both epidermal growth factor and platelet-derived growth factor receptor of RhPV 1 E5- and E7-containing cells responded to cytokine stimulation. Endogenous phosphatidylinositol-3′-kinase was up-regulated in NIH 3T3 cells transformed with the E5 genes of RhPV 1 and bovine papillomavirus 1. These results suggest that E5 genes of papillomaviruses play a major role in the regulation of transduction pathways.

Sequence variation in the Fanconi anemia gene FAA

Fanconi anemia (FA) is a genetically heterogeneous autosomal recessive syndrome associated with chromosomal instability, hypersensitivity to DNA crosslinking agents, and predisposition to malignancy. The gene for FA complementation group A (FAA) recently has been cloned. The cDNA is predicted to encode a polypeptide of 1,455 amino acids, with no homologies to any known protein that might suggest a function for FAA. We have used single-strand conformational polymorphism analysis to screen genomic DNA from a panel of 97 racially and ethnically diverse FA patients from the International Fanconi Anemia Registry for mutations in the FAA gene. A total of 85 variant bands were detected. Forty-five of the variants are probably benign polymorphisms, of which nine are common and can be used for various applications, including mapping studies for other genes in this region of chromosome 16q. Amplification refractory mutation system assays were developed to simplify their detection. Forty variants are likely to be pathogenic mutations. Seventeen of these are microdeletions/microinsertions associated with short direct repeats or homonucleotide tracts, a type of mutation thought to be generated by a mechanism of slipped-strand mispairing during DNA replication. A screening of 350 FA probands from the International Fanconi Anemia Registry for two of these deletions (1115–1118del and 3788–3790del) revealed that they are carried on about 2% and 5% of the FA alleles, respectively. 3788–3790del appears in a variety of ethnic groups and is found on at least two different haplotypes. We suggest that FAA is hypermutable, and that slipped-strand mispairing, a mutational mechanism recognized as important for the generation of germ-line and somatic mutations in a variety of cancer-related genes, including p53, APC, RB1, WT1, and BRCA1, may be a major mechanism for FAA mutagenesis.

Yeast microarrays for genome wide parallel genetic and gene expression analysis

We have developed high-density DNA microarrays of yeast ORFs. These microarrays can monitor hybridization to ORFs for applications such as quantitative differential gene expression analysis and screening for sequence polymorphisms. Automated scripts retrieved sequence information from public databases to locate predicted ORFs and select appropriate primers for amplification. The primers were used to amplify yeast ORFs in 96-well plates, and the resulting products were arrayed using an automated micro arraying device. Arrays containing up to 2,479 yeast ORFs were printed on a single slide. The hybridization of fluorescently labeled samples to the array were detected and quantitated with a laser confocal scanning microscope. Applications of the microarrays are shown for genetic and gene expression analysis at the whole genome level.

Transgenic cattle produced by reverse-transcribed gene transfer in oocytes

A critical requirement for integration of retroviruses, other than HIV and possibly related lentiviruses, is the breakdown of the nuclear envelope during mitosis. Nuclear envelope breakdown occurs during mitotic M-phase, the envelope reforming immediately after cell division, thereby permitting the translocation of the retroviral preintegration complex into the nucleus and enabling integration to proceed. In the oocyte, during metaphase II (MII) of the second meiosis, the nuclear envelope is also absent and the oocyte remains in MII arrest for a much longer period of time compared with M-phase in a somatic cell. Pseudotyped replication-defective retroviral vector was injected into the perivitelline space of bovine oocytes during MII. We show that reverse-transcribed gene transfer can take place in an oocyte in MII arrest of meiosis, leading to production of offspring, the majority of which are transgenic. We discuss the implications of this mechanism both as a means of production of transgenic livestock and as a model for naturally occurring recursive transgenesis.

Scrapie susceptibility-linked polymorphisms modulate the in vitro conversion of sheep prion protein to protease-resistant forms

Prion diseases are natural transmissible neurodegenerative disorders in humans and animals. They are characterized by the accumulation of a protease-resistant scrapie-associated prion protein (PrPSc) of the host-encoded cellular prion protein (PrPC) mainly in the central nervous system. Polymorphisms in the PrP gene are linked to differences in susceptibility for prion diseases. The mechanisms underlying these effects are still unknown. Here we describe studies of the influence of sheep PrP polymorphisms on the conversion of PrPC into protease-resistant forms. In a cell-free system, sheep PrPSc induced the conversion of sheep PrPC into protease-resistant PrP (PrP-res) similar or identical to PrPSc. Polymorphisms present in either PrPC or PrPSc had dramatic effects on the cell-free conversion efficiencies. The PrP variant associated with a high susceptibility to scrapie and short survival times of scrapie-affected sheep was efficiently converted into PrP-res. The wild-type PrP variant associated with a neutral effect on susceptibility and intermediate survival times was converted with intermediate efficiency. The PrP variant associated with scrapie resistance and long survival times was poorly converted. Thus the in vitro conversion characteristics of the sheep PrP variants reflect their linkage with scrapie susceptibility and survival times of scrapie-affected sheep. The modulating effect of the polymorphisms in PrPC and PrPSc on the cell-free conversion characteristics suggests that, besides the species barrier, polymorphism barriers play a significant role in the transmissibility of prion diseases.

Major recent and independent changes in levels and patterns of expression have occurred at the b gene, a regulatory locus in maize

The b locus encodes a transcription factor that regulates the expression of genes that produce purple anthocyanin pigment. Different b alleles are expressed in distinct tissues, causing tissue-specific anthocyanin production. Understanding how phenotypic diversity is produced and maintained at the b locus should provide models for how other regulatory genes, including those that influence morphological traits and development, evolve. We have investigated how different levels and patterns of pigmentation have evolved by determining the phenotypic and evolutionary relationships between 18 alleles that represent the diversity of b alleles in Zea mays. Although most of these alleles have few phenotypic differences, five alleles have very distinct tissue-specific patterns of pigmentation. Superimposing the phenotypes on the molecular phylogeny reveals that the alleles with strong and distinctive patterns of expression are closely related to alleles with weak expression, implying that the distinctive patterns have arisen recently. We have identified apparent insertions in three of the five phenotypically distinct alleles, and the fourth has unique upstream restriction fragment length polymorphisms relative to closely related alleles. The insertion in B-Peru has been shown to be responsible for its unique expression and, in the other two alleles, the presence of the insertion correlates with the phenotype. These results suggest that major changes in gene expression are probably the result of large-scale changes in DNA sequence and/or structure most likely mediated by transposable elements.

Genetic epidemiology of single-nucleotide polymorphisms

On the causal hypothesis, most genetic determinants of disease are single-nucleotide polymorphisms (SNPs) that are likely to be selected as markers for positional cloning. On the proximity hypothesis, most disease determinants will not be included among markers but may be detected through linkage disequilibrium with other SNPs. In that event, allelic association among SNPs is an essential factor in positional cloning. Recent simulation based on monotonic population expansion suggests that useful association does not usually extend beyond 3 kb. This is contradicted by significant disequilibrium at much greater distances, with corresponding reduction in the number of SNPs required for a cost-effective genome scan. A plausible explanation is that cyclical expansions follow population bottlenecks that establish new disequilibria. Data on more than 1,000 locus pairs indicate that most disequilibria trace to the Neolithic, with no apparent difference between haplotypes that are random or selected through a major disease gene. Short duration may be characteristic of alleles contributing to disease susceptibility and haplotypes characteristic of particular ethnic groups. Alleles that are highly polymorphic in all ethnic groups may be older, neutral, or advantageous, in weak disequilibrium with nearby markers, and therefore less useful for positional cloning of disease genes. Significant disequilibrium at large distance makes the number of suitably chosen SNPs required for genome screening as small as 30,000, or 1 per 100 kb, with greater density (including less common SNPs) reserved for candidate regions.

A serotonin transporter gene intron 2 polymorphic region, correlated with affective disorders, has allele-dependent differential enhancer-like properties in the mouse embryo

Polymorphic regions consisting of a variable number of tandem repeats within intron 2 of the gene coding for the serotonin transporter protein 5-HTT have been associated with susceptibility to affective disorders. We have cloned two of these intronic polymorphisms, Stin2.10 and Stin2.12, into an expression vector containing a heterologous minimal promoter and the bacterial LacZ reporter gene. These constructs were then used to produce transgenic mice. In embryonic day 10.5 embryos, both Stin2.10 and Stin2.12 produced consistent β-galactosidase expression in the embryonic midbrain, hindbrain, and spinal cord floor plate. However, we observed that the levels of β-galactosidase expression produced by both the Stin2.10 and Stin2.12 within the rostral hindbrain differed significantly at embryonic day 10.5. Our data suggest that these polymorphic variable number of tandem repeats regions act as transcriptional regulators and have allele-dependent differential enhancer-like properties within an area of the hindbrain where the 5-HTT gene is known to be transcribed at this stage of development.

A family of phase-variable restriction enzymes with differing specificities generated by high-frequency gene rearrangements

The hsd genes of Mycoplasma pulmonis encode restriction and modification enzymes exhibiting a high degree of sequence similarity to the type I enzymes of enteric bacteria. The S subunits of type I systems dictate the DNA sequence specificity of the holoenzyme and are required for both the restriction and the modification reactions. The M. pulmonis chromosome has two hsd loci, both of which contain two hsdS genes each and are complex, site-specific DNA inversion systems. Embedded within the coding region of each hsdS gene are a minimum of three sites at which DNA inversions occur to generate extensive amino acid sequence variations in the predicted S subunits. We show that the polymorphic hsdS genes produced by gene rearrangement encode a family of functional S subunits with differing DNA sequence specificities. In addition to creating polymorphisms in hsdS sequences, DNA inversions regulate the phase-variable production of restriction activity because the other genes required for restriction activity (hsdR and hsdM) are expressed only from loci that are oriented appropriately in the chromosome relative to the hsd promoter. These data cast doubt on the prevailing paradigms that restriction systems are either selfish or function to confer protection from invasion by foreign DNA.

Double muscling in cattle due to mutations in the myostatin gene

Myostatin (GDF-8) is a member of the transforming growth factor β superfamily of secreted growth and differentiation factors that is essential for proper regulation of skeletal muscle mass in mice. Here we report the myostatin sequences of nine other vertebrate species and the identification of mutations in the coding sequence of bovine myostatin in two breeds of double-muscled cattle, Belgian Blue and Piedmontese, which are known to have an increase in muscle mass relative to conventional cattle. The Belgian Blue myostatin sequence contains an 11-nucleotide deletion in the third exon which causes a frameshift that eliminates virtually all of the mature, active region of the molecule. The Piedmontese myostatin sequence contains a missense mutation in exon 3, resulting in a substitution of tyrosine for an invariant cysteine in the mature region of the protein. The similarity in phenotypes of double-muscled cattle and myostatin null mice suggests that myostatin performs the same biological function in these two species and is a potentially useful target for genetic manipulation in other farm animals.