Characterization and mapping of expressed sequence tags isolated from a subtracted cDNA library of Litopenaeus vannamei injected with white spot syndrome virus

A large number of polymorphic simple sequence repeats (SSRs) or microsatellites are needed to develop a genetic map for shrimp. However, developing an SSR map is very time-consuming, expensive, and most SSRs are not specifically linked to gene loci of immediate interest. We report here on our strategy to develop polymorphic markers using expressed sequence tags (ESTs) by designing primers flanking single or multiple SSRs with three or more repeats. A subtracted cDNA library was prepared using RNA from specific pathogen-free (SPF) Litopenaeus vannamei juveniles (similar to 1 g) collected before (0) and after (48 h) inoculation with the China isolate of white spot syndrome virus (WSSV). A total of 224 clones were sequenced, 194 of which were useful for homology comparisons against annotated genes in NCBI nonredundant (nr) and protein databases, providing 179 sequences encoded by nuclear DNA, 4 mitochondrial DNA, and 11 were similar to portions of WSSV genome. The nuclear sequences clustered in 43 groups, 11 of which were homologous to various ESTs of unknown function, 4 had no homology to any sequence, and 28 showed similarities to known genes of invertebrates and vertebrates, representatives of cellular metabolic processes such as calcium ion balance, cytoskeleton mRNAs, and protein synthesis. A few sequences were homologous to immune system-related (allergens) genes and two were similar to motifs of the sex-lethal gene of Drosophila. A large number of EST sequences were similar to domains of the EF-hand superfamily (Ca2+ binding motif and FRQ protein domain of myosin light chains). Single or multiple SSRs with three or more repeats were found in approximately 61 % of the 179 nuclear sequences. Primer sets were designed from 28 sequences representing 19 known or putative genes and tested for polymorphism (EST-SSR marker) in a small test panel containing 16 individuals. Ten (53%) of the 19 putative or unknown function genes were polymorphic, 4 monomorphic, and 3 either failed to satisfactorily amplify genomic DNA or the allele amplification conditions need to be further optimized. Five polymorphic ESTs were genotyped with the entire reference mapping family, two of them (actin, accession #CX535973 and shrimp allergen arginine kinase, accession #CX535999) did not amplify with all offspring of the IRMF panel suggesting presence of null alleles, and three of them amplified in most of the IRM F offspring and were used for linkage analysis. EF-hand motif of myosin light chain (accession #CX535935) was placed in ShrimpMap's linkage group 7, whereas ribosomal protein S5 (accession #CX535957) and troponin I (accession #CX535976) remained unassigned. Results indicate that (a) a large number of ESTs isolated from this cDNA library are similar to cytoskeleton mRNAs and may reflect a normal pathway of the cellular response after im infection with WSSV, and (b) primers flanking single or multiple SSRs with three or more repeats from shrimp ESTs could be an efficient approach to develop polymorphic markers useful for linkage mapping. Work is underway to map additional SSR-containing ESTs from this and other cDNA libraries as a plausible strategy to increase marker density in ShrimpMap.

Genetic mapping of size-related quantitative trait loci (QTL) in the bay scallop (Argopecten irradians) using AFLP and microsatellite markers

We constructed genetic linkage maps for the bay scallop Argopecten irradians using AFLP and microsatellite markers and conducted composite interval mapping (CIM) of body-size-related traits. Three hundred seventeen AFLP and 10 microsatellite markers were used for map construction. The female parent map contained 120 markers in 15 linkage groups, spanning 479.6 cM with an average interval of 7.0 cM. The male parent map had 190 markers in 17 linkage groups, covering 883.8 cM at 7.2 cM per marker. The observed coverage was 70.4% for the female and 81.1% for the male map. Markers that were distorted toward the same direction were closely linked to each other on the genetic maps, suggesting the presence of genes important for survival. Six size-related traits, shell length, shell height, shell width, total weight, soft tissue weight, and shell weight, were measured for QTL mapping. The size data were significantly correlated with each other. We subjected the data, log transformed firstly, to a principle component analysis and use the first principle component for QTL mapping. CIM analysis revealed one significant QTL (LOD=2.69, 1000 permutation, P<0.05) in linkage group 3 on the female parent map. The mapping of size-related QTL in this study raises the possibility of improving the growth of bay scallops through marker-assisted selection. (c) 2007 Published by Elsevier B.V.

Identification of Crassostrea ariakensis and related oysters by multiplex species-specific PCR

Genetic markers are needed for rapid and reliable identification of oysters. In this study, we developed multiplex genus- and species-specific PCR markers for the identification of oysters from China. We used the mitochondrial cytochrome oxidase I (COI) and nuclear 28S ribosomal RNA genes for marker development. DNA sequences from different species were obtained from GenBank or by direct sequencing. Sequences were aligned, and genus- and species-specific nucleotides were identified. Primers were designed for genus/species-specific amplification to generate fragments of different sizes. A multiplex set of genus- and species-specific primers from the 28S gene was able to separate C. ariakensis and C. hongkongensis from other species and assign oysters to four genera. A set of species-specific COI primers provided positive identification of all five Crassostrea species from China, C. ariakensis, C. hongkongensis, C. angulata, C. gigas, and C. sikamea in a single PCR. The multiplex PCR assays do not require fluorescence-labeling or post-PCR enzyme digestion, providing a simple, fast and reliable method for the identification of oysters from China.

The ecological effects of exotic disease resistance genes introgressed into British gooseberries

Warren, J. and James, P. (2006). The ecological effects of exotic disease resistance genes introgressed into British gooseberries. Oecologia 147(1),69-75. RAE2008

Genetic correlates of longevity and selected age-related phenotypes: a genome-wide association study in the Framingham Study

BACKGROUND: Family studies and heritability estimates provide evidence for a genetic contribution to variation in the human life span. METHODS:We conducted a genome wide association study (Affymetrix 100K SNP GeneChip) for longevity-related traits in a community-based sample. We report on 5 longevity and aging traits in up to 1345 Framingham Study participants from 330 families. Multivariable-adjusted residuals were computed using appropriate models (Cox proportional hazards, logistic, or linear regression) and the residuals from these models were used to test for association with qualifying SNPs (70, 987 autosomal SNPs with genotypic call rate [greater than or equal to]80%, minor allele frequency [greater than or equal to]10%, Hardy-Weinberg test p [greater than or equal to] 0.001).RESULTS:In family-based association test (FBAT) models, 8 SNPs in two regions approximately 500 kb apart on chromosome 1 (physical positions 73,091,610 and 73, 527,652) were associated with age at death (p-value < 10-5). The two sets of SNPs were in high linkage disequilibrium (minimum r2 = 0.58). The top 30 SNPs for generalized estimating equation (GEE) tests of association with age at death included rs10507486 (p = 0.0001) and rs4943794 (p = 0.0002), SNPs intronic to FOXO1A, a gene implicated in lifespan extension in animal models. FBAT models identified 7 SNPs and GEE models identified 9 SNPs associated with both age at death and morbidity-free survival at age 65 including rs2374983 near PON1. In the analysis of selected candidate genes, SNP associations (FBAT or GEE p-value < 0.01) were identified for age at death in or near the following genes: FOXO1A, GAPDH, KL, LEPR, PON1, PSEN1, SOD2, and WRN. Top ranked SNP associations in the GEE model for age at natural menopause included rs6910534 (p = 0.00003) near FOXO3a and rs3751591 (p = 0.00006) in CYP19A1. Results of all longevity phenotype-genotype associations for all autosomal SNPs are web posted at http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007. CONCLUSION: Longevity and aging traits are associated with SNPs on the Affymetrix 100K GeneChip. None of the associations achieved genome-wide significance. These data generate hypotheses and serve as a resource for replication as more genes and biologic pathways are proposed as contributing to longevity and healthy aging.

Polymorphisms in the ACE and ADRB2 genes and risks of aging-associated phenotypes: the case of myocardial infarction.

Multiple functions of the beta2-adrenergic receptor (ADRB2) and angiotensin-converting enzyme (ACE) genes warrant studies of their associations with aging-related phenotypes. We focus on multimarker analyses and analyses of the effects of compound genotypes of two polymorphisms in the ADRB2 gene, rs1042713 and rs1042714, and 11 polymorphisms of the ACE gene, on the risk of such an aging-associated phenotype as myocardial infarction (MI). We used the data from a genotyped sample of the Framingham Heart Study Offspring (FHSO) cohort (n = 1500) followed for about 36 years with six examinations. The ADRB2 rs1042714 (C-->G) polymorphism and two moderately correlated (r(2) = 0.77) ACE polymorphisms, rs4363 (A-->G) and rs12449782 (A-->G), were significantly associated with risks of MI in this aging cohort in multimarker models. Predominantly linked ACE genotypes exhibited opposite effects on MI risks, e.g., the AA (rs12449782) genotype had a detrimental effect, whereas the predominantly linked AA (rs4363) genotype exhibited a protective effect. This trade-off occurs as a result of the opposite effects of rare compound genotypes of the ACE polymorphisms with a single dose of the AG heterozygote. This genetic trade-off is further augmented by the selective modulating effect of the rs1042714 ADRB2 polymorphism. The associations were not altered by adjustment for common MI risk factors. The results suggest that effects of single specific genetic variants of the ADRB2 and ACE genes on MI can be readily altered by gene-gene or/and gene-environmental interactions, especially in large heterogeneous samples. Multimarker genetic analyses should benefit studies of complex aging-associated phenotypes.

Evolution of the sex-related locus and genomic features shared in microsporidia and fungi.

BACKGROUND: Microsporidia are obligate intracellular, eukaryotic pathogens that infect a wide range of animals from nematodes to humans, and in some cases, protists. The preponderance of evidence as to the origin of the microsporidia reveals a close relationship with the fungi, either within the kingdom or as a sister group to it. Recent phylogenetic studies and gene order analysis suggest that microsporidia share a particularly close evolutionary relationship with the zygomycetes. METHODOLOGY/PRINCIPAL FINDINGS: Here we expanded this analysis and also examined a putative sex-locus for variability between microsporidian populations. Whole genome inspection reveals a unique syntenic gene pair (RPS9-RPL21) present in the vast majority of fungi and the microsporidians but not in other eukaryotic lineages. Two other unique gene fusions (glutamyl-prolyl tRNA synthetase and ubiquitin-ribosomal subunit S30) that are present in metazoans, choanoflagellates, and filasterean opisthokonts are unfused in the fungi and microsporidians. One locus previously found to be conserved in many microsporidian genomes is similar to the sex locus of zygomycetes in gene order and architecture. Both sex-related and sex loci harbor TPT, HMG, and RNA helicase genes forming a syntenic gene cluster. We sequenced and analyzed the sex-related locus in 11 different Encephalitozoon cuniculi isolates and the sibling species E. intestinalis (3 isolates) and E. hellem (1 isolate). There was no evidence for an idiomorphic sex-related locus in this Encephalitozoon species sample. According to sequence-based phylogenetic analyses, the TPT and RNA helicase genes flanking the HMG genes are paralogous rather than orthologous between zygomycetes and microsporidians. CONCLUSION/SIGNIFICANCE: The unique genomic hallmarks between microsporidia and fungi are independent of sequence based phylogenetic comparisons and further contribute to define the borders of the fungal kingdom and support the classification of microsporidia as unusual derived fungi. And the sex/sex-related loci appear to have been subject to frequent gene conversion and translocations in microsporidia and zygomycetes.

Chromosomal organization of adrenergic receptor genes.

The adrenergic receptors (ARs) (subtypes alpha 1, alpha 2, beta 1, and beta 2) are a prototypic family of guanine nucleotide binding regulatory protein-coupled receptors that mediate the physiological effects of the hormone epinephrine and the neurotransmitter norepinephrine. We have previously assigned the genes for beta 2- and alpha 2-AR to human chromosomes 5 and 10, respectively. By Southern analysis of somatic cell hybrids and in situ chromosomal hybridization, we have now mapped the alpha 1-AR gene to chromosome 5q32----q34, the same position as beta 2-AR, and the beta 1-AR gene to chromosome 10q24----q26, the region where alpha 2-AR is located. In mouse, both alpha 2- and beta 1-AR genes were assigned to chromosome 19, and the alpha 1-AR locus was localized to chromosome 11. Pulsed field gel electrophoresis has shown that the alpha 1- and beta 2-AR genes in humans are within 300 kilobases (kb) and the distance between the alpha 2- and beta 1-AR genes is less than 225 kb. The proximity of these two pairs of AR genes and the sequence similarity that exists among all the ARs strongly suggest that they are evolutionarily related. Moreover, they likely arose from a common ancestral receptor gene and subsequently diverged through gene duplication and chromosomal duplication to perform their distinctive roles in mediating the physiological effects of catecholamines. The AR genes thus provide a paradigm for understanding the evolution of such structurally conserved yet functionally divergent families of receptor molecules.

Association between DNA damage response and repair genes and risk of invasive serous ovarian cancer.

BACKGROUND: We analyzed the association between 53 genes related to DNA repair and p53-mediated damage response and serous ovarian cancer risk using case-control data from the North Carolina Ovarian Cancer Study (NCOCS), a population-based, case-control study. METHODS/PRINCIPAL FINDINGS: The analysis was restricted to 364 invasive serous ovarian cancer cases and 761 controls of white, non-Hispanic race. Statistical analysis was two staged: a screen using marginal Bayes factors (BFs) for 484 SNPs and a modeling stage in which we calculated multivariate adjusted posterior probabilities of association for 77 SNPs that passed the screen. These probabilities were conditional on subject age at diagnosis/interview, batch, a DNA quality metric and genotypes of other SNPs and allowed for uncertainty in the genetic parameterizations of the SNPs and number of associated SNPs. Six SNPs had Bayes factors greater than 10 in favor of an association with invasive serous ovarian cancer. These included rs5762746 (median OR(odds ratio)(per allele) = 0.66; 95% credible interval (CI) = 0.44-1.00) and rs6005835 (median OR(per allele) = 0.69; 95% CI = 0.53-0.91) in CHEK2, rs2078486 (median OR(per allele) = 1.65; 95% CI = 1.21-2.25) and rs12951053 (median OR(per allele) = 1.65; 95% CI = 1.20-2.26) in TP53, rs411697 (median OR (rare homozygote) = 0.53; 95% CI = 0.35 - 0.79) in BACH1 and rs10131 (median OR( rare homozygote) = not estimable) in LIG4. The six most highly associated SNPs are either predicted to be functionally significant or are in LD with such a variant. The variants in TP53 were confirmed to be associated in a large follow-up study. CONCLUSIONS/SIGNIFICANCE: Based on our findings, further follow-up of the DNA repair and response pathways in a larger dataset is warranted to confirm these results.

Two Antarctic penguin genomes reveal insights into their evolutionary history and molecular changes related to the Antarctic environment.

BACKGROUND: Penguins are flightless aquatic birds widely distributed in the Southern Hemisphere. The distinctive morphological and physiological features of penguins allow them to live an aquatic life, and some of them have successfully adapted to the hostile environments in Antarctica. To study the phylogenetic and population history of penguins and the molecular basis of their adaptations to Antarctica, we sequenced the genomes of the two Antarctic dwelling penguin species, the Adélie penguin [Pygoscelis adeliae] and emperor penguin [Aptenodytes forsteri]. RESULTS: Phylogenetic dating suggests that early penguins arose ~60 million years ago, coinciding with a period of global warming. Analysis of effective population sizes reveals that the two penguin species experienced population expansions from ~1 million years ago to ~100 thousand years ago, but responded differently to the climatic cooling of the last glacial period. Comparative genomic analyses with other available avian genomes identified molecular changes in genes related to epidermal structure, phototransduction, lipid metabolism, and forelimb morphology. CONCLUSIONS: Our sequencing and initial analyses of the first two penguin genomes provide insights into the timing of penguin origin, fluctuations in effective population sizes of the two penguin species over the past 10 million years, and the potential associations between these biological patterns and global climate change. The molecular changes compared with other avian genomes reflect both shared and diverse adaptations of the two penguin species to the Antarctic environment.

Inhibitors of alpha-synuclein oligomerization and toxicity: a future therapeutic strategy for Parkinson's disease and related disorders.

An abundance of genetic, histopathological, and biochemical evidence has implicated the neuronal protein, alpha-synuclein (alpha-syn) as a key player in the development of several neurodegenerative diseases, the so-called synucleinopathies, of which Parkinson's disease (PD) is the most prevalent. Development of disease appears to be linked to events that increase the intracellular concentration of alpha-syn or cause its chemical modification, either of which can accelerate the rate at which it forms aggregates. Examples of such events include increased copy number of genes, decreased rate of degradation via the proteasome or other proteases, or altered forms of alpha-syn, such as truncations, missense mutations, or chemical modifications by oxidative reactions. Aggregated forms of the protein, especially newly formed soluble aggregates, are toxic to cells, so that one therapeutic strategy would be to reduce the rate at which such oligomerization occurs. We have therefore designed several peptides and also identified small molecules that can inhibit alpha-syn oligomerization and toxicity in vitro. These compounds could serve as lead compounds for the design of new drugs for the treatment of PD and related disorders in the future.

A novel diagnostic test detects a low frequency of the hemicentin Gln5345Arg variant among Northern Irish age related macular degeneration patients.

Purpose: Age related macular degeneration (AMD) is a common cause of severe vision loss. Identification of genes involved in AMD will facilitate early detection and ultimately help to identify pathways for treatment for this disorder. The A16,263G mutation in the HEMICENTIN-1 gene produces a non-conservative substitution of arginine for glutamine at codon 5345 which has been implicated in familial AMD. The aim of this study is to develop a rapid diagnostic assay for the detection of this mutation and to evaluate its frequency in a sample of AMD patients. Methods: A primer probe set was designed from exon 104 of the HEMICENTIN-1 gene to differentiate between mutant and wild type alleles. A region spanning the mutation was amplified by PCR using a LightCycler (Roche Diagnostic). The mutation was then detected by melt curve analysis of the hybrid formed between the PCR product and a specific fluorescent probe. The frequency of the mutation within the Northern Ireland population was evaluated by assaying 508 affected AMD patients, 25 possibly affected and 163 controls. Results: This assay clearly discriminates between the A16,263G mutant and wild type HEMICENTIN-1 alleles. The wild type sequence has a single base mismatch with the probe which decreases the stability of the hybrid, resulting in a lower TM (TM=51.27 °C) than that observed for the perfectly matched mutant allele (TM=59.9 °C). The mutant allele was detected in only one of the 696 subjects, an affected AMD patient. Conclusions: We describe a rapid assay for the genotyping of the Gln5345Arg mutation using real-time fluorescence PCR to facilitate rapid processing of samples through combined amplification and detection steps. These characteristics are suitable for a clinical setting where high throughput diagnostic procedures are required. The frequency of this mutation within the Northern Ireland population has been estimated at 0.2%, concurring with previous findings that this mutation is a rare variant associated with AMD. A rapid diagnostic assay will facilitate a reliable and convenient evaluation of the frequency of the Gln5345Arg mutation and its association with AMD within other populations.

The complex array of bradykinin-related peptides (BRPs) in the peptidome of pickerel frog (Rana palustris) skin secretion is the product of transcriptional economy.

Previous peptidomic analyses of the defensive skin secretion from the North American pickerel frog, Rana palustris, have established the presence of canonical bradykinin and multiple bradykinin-related peptides (BRPs). As a consequence of the multiplicity of peptides identified and their diverse primary structures, it was speculated that they must represent the products of expression of multiple genes. Here, we present unequivocal evidence that the majority of BRPs (11/13) identified in skin secretion by the peptidomic approach can be generated by differential site-specific protease cleavage from a single common precursor of 321 amino acid residues, named skin kininogen 1, whose primary structure was deduced from cloned skin secretion-derived cDNA. The organization of skin kininogen 1 consists of a hydrophobic signal peptide followed by eight non-identical domains each encoding a single copy of either canonical bradykinin or a BRP. Two additional splice variants, encoding precursors of 233 (skin kininogen 2) or 189 amino acid residues (skin kininogen 3), were also cloned and were found to lack BRP-encoding domains 5 and 6 or 4, 5 and 6, respectively. Thus, generation of peptidome diversity in amphibian defensive skin secretions can be achieved in part by differential protease cleavage of relatively large and multiple-encoding domain precursors reflecting a high degree of transcriptional economy.

Further assessment of the complement component 2 and factor B region associated with age-related macular degeneration

PURPOSE. Polymorphic variation in genes involved in regulation of the complement system has been implicated as a major cause of genetic risk, in addition to the LOC387715/HTRA1 locus and other environmental influences. Previous studies have identified polymorphisms in the complement component 2 (CC2) and factor B (CFB) genes, as potential functional variants associated with AMD, in particular CFB R32Q and CC2 rs547154, both of which share strong linkage disequilibrium (LD). METHODS. Data derived from the HapMap Project were used to select 18 haplotype-tagging SNPs across the extended CC2/ CFB region for genotyping, to measure the strength of LD in 318 patients with neovascular AMD and 243 age-matched control subjects to identify additional potential functional variants in addition to those originally reported. RESULTS. Strong LD was measured across this region as far as the superkiller viralicidic activity 2-like gene (SKIV2L). Nine SNPs were identified to be significantly associated with the genetic effect observed at this locus. Of these, a nonsynonymous coding variant SKIV2L R151Q (rs438999; OR, 0.48; 95% confidence interval [CI], 0.31- 0.74; P < 0.001), was in strong LD with CFB R32Q, rs641153 (r2 = 0.95) and may exert a functional effect. When assessed within a logistic regression model measuring the effects of genetic variation at the CFH and LOC387715/HTRA1 loci and smoking, the effect remained significant (OR, 0.38; 95% CI, 0.22- 0.65; P < 0.001). Additional variation identified within this region may also confer a weaker but independent effect and implicate additional genes within the pathogenesis of AMD. CONCLUSIONS. Because of the high level of LD within the extended CC2/CFB region, variation within SKIV2L may exert a functional effect in AMD. Copyright © Association for Research in Vision and Ophthalmology.

VEGF-induced retinal angiogenic signalling is critically dependent on Ca2+ signalling via Ca2+/calmodulin-dependent protein kinase II

Purpose. The authors conducted an in vitro investigation of the role of Ca2+-dependent signaling in vascular endothelial growth factor (VEGF)-induced angiogenesis in the retina.

Methods. Bovine retinal endothelial cells (BRECs) were stimulated with VEGF in the presence or absence of 1,2-bis(2-aminophenoxy)ethane-N,N,N',N'-tetraacetic acid-acetoxymethyl ester (BAPTA-AM; intracellular Ca2+ chelator), U73122 (phospholipase C (PLC) inhibitor), xestospongin C (Xe-C), and 2-aminoethoxydiphenyl borate (2APB) (inhibitors of inositol-1,4,5 triphosphate (IP3) signaling). Intracellular Ca2+ concentration ([Ca2+]i) was estimated using fura-2 Ca2+ microfluorometry, Akt phosphorylation quantified by Western blot analysis, and angiogenic responses assessed using cell migration, proliferation, tubulogenesis, and sprout formation assays. The effects of the Ca2+/calmodulin-dependent protein kinase II (CaMKII) inhibitor KN93 were also evaluated on VEGF-induced Akt signaling and angiogenic activity.

Results. Stimulation of BRECs with 25 ng/mL VEGF induced a biphasic increase in [Ca2+]i, with an initial transient peak followed by a sustained plateau phase. VEGF-induced [Ca2+]i increases were almost completely abolished by pretreating the cells with BAPTA-AM, U73122, Xe-C, or 2APB. These agents also inhibited VEGF-induced phosphorylation of Akt, cell migration, proliferation, tubulogenesis, and sprouting angiogenesis. KN93 was similarly effective at blocking the VEGF-induced activation of Akt and angiogenic responses.

Conclusions. VEGF increases [Ca2+]i in BRECs through activation of the PLC-IP3 signal transduction pathway. VEGF-induced phosphorylation of the proangiogenic protein Akt is critically dependent on this increase in [Ca2+]i and the subsequent activation of CaMKII. Pharmacologic inhibition of Ca2+-mediated signaling in retinal endothelial cells blocks VEGF-induced angiogenic responses. These results suggest that the PLC/IP3/Ca2+/CaMKII signaling pathway may be a rational target for the treatment of angiogenesis-related disorders of the eye.