Evolutionary Genomics and Conservation of the Endangered Przewalski's Horse.

Przewalski's horses (PHs, Equus ferus ssp. przewalskii) were discovered in the Asian steppes in the 1870s and represent the last remaining true wild horses. PHs became extinct in the wild in the 1960s but survived in captivity, thanks to major conservation efforts. The current population is still endangered, with just 2,109 individuals, one-quarter of which are in Chinese and Mongolian reintroduction reserves [1]. These horses descend from a founding population of 12 wild-caught PHs and possibly up to four domesticated individuals [2-4]. With a stocky build, an erect mane, and stripped and short legs, they are phenotypically and behaviorally distinct from domesticated horses (DHs, Equus caballus). Here, we sequenced the complete genomes of 11 PHs, representing all founding lineages, and five historical specimens dated to 1878-1929 CE, including the Holotype. These were compared to the hitherto-most-extensive genome dataset characterized for horses, comprising 21 new genomes. We found that loci showing the most genetic differentiation with DHs were enriched in genes involved in metabolism, cardiac disorders, muscle contraction, reproduction, behavior, and signaling pathways. We also show that DH and PH populations split ∼45,000 years ago and have remained connected by gene-flow thereafter. Finally, we monitor the genomic impact of ∼110 years of captivity, revealing reduced heterozygosity, increased inbreeding, and variable introgression of domestic alleles, ranging from non-detectable to as much as 31.1%. This, together with the identification of ancestry informative markers and corrections to the International Studbook, establishes a framework for evaluating the persistence of genetic variation in future reintroduced populations.

Thrombophilia and risk of VTE recurrence according to the age at the time of first VTE manifestation

BACKGROUND Whether screening for thrombophilia is useful for patients after a first episode of venous thromboembolism (VTE) is a controversial issue. However, the impact of thrombophilia on the risk of recurrence may vary depending on the patient's age at the time of the first VTE. PATIENTS AND METHODS Of 1221 VTE patients (42 % males) registered in the MAISTHRO (MAin-ISar-THROmbosis) registry, 261 experienced VTE recurrence during a 5-year follow-up after the discontinuation of anticoagulant therapy. RESULTS Thrombophilia was more common among patients with VTE recurrence than those without (58.6 % vs. 50.3 %; p = 0.017). Stratifying patients by the age at the time of their initial VTE, Cox proportional hazards analyses adjusted for age, sex and the presence or absence of established risk factors revealed a heterozygous prothrombin (PT) G20210A mutation (hazard ratio (HR) 2.65; 95 %-confidence interval (CI) 1.71 - 4.12; p < 0.001), homozygosity/double heterozygosity for the factor V Leiden and/or PT mutation (HR 2.35; 95 %-CI 1.09 - 5.07, p = 0.030), and an antithrombin deficiency (HR 2.12; 95 %-CI 1.12 - 4.10; p = 0.021) to predict recurrent VTE in patients aged 40 years or older, whereas lupus anticoagulants (HR 3.05; 95%-CI 1.40 - 6.66; p = 0.005) increased the risk of recurrence in younger patients. Subgroup analyses revealed an increased risk of recurrence for a heterozygous factor V Leiden mutation only in young females without hormonal treatment whereas the predictive value of a heterozygous PT mutation was restricted to males over the age of 40 years. CONCLUSIONS Our data do not support a preference of younger patients for thrombophilia testing after a first venous thromboembolic event.

The Vacuolar H+-ATPase B1 Subunit Polymorphism p.E161K Associates with Impaired Urinary Acidification in Recurrent Stone Formers

Mutations in the vacuolar–type H+-ATPase B1 subunit gene ATP6V1B1 cause autosomal–recessive distal renal tubular acidosis (dRTA). We previously identified a single-nucleotide polymorphism (SNP) in the human B1 subunit (c.481G.A; p.E161K) that causes greatly diminished pump function in vitro. To investigate the effect of this SNP on urinary acidification, we conducted a genotype-phenotype analysis of recurrent stone formers in theDallas and Bern kidney stone registries. Of 555 patients examined, 32 (5.8%) were heterozygous for the p.E161K SNP, and the remaining 523 (94.2%) carried two wild–type alleles. After adjustment for sex, age, body mass index, and dietary acid and alkali intake, p.E161K SNP carriers had a nonsignificant tendency to higher urinary pH on a random diet (6.31 versus 6.09; P=0.09). Under an instructed low–Ca and low–Na diet, urinary pH was higher in p.E161K SNP carriers (6.56 versus 6.01; P,0.01). Kidney stones of p.E161K carriers were more likely to contain calcium phosphate than stones of wild-type patients. In acute NH4Cl loading, p.E161K carriers displayed a higher trough urinary pH (5.34 versus 4.89; P=0.01) than wild-type patients. Overall, 14.6% of wild-type patients and 52.4% of p.E161K carriers were unable to acidify their urine below pH 5.3 and thus, can be considered to have incomplete dRTA. In summary, our data indicate that recurrent stone formers with the vacuolar H+-ATPase B1 subunit p.E161K SNP exhibit a urinary acidification deficit with an increased prevalence of calcium phosphate– containing kidney stones. The burden of E161K heterozygosity may be a forme fruste of dRTA.

A significant effect of the KIR ligand HLA-C on fibrosis progression in chronic C hepatitis with or without liver transplantation.

BACKGROUND & AIMS The interaction of KIR with their HLA ligands drives the activation and inhibition of natural killer (NK) cells. NK cells could be implicated in the development of liver fibrosis in chronic hepatitis C. METHODS We analysed 206 non-transplanted and 53 liver transplanted patients, selected according to their Metavir fibrosis stage. Several variables such as the number of activator KIR or the HLA ligands were considered in multinomial and logistic regression models. Possible confounding variables were also investigated. RESULTS The KIRs were not significant predictors of the fibrosis stage. Conversely, a significant reduction of the HLA-C1C2 genotype was observed in the most advanced fibrosis stage group (F4) in both cohorts. Furthermore, the progression rate of fibrosis was almost 10 times faster in the subgroup of patients after liver transplantation and HLA-C1C2 was significantly reduced in this cohort compared to non-transplanted patients. CONCLUSION This study suggests a possible role of KIR and their ligands in the development of liver damage. The absence of C1 and C2 ligands heterozygosity could lead to less inhibition of NK cells and a quicker progression to a high level of fibrosis in patients infected by HCV, especially following liver transplantation. This article is protected by copyright. All rights reserved.

IGHD II: A Novel GH-1 Gene Mutation (GH-L76P) Severely Affects GH Folding, Stability, and Secretion

CONTEXT The autosomal dominant form of GH deficiency (IGHD II) is characterized by markedly reduced GH secretion combined with low concentrations of IGF-1 leading to short stature. OBJECTIVE Structure-function analysis of a missense mutation in the GH-1 gene converting codon 76 from leucine (L) to proline (P) yielding a mutant GH-L76P peptide. DESIGN, SETTINGS, AND PATIENTS Heterozygosity for GH-L76P/wt-GH was identified in a nonconsanguineous Spanish family. The index patients, two siblings, a boy and a girl, were referred for assessment of their short stature (-3.2 and -3.8 SD). Their grandmother, father, and aunt were also carrying the same mutation and showed severe short stature; therefore, IGHD II was diagnosed. INTERVENTIONS AND RESULTS AtT-20 cells coexpressing both wt-GH and GH-L76P showed a reduced GH secretion (P < .001) after forskolin stimulation when compared with the cells expressing only wt-GH. In silico mutagenesis and molecular dynamics simulations presented alterations of correct folding and mutant stability compared with wt-GH. Therefore, further structural analysis of the GH-L76P mutant was performed using expressed and purified proteins in Escherichia coli by thermofluor assay and fast degradation proteolysis assay. Both assays revealed that the GH-L76P mutant is unstable and misfolded compared to wt-GH confirming the bioinformatic model prediction. CONCLUSIONS This is the first report of a family suffering from short stature caused by IGHD II, which severely affects intracellular GH folding and stability as well as secretion, highlighting the necessity of functional analysis of any GH variant for defining new mechanisms as a cause for IGHD II.

Expansion load: recessive mutations and the role of standing genetic variation

Expanding populations incur a mutation burden – the so-called expansion load. Previous studies of expansion load have focused on codominant mutations. An important consequence of this assumption is that expansion load stems exclusively from the accumulation of new mutations occurring in individuals living at the wave front. Using individual-based simulations, we study here the dynamics of standing genetic variation at the front of expansions, and its consequences on mean fitness if mutations are recessive. We find that deleterious genetic diversity is quickly lost at the front of the expansion, but the loss of deleterious mutations at some loci is compensated by an increase of their frequencies at other loci. The frequency of deleterious homozygotes therefore increases along the expansion axis, whereas the average number of deleterious mutations per individual remains nearly constant across the species range. This reveals two important differences to codominant models: (i) mean fitness at the front of the expansion drops much faster if mutations are recessive, and (ii) mutation load can increase during the expansion even if the total number of deleterious mutations per individual remains constant. We use our model to make predictions about the shape of the site frequency spectrum at the front of range expansion, and about correlations between heterozygosity and fitness in different parts of the species range. Importantly, these predictions provide opportunities to empirically validate our theoretical results. We discuss our findings in the light of recent results on the distribution of deleterious genetic variation across human populations and link them to empirical results on the correlation of heterozygosity and fitness found in many natural range expansions.

Implications of Longterm Diameter-Limit Harvesting: Effects on Radial Growth of Red Spruce (Picea rubens) and Genetic Diversity of White Pine (Pinus strobus)

For over 3 centuries, diameter-limit harvesting has been a predominant logging method in the northeastern United States. Silvicultural theory asserts that such intensively selective harvesting can lead to genetic degradation. A decrease in softwood productivity has recently been reported in Maine - has a long history of dysgenic selection degraded the genetic resources of Maine softwoods, contributing to a decrease in growth and productivity? This study examines two aspects of potential implications of diameter-limit harvesting: effects on residual phenotypes of red spruce and impacts on genetic diversity of white pine. Radial growth of residual red spruce trees in stands experiencing 50 years of fixed diameter-limit harvesting was measured using annual increment rings and compared with residual red spruce trees in positive selection stands. Trees remaiaing after several rounds of diameter-limit harvesting exhibited sigdicantl y smaller radial sizes throughout their lives, and displayed significantly slower growth rates for the first 80 years of measured growth. These results strongly suggest that the largest and fastest-growing genotypes and their respective gene complexes determining good radial growth have been removed from the diameter-limit stand. Dysgenic selection can be observed in fixed diarneter-limit stands, resulting in a diminished genetic resource and decreased residual stand value. To examine more direct genetic implications of long-term diameter-limit harvesting, microsatellite DNA markers were implemented to study genetic diversity of eastern white pine in Maine. Three age groups of trees were studied: mature trees older than 200 years, juvenile trees 5-30 years old, and embryos. Trees were genotyped at 10 microsatellite loci. Overall genetic diversity levels of eastern white pine in Maine were extremely high, with an average observed heterozygosity of 0.762. Genetic differentiation was minimal among and between all three age groups, although an excess of heterozygotes was shown in the mature and juvenile groups that was not reflected in the embryo group, which actually had a slight heterozygote deficiency. Allele frequencies did not differ significantly between age groups, but did reveal more rare and low frequency alleles in the embryo groups than in the mature group. Overall, low frequency alleles comprise the largest portion of alleles in the sample population, with no common alleles evident overall. These results suggest that significant genetic degradation has either not occurred for white pine, or that the results of dysgenic selection have not yet emerged. It is clear, however, that selective harvesting could result in a loss of low frequency alleles, which are a primary reserve of evolutionary potential in a species. Implications of these studies affect industrial forestry, regional economics, and ecological concerns for the northeast. Long-term diameter-limit harvesting can lead to a degradation of residual phenotypes, and an overall decrease in stand quality. Potentially, a loss of low frequency, locally adapted alleles could result in a decrease of allelic richness and degradation of the regidnal genetic resource. Decreased genetic variation can lead to seriously limited evolutionary potential of species and ecosystems, particularly in rapidly changing environments. Based on these findings, I recommend a reassessment of any harvesting prescription that includes fixed diameter-limit removals, particularly for species that have low natural genetic diversity levels or a limited natural range, such as red spruce. Maintenance of a healthy genetic reserve can avoid effects of dysgenic harvesting.

The localization and identification of candidate tumor suppressor genes involved in prostate cancer

Frequent loss of heterozygosity (LOH) at specific chromosomal regions are highly associated with the inactivation of tumor suppressor genes (TSGs) (Weinberg, 1991; Bishop, 1989). Chromosome 8p is the most frequently reported site of LOH (∼60%) in prostate cancer (PC), suggesting that there may be inactivated TSG(s) involved in PC on chromosome 8p. (Bergerheim et. al., 1991; Kagan et. al., 1995). In order to identify the smallest common regions of frequent LOH (SCLs) on chromosome 8, we screened 52 PC patient/tumor samples with 39 polymorphic markers in successive screenings. In the course of refining the SCLs, we identified 3 tumors with >6 Mb homozygous deletions (HZDs) at 8p22 and 8p21, suggesting the presence of candidate TSGs at both loci. These HZDs spanned the two SCLs at 8p22 (46%) and 8p21 (45%). The SCLs were narrowed to 3.2 cM at 8p22 and less than 3 cM at 8p21. ^ In order to identify candidate TSGs within the SCLs on 8p, two approaches were used. In the candidate gene approach, thirty genes that mapped to the SCLs were evaluated for expression in normal prostate and in PC cell lines. One of the candidate genes, Clusterin, showed decreased expression in 4/7 (57%) prostate cancer cell lines by Northern blot analysis. Clusterin will be further examined as a candidate TSG. ^ The second approach involved utilizing subtractive hybridization and hybrid affinity capture to generate pools of expressed sequence tags (ESTs) enriched for genes that are downregulated or deleted in PC and that map to specific regions of interest. We took advantage of a prostate cancer cell line (PC3) with a known HZD of a candidate TSG, CTNNA1 on 5q31, to develop and validate a model system. We then developed subtracted libraries enriched for 8p22 and 8p21 ESTs by this method, using two cell lines, MDAPCa-2b and PC3. The ESTs were cloned, and 40 were sequenced and evaluated for expression in normal prostate and PC cell lines. Three ESTs from the subtracted libraries, C2, C17 and F12, showed decreased expression in 29–57% of the prostate tumor cell lines studied, and will be further examined as candidate TSGs. ^

Evaluation of candidate genes for susceptibility to spina bifida

Neural tube defects (NTDs) are the most common severely disabling birth defects in the United States, with a frequency of approximately 1–2 of every 1,000 births. This text includes the identification and evaluation of candidate susceptibility genes that confer risk for the development of neural tube defects (NTDs). The project focused on isolated meningomyelocele, also termed spina bifida (SB). ^ Spina bifida is a complex disease with multifactorial inheritance, therefore the subject population (consisting of North American Caucasians and Hispanics of Mexicali-American descent) was composed of 459 simplex SB families who were tested for genetic associations utilizing the transmission disequilibrium test (TDT), a nonparametric linkage technique. Three categories of candidate genes were studied, including (1) human equivalents of genes determined in mouse models to cause NTDs, (2) HOX and PAX genes, and (3) the MTHFR gene involved in the metabolic pathway of folate. ^ The C677T variant of the 5,10-methylenetetrahydrofolate reductase (MTHFR) gene was the first mutation in this gene to be implicated as a risk factor for NTDs. Our evaluation of the MTHFR gene provides evidence that maternal C677T homozygosity is a risk factor for upper level spina bifida defects in Hispanics [OR = 2.3, P = 0.02]. This observed risk factor is of great importance due to the high prevalence of this homozygous genotype in the Hispanic population. Additionally, maternal C677T/A1298C compound heterozygosity is a risk factor for upper level spina bifida defects in non-Hispanic whites [OR = 3.6, P = 0.03]. ^ For TDT analysis, our total population of 1128 subjects were genotyped for 54 markers from within and/or flanking the 20 candidate genes/gene regions of interest. Significant TDT findings were obtained for 3 of the 54 analyzed markers: d20s101 flanking the PAX1 gene (P = 0.019), d1s228 within the PAX7 gene (P = 0.011), and d2s110 within the PAX8 gene (P = 0.013). These results were followed-up by testing the genes directly for mutations utilizing single-strand conformational analysis (SSCA) and direct sequencing. Multiple variations were detected in each of these PAX genes; however, these variations were not passed from parent to child in phase with the positively transmitted alleles. Therefore, these variations do not contribute to the susceptibility of spina bifida, but rather are previously unreported single nucleotide polymorphisms. ^

Physical and functional mapping of a novel tumor suppressor locus for prostate cancer within chromosome 10p12.31-q11

Prostate cancer remains the second leading cause of male cancer deaths in the United States, yet the molecular mechanisms underlying this disease remain largely unknown. Cytogenetic and molecular analyses of prostate tumors suggest a consistent association with the loss of chromosome 10. Previously, we have defined a novel tumor suppressor locus PAC-1 within chromosome 10pter-q11. Introduction of the short arm of chromosome 10 into a prostatic adenocarcinoma cell line PC-3H resulted in dramatic tumor suppression and restoration of a programmed cell death pathway. Using a combined approach of comparative genomic hybridization and microsatellite analysis of PC-3H, I have identified a region of hemizygosity within 10p12-p15. This region has been shown to be involved in frequent loss of heterozygosity in gliomas and melanoma. To functionally dissect the region within chromosome 10p containing PAC-1, we developed a strategy of serial microcell fusion, a technique that allows the transfer of defined fragments of chromosome 10p into PC-3H. Serial microcell fusion was used to transfer defined 10p fragments into a mouse A9 fibrosarcoma cell line. Once characterized by FISH and microsatellite analyses, the 10p fragments were subsequently transferred into PC-3H to generate a panel of microcell hybrid clones containing overlapping deletions of chromosome 10p. In vivo and microsatellite analyses of these PC hybrids identified a small chromosome 10p fragment (an estimated 31 Mb in size inclusive of the centromere) that when transferred into the PC-3H background, resulted in significant tumor suppression and limited a region of functional tumor suppressor activity to chromosome 10p12.31-q11. This region coincides with a region of LOH demonstrated in prostate cancer. These studies demonstrate the utility of this approach as a powerful tool to limit regions of functional tumor suppressor activity. Furthermore, these data used in conjunction with data generated by the Human Genome Project lent a focused approach to identify candidate tumor suppressor genes involved in prostate cancer. ^

Towards the identification of a tumor suppressor gene on chromosome 3p12: Physical mapping and candidate gene screening

Renal cell carcinoma (RCC) is the most common malignant tumor of the kidney. Characterization of RCC tumors indicates that the most frequent genetic event associated with the initiation of tumor formation involves a loss of heterozygosity or cytogenetic aberration on the short arm of human chromosome 3. A tumor suppressor locus Nonpapillary Renal Carcinoma-1 (NRC-1, OMIM ID 604442) has been previously mapped to a 5–7 cM region on chromosome 3p12 and shown to induce rapid tumor cell death in vivo, as demonstrated by functional complementation experiments. ^ To identify the gene that accounts for the tumor suppressor activities of NRC-1, fine-scale physical mapping was conducted with a novel real-time quantitative PCR based method developed in this study. As a result, NRC-1 was mapped within a 4.6-Mb region defined by two unique sequences within UniGene clusters Hs.41407 and Hs.371835 (78,545Kb–83,172Kb in the NCBI build 31 physical map). The involvement of a putative tumor suppressor gene Robo1/Dutt1 was excluded as a candidate for NRC-1. Furthermore, a transcript map containing eleven candidate genes was established for the 4.6-Mb region. Analyses of gene expression patterns with real-time quantitative RT-PCR assays showed that one of the eleven candidate genes in the interval (TSGc28) is down-regulated in 15 out of 20 tumor samples compared with matched normal samples. Three exons of this gene have been identified by RACE experiments, although additional exon(s) seem to exist. Further gene characterization and functional studies are required to confirm the gene as a true tumor suppressor gene. ^ To study the cellular functions of NRC-1, gene expression profiles of three tumor suppressive microcell hybrids, each containing a functional copy of NRC-1, were compared with those of the corresponding parental tumor cell lines using 16K oligonucleotide microarrays. Differentially expressed genes were identified. Analyses based on the Gene Ontology showed that introduction of NRC-1 into tumor cell lines activates genes in multiple cellular pathways, including cell cycle, signal transduction, cytokines and stress response. NRC-1 is likely to induce cell growth arrest indirectly through WEE1. ^

ARD1 stabilization of TSC2 suppresses tumorigenesis via the mTOR signaling pathway

Mammalian target of rapamycin (mTOR) plays an important role in regulating various cellular functions, and the tuberous sclerosis 1 (TSC1)/TSC2 complex serves as a major repressor of the mTOR pathway. Here we demonstrated that arrest-defective protein 1 (ARD1) physically interacts with, acetylates, and stabilizes TSC2, thereby reducing mTOR activity. The inhibition of mTOR by ARD1 suppresses cell proliferation and increases autophagy, which further impairs tumorigenicity. Correlation between the levels of ARD1 and TSC2 was found in multiple tumor types, suggesting the physiological importance of ARD1 in stabilizing TSC2. Moreover, evaluation of loss of heterozygosity (LOH) at Xq28 revealed allelic loss in 31% of tested breast cancer cell lines and tumor samples. Together, our findings suggest that ARD1 functions as a negative regulator of the mTOR pathway and that dysregulation of the ARD1/TSC2/mTOR axis may contribute to cancer development. To further explore the signaling pathway of ARD1, we provided evidence showing the phosphorylation of ARD1 by IKKβ, which mediated the destabilization of ARD1. Future work may be needed to study the biological effect of this post-translational modification. ^

Caracterización molecular de variedades de vid ( Vitis vinifera L.) de calidad enológica por marcadores microsatelitales

Siete variedades de vid empleadas en Argentina para la elaboración de vinos de alta gama fueron caracterizadas molecularmente a través del empleo de microsatélites. Seis de los 8 pares de cebadores usados amplificaron patrones de bandas reproducibles. El número de alelos detectados por locus varió entre 5 y 10, con un total de 42 alelos, registrándose desde 1 a 7 alelos únicos. El número de genotipos microsatélites encontrados para cada locus osciló entre 3 y 7, con un total de 28. Las variedades Tempranillo y Chenin mostraron 6 y 5 genotipos «únicos», respectivamente, con los 6 loci analizados. La heterocigosidad observada varió entre 42,9 y 100 %, mientras que la heterocigosidad esperada osciló entre 64,3 y 87,8 %. El locus más informativo fue VrZAG79 (con un contenido de información polimórfica de 84,9 % y un número de alelos efectivos de 8,17) mientras que el menos informativo fue VrZAG62. La probabilidad acumulada de obtener genotipos idénticos fue de 7,68 x 10-05 indicando que, mediante el uso combinado de estos 6 cebadores microsatélites se pueden discriminar todas las variedades ensayadas, con alto grado de certeza. El análisis de agrupamiento por el método UPGMA separó claramente las variedades francesas (Merlot, Pinot Noir y Cabernet Sauvignon) y Syrah de los cepajes Tempranillo y Bonarda. Esta técnica podría ofrecerse como servicio a viveristas y productores vitícolas interesados en garantizar la identidad genética de sus materiales comercializados.

Description of two mass species of Pacific copepods in terms of two polymorphic loci

Population genetics of two species of mass copepods Undinula darwini and Calanus australis, with different range types, is investigated. Both species exhibit considerable genetic diversity, especially C. australis (observed heterozygoticity = 0.36), which inhabits a variable biotope in the zone of the Peru current. Samples of both species exhibited highly significant genetic heterogeneity as well as heterozygote deficiency compared with the situation expected from the Hardy-Weinberg law. Contribution of distance isolation to genetic differentiation of populations is estimated. Gene drift is discussed as a source of heterogeneity in populations of planktic copepods. Possible aspects of population genetic research on marine plank-tic crustaceans are discussed.