Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55-85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or -4.0 to -1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD-associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies. © 2011 Duncan et al.

Genetic determinants of bone density and fracture risk: State of the art and future directions

Context: Osteoporosis is a common, highly heritable condition that causes substantial morbidity and mortality, the etiopathogenesis of which is poorly understood. Genetic studies are making increasingly rapid progress in identifying the genes involved. Evidence Acquisition and Synthesis: In this review, we will summarize the current understanding of the genetics of osteoporosis based on publications from PubMed from the year 1987 onward. Conclusions: Most genes involved in osteoporosis identified to date encode components of known pathways involved in bone synthesis or resorption, but as the field progresses, new pathways are being identified. Only a small proportion of the total genetic variation involved in osteoporosis has been identified, and new approaches will be required to identify most of the remaining genes.

Whole-genome sequencing reveals untapped genetic potential in Africa’s indigenous cereal crop sorghum

Sorghum is a food and feed cereal crop adapted to heat and drought and a staple for 500 million of the world’s poorest people. Its small diploid genome and phenotypic diversity make it an ideal C4 grass model as a complement to C3 rice. Here we present high coverage (16-45 × ) resequenced genomes of 44 sorghum lines representing the primary gene pool and spanning dimensions of geographic origin, end-use and taxonomic group. We also report the first resequenced genome of S. propinquum, identifying 8 M high-quality SNPs, 1.9 M indels and specific gene loss and gain events in S. bicolor. We observe strong racial structure and a complex domestication history involving at least two distinct domestication events. These assembled genomes enable the leveraging of existing cereal functional genomics data against the novel diversity available in sorghum, providing an unmatched resource for the genetic improvement of sorghum and other grass species.

Phycobilisomes linker family in cyanobacterial genomes: divergence and evolution

Cyanobacteria are the oldest life form making important contributions to global CO2 fixation on the Earth. Phycobilisomes (PBSs) are the major light harvesting systems of most cyanobacteria species. Recent availability of the whole genome database of cyanobacteria provides us a global and further view on the complex structural PBSs. A PBSs linker family is crucial in structure and function of major light-harvesting PBSs complexes. Linker polypeptides are considered to have the same ancestor with other phycobiliproteins (PBPs), and might have been diverged and evolved under particularly selective forces together. In this paper, a total of 192 putative linkers including 167 putative PBSs-associated linker genes and 25 Ferredoxin-NADP oxidoreductase (FNR) genes were detected through whole genome analysis of all 25 cyanobacterial genomes (20 finished and 5 in draft state). We compared the PBSs linker family of cyanobacteria in terms of gene structure, chromosome location, conservation domain, and polymorphic variants, and discussed the features and functions of the PBSs linker family. Most of PBSs-associated linkers in PBSs linker family are assembled into gene clusters with PBPs. A phylogenetic analysis based on protein data demonstrates a possibility of six classes of the linker family in cyanobacteria. Emergence, divergence, and disappearance of PBSs linkers among cyanobacterial species were due to speciation, gene duplication, gene transfer, or gene loss, and acclimation to various environmental selective pressures especially light.

RNAi Effector Diversity in Nematodes

While RNA interference (RNAi) has been deployed to facilitate gene function studies in diverse helminths, parasitic nematodes appear variably susceptible. To test if this is due to inter-species differences in RNAi effector complements, we performed a primary sequence similarity survey for orthologs of 77 Caenorhabditis elegans RNAi pathway proteins in 13 nematode species for which genomic or transcriptomic datasets were available, with all outputs subjected to domain-structure verification. Our dataset spanned transcriptomes of Ancylostoma caninum and Oesophagostomum dentatum, and genomes of Trichinella spiralis, Ascaris suum, Brugia malayi, Haemonchus contortus, Meloidogyne hapla, Meloidogyne incognita and Pristionchus pacificus, as well as the Caenorhabditis species C. brenneri, C. briggsae, C. japonica and C. remanei, and revealed that: (i) Most of the C. elegans proteins responsible for uptake and spread of exogenously applied double stranded (ds)RNA are absent from parasitic species, including RNAi-competent plant-nematodes; (ii) The Argonautes (AGOs) responsible for gene expression regulation in C. elegans are broadly conserved, unlike those recruited during the induction of RNAi by exogenous dsRNA; (iii) Secondary Argonautes (SAGOs) are poorly conserved, and the nuclear AGO NRDE-3 was not identified in any parasite; (iv) All five Caenorhabditis spp. possess an expanded RNAi effector repertoire relative to the parasitic nematodes, consistent with the propensity for gene loss in nematode parasites; (v) In spite of the quantitative differences in RNAi effector complements across nematode species, all displayed qualitatively similar coverage of functional protein groups. In summary, we could not identify RNAi effector deficiencies that associate with reduced susceptibility in parasitic nematodes. Indeed, similarities in the RNAi effector complements of RNAi refractory and competent nematode parasites support the broad applicability of this research genetic tool in nematodes.

The evolutionary history of the CD209 (DC-SIGN) family in humans and non-human primates

The CD209 gene family that encodes C-type lectins in primates includes CD209 (DC-SIGN), CD209L (L-SIGN) and CD209L2. Understanding the evolution of these genes can help understand the duplication events generating this family, the process leading to the repeated neck region and identify protein domains under selective pressure. We compiled sequences from 14 primates representing 40 million years of evolution and from three non-primate mammal species. Phylogenetic analyses used Bayesian inference, and nucleotide substitutional patterns were assessed by codon-based maximum likelihood. Analyses suggest that CD209 genes emerged from a first duplication event in the common ancestor of anthropoids, yielding CD209L2 and an ancestral CD209 gene, which, in turn, duplicated in the common Old World primate ancestor, giving rise to CD209L and CD209. K(A)/K(S) values averaged over the entire tree were 0.43 (CD209), 0.52 (CD209L) and 0.35 (CD209L2), consistent with overall signatures of purifying selection. We also assessed the Toll-like receptor (TLR) gene family, which shares with CD209 genes a common profile of evolutionary constraint. The general feature of purifying selection of CD209 genes, despite an apparent redundancy (gene absence and gene loss), may reflect the need to faithfully recognize a multiplicity of pathogen motifs, commensals and a number of self-antigens

Approches algorithmiques pour l’inférence d’histoires de duplication en tandem avec inversions et délétions pour des familles multigéniques

[Français] Une fraction importante des génomes eucaryotes est constituée de Gènes Répétés en Tandem (GRT). Un mécanisme fondamental dans l’évolution des GRT est la recombinaison inégale durant la méiose, entrainant la duplication locale (en tandem) de segments chromosomiques contenant un ou plusieurs gènes adjacents. Différents algorithmes ont été proposés pour inférer une histoire de duplication en tandem pour un cluster de GRT. Cependant, leur utilisation est limitée dans la pratique, car ils ne tiennent pas compte d’autres événements évolutifs pourtant fréquents, comme les inversions, les duplications inversées et les délétions. Cette thèse propose différentes approches algorithmiques permettant d’intégrer ces événements dans le modèle de duplication en tandem classique. Nos contributions sont les suivantes: • Intégrer les inversions dans un modèle de duplication en tandem simple (duplication d’un gène à la fois) et proposer un algorithme exact permettant de calculer le nombre minimal d’inversions s’étant produites dans l’évolution d’un cluster de GRT. • Généraliser ce modèle pour l’étude d’un ensemble de clusters orthologues dans plusieurs espèces. • Proposer un algorithme permettant d’inférer l’histoire évolutive d’un cluster de GRT en tenant compte des duplications en tandem, duplications inversées, inversions et délétions de segments chromosomiques contenant un ou plusieurs gènes adjacents.

Algorithmes pour la réconciliation d’un arbre de gènes avec un arbre d’espèces

Une réconciliation entre un arbre de gènes et un arbre d’espèces décrit une histoire d’évolution des gènes homologues en termes de duplications et pertes de gènes. Pour inférer une réconciliation pour un arbre de gènes et un arbre d’espèces, la parcimonie est généralement utilisée selon le nombre de duplications et/ou de pertes. Les modèles de réconciliation sont basés sur des critères probabilistes ou combinatoires. Le premier article définit un modèle combinatoire simple et général où les duplications et les pertes sont clairement identifiées et la réconciliation parcimonieuse n’est pas la seule considérée. Une architecture de toutes les réconciliations est définie et des algorithmes efficaces (soit de dénombrement, de génération aléatoire et d’exploration) sont développés pour étudier les propriétés combinatoires de l’espace de toutes les réconciliations ou seulement les plus parcimonieuses. Basée sur le processus classique nommé naissance-et-mort, un algorithme qui calcule la vraisemblance d’une réconciliation a récemment été proposé. Le deuxième article utilise cet algorithme avec les outils combinatoires décrits ci-haut pour calculer efficacement (soit approximativement ou exactement) les probabilités postérieures des réconciliations localisées dans le sous-espace considéré. Basé sur des taux réalistes (selon un modèle probabiliste) de duplication et de perte et sur des données réelles/simulées de familles de champignons, nos résultats suggèrent que la masse probabiliste de toute l’espace des réconciliations est principalement localisée autour des réconciliations parcimonieuses. Dans un contexte d’approximation de la probabilité d’une réconciliation, notre approche est une alternative intéressante face aux méthodes MCMC et peut être meilleure qu’une approche sophistiquée, efficace et exacte pour calculer la probabilité d’une réconciliation donnée. Le problème nommé Gene Tree Parsimony (GTP) est d’inférer un arbre d’espèces qui minimise le nombre de duplications et/ou de pertes pour un ensemble d’arbres de gènes. Basé sur une approche qui explore tout l’espace des arbres d’espèces pour les génomes considérés et un calcul efficace des coûts de réconciliation, le troisième article décrit un algorithme de Branch-and-Bound pour résoudre de façon exacte le problème GTP. Lorsque le nombre de taxa est trop grand, notre algorithme peut facilement considérer des relations prédéfinies entre ensembles de taxa. Nous avons testé notre algorithme sur des familles de gènes de 29 eucaryotes.

Caractérisation du rôle et des mécanismes d’action des gènes Hoxa dans l’hématopoïèse adulte

Chez les humains, un large pourcentage de leucémies myéloïdes et lymphoïdes exprime des gènes Homéobox (Hox) de façon aberrante, principalement ceux du groupe des gènes Hoxa. Cette dérégulation de l’expression des gènes Hox peut provenir directement des translocations impliquant des gènes Hox ou indirectement par d’autres protéines ayant un potentiel oncogénique. De plus, plusieurs études indiquent que les gènes Hox jouent un rôle essentiel dans l'initiation de diverses leucémies. Comprendre le fonctionnement des gènes Hox dans l'hématopoïèse normale est donc une condition préalable pour élucider leurs fonctions dans les leucémies, ce qui pourrait éventuellement conduire à l’élaboration de nouveaux traitements contre cette maladie. Plusieurs études ont tenté d’élucider les rôles exacts des gènes Hox dans l'hématopoïèse via l’utilisation de souris mutantes pour un seul gène Hox. Or, en raison du phénomène de redondance fonctionnelle chez cette famille de gènes, ces études ont été peu concluantes. Il a été précédemment démontré que dans une population de cellules enrichies en cellules souches hématopoïétiques (CSH), les gènes du cluster Hoxa sont plus exprimés que les gènes Hox des autres clusters. Aussi, il a été établi que les gènes du cluster Hoxb sont non essentiels à l’hématopoïèse définitive puisque les CSH mutantes pour les gènes Hoxb1-9 conservent leur potentiel de reconstitution à long terme. En nous basant sur ces données, nous avons émis l'hypothèse suivante : les gènes Hoxa sont essentiels pour l'hématopoïèse normale adulte. Pour tester notre hypothèse, nous avons choisi d’utiliser un modèle de souris comportant une délétion pour l’ensemble des gènes Hoxa. Dans le cadre de cette recherche, nous avons démontré que les CSH, les progéniteurs primitifs et les progéniteurs des cellules B sont particulièrement sensibles au niveau d'expression des gènes Hoxa. Plus particulièrement, une baisse de la survie et une différenciation prématurée semblent être à l’origine de la perte des CSH Hoxa-/- dans la moelle osseuse. L’analyse du profil transcriptionnel des CSH par séquençage de l'ARN a révélé que les gènes Hoxa sont capables de réguler un vaste réseau de gènes impliqués dans divers processus biologiques. En effet, les gènes Hoxa régulent l’expression de plusieurs gènes codant pour des récepteurs de cytokine. De plus, les gènes Hoxa influencent l’expression de gènes jouant une fonction dans l’architecture de la niche hématopoïétique. L’expression de plusieurs molécules d’adhésion est aussi modulée par les gènes Hoxa, ce qui peut affecter la relation des CSH avec la niche hématopoïétique. L’ensemble de ces résultats démontre que les gènes Hoxa sont d'importants régulateurs de l'hématopoïèse adulte puisqu’ils sont nécessaires au maintien des CSH et des progéniteurs grâce à leurs effets sur plusieurs processus biologiques comme l'apoptose, le cycle cellulaire et les interactions avec la niche.

Complex history of a chromosomal paralogy region: insights from amphioxus aromatic amino acid hydroxylase genes and insulin-related genes

Aromatic amino acid hydroxylase (AAAH) genes and insulin-like genes form part of an extensive paralogy region shared by human chromosomes 11 and 12, thought to have arisen by tetraploidy in early vertebrate evolution. Cloning of a complementary DNA (cDNA) for an amphioxus (Branchiostoma floridae) hydroxylase gene (AmphiPAH) allowed us to investigate the ancestry of the human chromosome 11/12 paralogy region. Molecular phylogenetic evidence reveals that AmphiPAH is orthologous to vertebrate phenylalanine (PAH) genes; the implication is that all three vertebrate AAAH genes arose early in metazoan evolution, predating vertebrates. In contrast, our phylogenetic analysis of amphioxus and vertebrate insulin-related gene sequences is consistent with duplication of these genes during early chordate ancestry. The conclusion is that two tightly linked gene families on human chromosomes 11 and 12 were not duplicated coincidentally. We rationalize this paradox by invoking gene loss in the AAAH gene family and conclude that paralogous genes shared by paralogous chromosomes need not have identical evolutionary histories.

Divergent evolutionary pattern of starch biosynthetic pathway genes in grasses and dicots

Starch is the most widespread and abundant storage carbohydrate in crops and its production is critical to both crop yield and quality. As regards the starch content in the seeds of crop plants, there are distinct difference between grasses (Poaceae) and dicots. However, few studies have described the evolutionary pattern of genes in the starch biosynthetic pathway in these two groups of plants. In this study, therefore, an attempt was made to compare the evolutionary rate, gene duplication and selective pattern of the key genes involved in this pathway between the two groups, using five grasses and five dicots as materials. The results showed (i) distinct differences in patterns of gene duplication and loss between grasses and dicots; duplication in grasses mainly occurred prior to the divergence of grasses, whereas duplication mostly occurred in individual species within the dicots; there is less gene loss in grasses than in dicots; (ii) a considerably higher evolutionary rate in grasses than in dicots in most gene families analyzed; (iii) evidence of a different selective pattern between grasses and dicots; positive selection may have occurred asymmetrically in grasses in some gene families, e.g. AGPase small subunit. Therefore, we deduced that gene duplication contributes to, and a higher evolutionary rate is associated with, the higher starch content in grasses. In addition, two novel aspects of the evolution of the starch biosynthetic pathway were observed.

Evaluation of monocot and eudicot divergence using the sugarcane transcriptome(1[w])

Over 40,000 sugarcane (Saccharum officinarum) consensus sequences assembled from 237,954 expressed sequence tags were compared with the protein and DNA sequences from other angiosperms, including the genomes of Arabidopsis and rice (Oryza sativa). Approximately two-thirds of the sugarcane transcriptome have similar sequences in Arabidopsis. These sequences may represent a core set of proteins or protein domains that are conserved among monocots and eudicots and probably encode for essential angiosperm. functions. The remaining sequences represent putative monocot-specific genetic material, one-half of which were found only in sugarcane. These monocot-specific cDNAs represent either novelties or, in many cases, fast-evolving sequences that diverged substantially from their eudicot homologs. The wide comparative genome analysis presented here provides information on the evolutionary changes that underlie the divergence of monocots and eudicots. Our comparative analysis also led to the identification of several not yet annotated putative genes and possible gene loss events in Arabidopsis.

Quantitative real-time PCR identifies a critical region of deletion on 22q13 related to prognosis in oral cancer

Quantitative real time PCR was performed on genomic DNA from 40 primary oral carcinomas and the normal adjacent tissues. The target genes ECGFB, DIA1, BIK, and PDGFB and the microsatellite markers D22S274 and D22S277, mapped on 22q13, were selected according to our previous loss of heterozygosity findings in head and neck tumors. Quantitative PCR relies on the comparison of the amount of product generated from a target gene and that generated from a disomic reference gene (GAPDH-housekeeping gene). Reactions have been performed with normal control in triplicates, using the 7700 Sequence Detection System (PE Applied Biosystems). Losses in the sequences D22S274 (22q13.31) and in the DIA1 (22q13.2-13.31) gene were detected in 10 out of 40 cases (25%) each. Statistically significant correlations were observed for patients with relative copy number loss of the marker D22S274 and stages T3-T4 of disease (P=0.025), family history of cancer (P = 0.001), and death (P = 0.021). Relative copy number loss involving the DIA1 gene was correlated to family history of cancer (P<0.001), death (P=0.002), and consumption of alcohol (P=0.026). Log-rank test revealed a significant decrease in survival (P=0.0018) for patients with DIA1 gene loss. Relative copy number losses detected in these sequences may be related to disease progression and a worse prognosis in patients with oral cancer.

Diversidade genética de isolados ambientais de Vibrio cholerae da Amazônia brasileira

Vibrio cholerae, agente etiológico da cólera, é uma bactéria nativa de ambientes aquáticos de regiões temperadas e tropicais em todo o mundo. A cólera é endemica e epidemica em países da África, Ásia e Americas Central e do Sul. Neste trabalho o objetivo foi estudar a diversidade genética de isolados desta espécie, de ambientes aquáticos da Amazônia brasileira. Um total de 148 isolados de V.cholerae não-O1 e não-O139 (NAGs) e O1 ambientais da Amazônia, obtidos entre 1977 e 2007, foram caracterizados e comparados a linhagens clínicas de V.cholerae O1 da sexta e sétima pandemias. Utilizou-se os perfis de macrorestrição definidos em eletroforese em gel de agarose em campo pulsado (PFGE), para determinar a relação clonal entre V.cholerae non-O1 e O1 ambientais e clínicos. A presença de genes de virulência (hlyA/hem, hlyB, hlyC, rtxA, rtxC, tcp, ctx, zot, ace, stn/sto) e integrons de classe 1, 2 e 3 (intI 1, 2 e 3), foi analisada utilizando-se a reação em cadeia da polimerase. A análise dos perfis de macrorestrição revelou que os NAGs apresentaram uma grande diversidade genética comparada aos V.cholerae O1. Isolados de NAGs e O1 segregaram em distintos grupos e a maioria dos O1 ambientais apresentou relação clonal com isolados clínicos da sétima pandemia de cólera. A distribuição dos genes de virulência entre os NAGs é diferente a dos O1, os quais, em geral, foram positivos para todos os genes de virulência estudados exceto stn/sto e integrons de classe 1, 2 e 3. Alguns V.cholerae O1 ambientais pertencentes a linhagem da sétima pandemia, apresentaram uma extensiva perda de genes. Diferentes NAGs foram stn/sto+ e intI 1+. Dois alelos do gene aadA foram encontrados: aadA2 e aadA7. De modo interessante os V.cholerae O1 ambientais pertencentes à linhagem pandêmica, só foram isolados durante o período da última epidemia de cólera na região Amazônica brasileira (1991-1996).

Different patterns of evolution for duplicated DNA repair genes in bacteria of the Xanthomonadales group

Abstract Background DNA repair genes encode proteins that protect organisms against genetic damage generated by environmental agents and by-products of cell metabolism. The importance of these genes in life maintenance is supported by their high conservation, and the presence of duplications of such genes may be easily traced, especially in prokaryotic genomes. Results The genome sequences of two Xanthomonas species were used as the basis for phylogenetic analyses of genes related to DNA repair that were found duplicated. Although 16S rRNA phylogenetic analyses confirm their classification at the basis of the gamma proteobacteria subdivision, differences were found in the origin of the various genes investigated. Except for lexA, detected as a recent duplication, most of the genes in more than one copy are represented by two highly divergent orthologs. Basically, one of such duplications is frequently positioned close to other gamma proteobacteria, but the second is often positioned close to unrelated bacteria. These orthologs may have occurred from old duplication events, followed by extensive gene loss, or were originated from lateral gene transfer (LGT), as is the case of the uvrD homolog. Conclusions Duplications of DNA repair related genes may result in redundancy and also improve the organisms' responses to environmental challenges. Most of such duplications, in Xanthomonas, seem to have arisen from old events and possibly enlarge both functional and evolutionary genome potentiality.