Enlarging a training set for genomic selction by imputation of un-genotyped animals in populations of varying genetic architecture

Background: The most common application of imputation is to infer genotypes of a high-density panel of markers on animals that are genotyped for a low-density panel. However, the increase in accuracy of genomic predictions resulting from an increase in the number of markers tends to reach a plateau beyond a certain density. Another application of imputation is to increase the size of the training set with un-genotyped animals. This strategy can be particularly successful when a set of closely related individuals are genotyped. ----- Methods: Imputation on completely un-genotyped dams was performed using known genotypes from the sire of each dam, one offspring and the offspring’s sire. Two methods were applied based on either allele or haplotype frequencies to infer genotypes at ambiguous loci. Results of these methods and of two available software packages were compared. Quality of imputation under different population structures was assessed. The impact of using imputed dams to enlarge training sets on the accuracy of genomic predictions was evaluated for different populations, heritabilities and sizes of training sets. ----- Results: Imputation accuracy ranged from 0.52 to 0.93 depending on the population structure and the method used. The method that used allele frequencies performed better than the method based on haplotype frequencies. Accuracy of imputation was higher for populations with higher levels of linkage disequilibrium and with larger proportions of markers with more extreme allele frequencies. Inclusion of imputed dams in the training set increased the accuracy of genomic predictions. Gains in accuracy ranged from close to zero to 37.14%, depending on the simulated scenario. Generally, the larger the accuracy already obtained with the genotyped training set, the lower the increase in accuracy achieved by adding imputed dams. ----- Conclusions: Whenever a reference population resembling the family configuration considered here is available, imputation can be used to achieve an extra increase in accuracy of genomic predictions by enlarging the training set with completely un-genotyped dams. This strategy was shown to be particularly useful for populations with lower levels of linkage disequilibrium, for genomic selection on traits with low heritability, and for species or breeds for which the size of the reference population is limited.

Polimorfismos del gen Butirilcolinesterasa responsables de reacciones adversas en pacientes consumidores de “cocaína”

La butirilcolinesterasa humana (BChE; EC 3.1.1.8) es una enzima polimórfica sintetizada en el hígado y en el tejido adiposo, ampliamente distribuida en el organismo y encargada de hidrolizar algunos ésteres de colina como la procaína, ésteres alifáticos como el ácido acetilsalicílico, fármacos como la metilprednisolona, el mivacurium y la succinilcolina y drogas de uso y/o abuso como la heroína y la cocaína. Es codificada por el gen BCHE (OMIM 177400), habiéndose identificado más de 100 variantes, algunas no estudiadas plenamente, además de la forma más frecuente, llamada usual o silvestre. Diferentes polimorfismos del gen BCHE se han relacionado con la síntesis de enzimas con niveles variados de actividad catalítica. Las bases moleculares de algunas de esas variantes genéticas han sido reportadas, entre las que se encuentra las variantes Atípica (A), fluoruro-resistente del tipo 1 y 2 (F-1 y F-2), silente (S), Kalow (K), James (J) y Hammersmith (H). En este estudio, en un grupo de pacientes se aplicó el instrumento validado Lifetime Severity Index for Cocaine Use Disorder (LSI-C) para evaluar la gravedad del consumo de “cocaína” a lo largo de la vida. Además, se determinaron Polimorfismos de Nucleótido Simple (SNPs) en el gen BCHE conocidos como responsables de reacciones adversas en pacientes consumidores de “cocaína” mediante secuenciación del gen y se predijo el efecto delos SNPs sobre la función y la estructura de la proteína, mediante el uso de herramientas bio-informáticas. El instrumento LSI-C ofreció resultados en cuatro dimensiones: consumo a lo largo de la vida, consumo reciente, dependencia psicológica e intento de abandono del consumo. Los estudios de análisis molecular permitieron observar dos SNPs codificantes (cSNPs) no sinónimos en el 27.3% de la muestra, c.293A>G (p.Asp98Gly) y c.1699G>A (p.Ala567Thr), localizados en los exones 2 y 4, que corresponden, desde el punto de vista funcional, a la variante Atípica (A) [dbSNP: rs1799807] y a la variante Kalow (K) [dbSNP: rs1803274] de la enzima BChE, respectivamente. Los estudios de predicción In silico establecieron para el SNP p.Asp98Gly un carácter patogénico, mientras que para el SNP p.Ala567Thr, mostraron un comportamiento neutro. El análisis de los resultados permite proponer la existencia de una relación entre polimorfismos o variantes genéticas responsables de una baja actividad catalítica y/o baja concentración plasmática de la enzima BChE y algunas de las reacciones adversas ocurridas en pacientes consumidores de cocaína.

Búsqueda de selección en el polimorfismo 677C>T (c.665C>T) del gen de la metilentetrahidrofolato reductasa (MTHFR) en una población Colombiana

Se realizó un estudio genético – poblacional en dos grupos etarios de población colombiana con la finalidad de evaluar las diferencias genéticas relacionadas con el polimorfismo MTHFR 677CT en busca de eventos genéticos que soporten la persistencia de este polimorfismo en la especie humana debido que este ha sido asociado con múltiples enfermedades. De esta manera se genotipificaron los individuos, se analizaron los genotipos, frecuencias alélicas y se realizaron diferentes pruebas genéticas-poblacionales. Contrario a lo observado en poblaciones Colombianas revisadas se identificó la ausencia del Equilibrio Hardy-Weinberg en el grupo de los niños y estructuras poblacionales entre los adultos lo que sugiere diferentes historias demográficas y culturales entre estos dos grupos poblacionales al tiempo, lo que soporta la hipótesis de un evento de selección sobre el polimorfismo en nuestra población. De igual manera nuestros datos fueron analizados junto con estudios previos a nivel nacional y mundial lo cual sustenta que el posible evento selectivo es debido a que el aporte de ácido fólico se ha incrementado durante las últimas dos décadas como consecuencia de las campañas de fortificación de las harinas y suplementación a las embarazadas con ácido fólico, por lo tanto aquí se propone un modelo de selección que se ajusta a los datos encontrados en este trabajo se establece una relación entre los patrones nutricionales de la especie humana a través de la historia que explica las diferencias en frecuencias de este polimorfismo a nivel espacial y temporal.  

Análisis del gen adra2a receptor alfa 2a adrenergico en pacientes con trastorno de hiperactividad y déficit de atención

El trastorno de hiperactividad y déficit de atención (THDA), es definido clínicamente como una alteración en el comportamiento, caracterizada por inatención, hiperactividad e impulsividad. Estos aspectos son clasificados en tres subtipos, que son: Inatento, hiperactivo impulsivo y mixto. Clínicamente se describe un espectro amplio que incluye desordenes académicos, trastornos de aprendizaje, déficit cognitivo, trastornos de conducta, personalidad antisocial, pobres relaciones interpersonales y aumento de la ansiedad, que pueden continuar hasta la adultez. A nivel global se ha estimado una prevalencia entre el 1% y el 22%, con amplias variaciones, dadas por la edad, procedencia y características sociales. En Colombia, se han realizado estudios en Bogotá y Antioquia, que han permitido establecer una prevalencia del 5% y 15%, respectivamente. La causa específica no ha sido totalmente esclarecida, sin embargo se ha calculado una heredabilidad cercana al 80% en algunas poblaciones, demostrando el papel fundamental de la genética en la etiología de la enfermedad. Los factores genéticos involucrados se relacionan con cambios neuroquímicos de los sistemas dopaminérgicos, serotoninérgicos y noradrenérgicos, particularmente en los sistemas frontales subcorticales, corteza cerebral prefrontal, en las regiones ventral, medial, dorsolateral y la porción anterior del cíngulo. Basados en los datos de estudios previos que sugieren una herencia poligénica multifactorial, se han realizado esfuerzos continuos en la búsqueda de genes candidatos, a través de diferentes estrategias. Particularmente los receptores Alfa 2 adrenérgicos, se encuentran en la corteza cerebral, cumpliendo funciones de asociación, memoria y es el sitio de acción de fármacos utilizados comúnmente en el tratamiento de este trastorno, siendo esta la principal evidencia de la asociación de este receptor con el desarrollo del THDA. Hasta la fecha se han descrito más de 80 polimorfismos en el gen (ADRA2A), algunos de los cuales se han asociado con la entidad. Sin embargo, los resultados son controversiales y varían según la metodología diagnóstica empleada y la población estudiada, antecedentes y comorbilidades. Este trabajo pretende establecer si las variaciones en la secuencia codificante del gen ADRA2A, podrían relacionarse con el fenotipo del Trastorno de Hiperactividad y el Déficit de Atención.

Bacterial evolution by genomic island transfer occurs via DNA transformation in planta

Our understanding of the evolution of microbial pathogens has been advanced by the discovery of "islands" of DNA that differ from core genomes and contain determinants of virulence [1, 2]. The acquisition of genomic islands (GIs) by horizontal gene transfer (HGT) is thought to have played a major role in microbial evolution. There are, however, few practical demonstrations of the acquisition of genes that control virulence, and, significantly, all have been achieved outside the animal or plant host. Loss of a GI from the bean pathogen Pseudomonas syringae pv. phaseolicola (Pph) is driven by exposure to the stress imposed by the plant's resistance response [3]. Here, we show that the complete episomal island, which carries pathogenicity genes including the effector avrPphB, transfers between strains of Pph by transformation in planta and inserts at a specific att site in the genome of the recipient. Our results show that the evolution of bacterial pathogens by HGT may be achieved via transformation, the simplest mechanism of DNA exchange. This process is activated by exposure to plant defenses, when the pathogen is in greatest need of acquiring new genetic traits to alleviate the antimicrobial stress imposed by plant innate immunity [4].

Genomic and genetic analyses of diversity and plant interactions of Pseudomonas fluorescens

Background: Pseudomonas fluorescens are common soil bacteria that can improve plant health through nutrient cycling, pathogen antagonism and induction of plant defenses. The genome sequences of strains SBW25 and Pf0-1 were determined and compared to each other and with P. fluorescens Pf-5. A functional genomic in vivo expression technology (IVET) screen provided insight into genes used by P. fluorescens in its natural environment and an improved understanding of the ecological significance of diversity within this species. Results: Comparisons of three P. fluorescens genomes (SBW25, Pf0-1, Pf-5) revealed considerable divergence: 61% of genes are shared, the majority located near the replication origin. Phylogenetic and average amino acid identity analyses showed a low overall relationship. A functional screen of SBW25 defined 125 plant-induced genes including a range of functions specific to the plant environment. Orthologues of 83 of these exist in Pf0-1 and Pf-5, with 73 shared by both strains. The P. fluorescens genomes carry numerous complex repetitive DNA sequences, some resembling Miniature Inverted-repeat Transposable Elements (MITEs). In SBW25, repeat density and distribution revealed 'repeat deserts' lacking repeats, covering approximately 40% of the genome. Conclusions: P. fluorescens genomes are highly diverse. Strain-specific regions around the replication terminus suggest genome compartmentalization. The genomic heterogeneity among the three strains is reminiscent of a species complex rather than a single species. That 42% of plant-inducible genes were not shared by all strains reinforces this conclusion and shows that ecological success requires specialized and core functions. The diversity also indicates the significant size of genetic information within the Pseudomonas pan genome.

Non-genomic signalling of the retinoic X receptor through inhibition of Gq signalling in human platelets

Retinoid X receptors (RXRs) are important transcriptional nuclear hormone receptors, acting as either homodimers or the binding partner for at least one fourth of all the known human nuclear receptors. Functional nongenomic effects of nuclear receptors are poorly understood; however, recently peroxisome proliferator-activated receptor (PPAR) gamma, PPARbeta, and the glucocorticoid receptor have all been found active in human platelets. Human platelets express RXRalpha and RXRbeta. RXR ligands inhibit platelet aggregation and TXA(2) release to ADP and the TXA(2) receptors, but only weakly to collagen. ADP and TXA(2) both signal via the G protein, Gq. RXR rapidly binds Gq but not Gi/z/o/t/gust in a ligand-dependent manner and inhibits Gq-induced Rac activation and intracellular calcium release. We propose that RXR ligands may have beneficial clinical actions through inhibition of platelet activation. Furthermore, our results demonstrate a novel nongenomic mode for nuclear receptor action and a functional cross-talk between G-protein and nuclear receptor signaling families.

Genomic analysis of recombinant Sabin clinical isolates

Recombination in Poliovirus vaccine strains is a very frequent phenomenon. In this report 23 polio/Sabin strains isolated from healthy vaccinees or from VAPP patients after OPV administration, were investigated in order to identify recombination sites from 2C to 3D regions of the poliovirus genome. RT-PCR, followed by Restriction Fragment Length Polymorphism (RFLP) screening analysis were applied in four distant genomic regions (5' UTR, VP1, 2C and 3C-3D) in order to detect any putative recombinant. The detected recombinants were sequenced from 2C to the end of the genome (3' UTR) and the exact recombination sites were determined with computational analysis. Five of the 23 isolated strains were recombinant in one genomic region, two of them in 2C, isolates EP16:S3/S2, EP23:S3/S1, two in 3D isolates EP6:S2/S1, EP12:S2/S1 and one in 3A isolate EP9:S2/Sl. Point mutations were found in strains EP3, EP6, EP9 and EP12. Recombination specific types and sites re-occurrence along with point mutations are discussed concerning the polioviruses evolution.

On the structural differences between markers and genomic AC microsatellites

AC microsatellites have proved particularly useful as genetic markers. For some purposes, such as in population biology, the inferences drawn depend on the quantitative values of their mutation rates. This, together with intrinsic biological interest, has led to widespread study of microsatellite mutational mechanisms. Now, however, inconsistencies are appearing in the results of marker-based versus non-marker-based studies of mutational mechanisms. The reasons for this have not been investigated, but one possibility, pursued here, is that the differences result from structural differences between markers and genomic microsatellites. Here we report a comparison between the CEPH AC marker microsatellites and the global population of AC microsatellites in the human genome. AC marker microsatellites are longer than the global average. Controlling for length, marker microsatellites contain on average fewer interruptions, and have longer segments, than their genomic counterparts. Related to this, marker microsatellites show a greater tendency to concentrate the majority of their repeats into one segment. These differences plausibly result from scientists selecting markers for their high polymorphism. In addition to the structural differences, there are differences in the base composition of flanking sequences, marker flanking regions being richer in C and G and poorer in A and T. Our results indicate that there are profound differences between marker and genomic microsatellites that almost certainly affect their mutation rates. There is a need for a unified model of mutational mechanisms that accounts for both marker-derived and genomic observations. A suggestion is made as to how this might be done.

The M1 matrix protein controls the filamentous phenotype of influenza A virus

We show that most isolates of influenza A induce filamentous changes in infected cells in contrast to A/WSN/33 and A/PR8/34 strains which have undergone extensive laboratory passage and are mouse-adapted. Using reverse genetics, we created recombinant viruses in the naturally filamentous genetic background of A/Victoria/3/75 and established that this property is regulated by the M1 protein sequence, but that the phenotype is complex and several residues are involved. The filamentous phenotype was lost when the amino acid at position 41 was switched from A to V, at the same time, this recombinant virus also became insensitive to the antibody 14C2. On the other hand, the filamentous phenotype could be fully transferred to a virus containing RNA segment 7 of the A/WSN/33 virus by a combination of three mutations in both the amino and carboxy regions of the M1 protein. This observation suggests that an interaction among these regions of M1 may occur during assembly. (C) 2004 Elsevier Inc. All rights reserved.

The chromosomal assessment of salt tolerant substituted tritipyrum using genomic fluorescent in situ hybridiization (FISH)

Wheat, although moderately tolerant to salt, can not be cultivated in many areas. However, in the triticeae tribe, some of the wild wheat relatives are highly tolerant, e.g. Thinopyrum bessarabicum, which grows on the sea shore. Eight primary hexaploid tritipyrum lines, amphiploids between Triticum durum and Thinopyrum bessarabicum have been produced which can set seed in at least 250 mM NaCl. These tritipyrums (2n=6x=42, AABBEbEb) due to reasons such as brittle rachis, continuous production of tillers, late maturity, tall stature and meiotic instability will not fulfill the requirements of a successful commercial salt tolerant crop. To overcome such problems the substituted tritipyrum, in which selected Eb chromosomes are replaced by D genome chromosomes of 6x wheat, was produced from 6x tritipyrum x 6x wheat hybrids (F1: 2n=6x=42, AABBDEb) followed by selfing and backcrossing with 6x tritipyrum. The fertile plants among the above progenies were screened by the genomic fluorescent in situ hybridization technique to identify their Eb and D chromosome constitution. This study showed that producing tritiprum with variable numbers of Eb and D genome chromosomes is feasible and that FISH is a useful technique for determining the number of Eb chromosomes present.

Confocal imaging of pseudomonas syringae pv. phaseolicola colony development in bean reveals reduced multiplication of strains containing the genomic island PPHGI-1.

Pseudomonas syringae pv. phaseolicola is the seed borne causative agent of halo blight in the common bean Phaseolus vulgaris. Pseudomonas syringae pv. phaseolicola race 4 strain 1302A contains the avirulence gene hopAR1 (located on a 106-kb genomic island, PPHGI-1, and earlier named avrPphB), which matches resistance gene R3 in P. vulgaris cultivar Tendergreen (TG) and causes a rapid hypersensitive reaction (HR). Here, we have fluorescently labeled selected Pseudomonas syringae pv. phaseolicola 1302A and 1448A strains (with and without PPHGI-1) to enable confocal imaging of in-planta colony formation within the apoplast of resistant (TG) and susceptible (Canadian Wonder [CW]) P. vulgaris leaves. Temporal quantification of fluorescent Pseudomonas syringae pv. phaseolicola colony development correlated with in-planta bacterial multiplication (measured as CFU/ml) and is, therefore, an effective means of monitoring Pseudomonas syringae pv. phaseolicola endophytic colonization and survival in P. vulgaris. We present advances in the application of confocal microscopy for in-planta visualization of Pseudomonas syringae pv. phaseolicola colony development in the leaf mesophyll to show how the HR defense response greatly affects colony morphology and bacterial survival. Unexpectedly, the presence of PPHGI-1 was found to cause a reduction of colony development in susceptible P. vulgaris CW leaf tissue. We discuss the evolutionary consequences that the acquisition and retention of PPHGI-1 brings to Pseudomonas syringae pv. phaseolicola in planta.