Study of gene flow through a hybrid zone in the common shrew (Sorex araneus) using microsatellites

The common shrew (Sorer araneus) is subdivided into several chromosomal races. As hybrid zones between them have been characterized, this organism is of particular interest in studying the role of chromosomes in speciation. Six microsatellite loci were used to evaluate the level of gene how in the S. araneus hybrid zone between the Cordon and Valais races. Most of these loci were very polymorphic, the total number of alleles detected per locus ranging from 3 to 20. Using Mantel tests, we showed that the effect of rivers as barriers to gene flow is less important at this sampling scale. The effect of the chromosomal race is of particular importantance in diminishing gene flow.

Molecular Modeling Approaches for Determining Gene Function: application to a Putative Poly-A Binding Protein from Leishmania amazonensis (LaPABP)

The great expansion in the number of genome sequencing projects has revealed the importance of computational methods to speed up the characterization of unknown genes. These studies have been improved by the use of three dimensional information from the predicted proteins generated by molecular modeling techniques. In this work, we disclose the structure-function relationship of a gene product from Leishmania amazonensis by applying molecular modeling and bioinformatics techniques. The analyzed sequence encodes a 159 aminoacids polypeptide (estimated 18 kDa) and was denoted LaPABP for its high homology with poly-A binding proteins from trypanosomatids. The domain structure, clustering analysis and a three dimensional model of LaPABP, basically obtained by homology modeling on the structure of the human poly-A binding protein, are described. Based on the analysis of the electrostatic potential mapped on the model's surface and conservation of intramolecular contacts responsible for folding stabilization we hypothesize that this protein may have less avidity to RNA than it's L. major counterpart but still account for a significant functional activity in the parasite. The model obtained will help in the design of mutagenesis experiments aimed to elucidate the mechanism of gene expression in trypanosomatids and serve as a starting point for its exploration as a potential source of targets for a rational chemotherapy.

Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study.

The identification of genetically homogeneous groups of individuals is a long standing issue in population genetics. A recent Bayesian algorithm implemented in the software STRUCTURE allows the identification of such groups. However, the ability of this algorithm to detect the true number of clusters (K) in a sample of individuals when patterns of dispersal among populations are not homogeneous has not been tested. The goal of this study is to carry out such tests, using various dispersal scenarios from data generated with an individual-based model. We found that in most cases the estimated 'log probability of data' does not provide a correct estimation of the number of clusters, K. However, using an ad hoc statistic DeltaK based on the rate of change in the log probability of data between successive K values, we found that STRUCTURE accurately detects the uppermost hierarchical level of structure for the scenarios we tested. As might be expected, the results are sensitive to the type of genetic marker used (AFLP vs. microsatellite), the number of loci scored, the number of populations sampled, and the number of individuals typed in each sample.

Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis.

Wounding in multicellular eukaryotes results in marked changes in gene expression that contribute to tissue defense and repair. Using a cDNA microarray technique, we analyzed the timing, dynamics, and regulation of the expression of 150 genes in mechanically wounded leaves of Arabidopsis. Temporal accumulation of a group of transcripts was correlated with the appearance of oxylipin signals of the jasmonate family. Analysis of the coronatine-insensitive coi1-1 Arabidopsis mutant that is also insensitive to jasmonate allowed us to identify a large number of COI1-dependent and COI1-independent wound-inducible genes. Water stress was found to contribute to the regulation of an unexpectedly large fraction of these genes. Comparing the results of mechanical wounding with damage by feeding larvae of the cabbage butterfly (Pieris rapae) resulted in very different transcript profiles. One gene was specifically induced by insect feeding but not by wounding; moreover, there was a relative lack of water stress-induced gene expression during insect feeding. These results help reveal a feeding strategy of P. rapae that may minimize the activation of a subset of water stress-inducible, defense-related genes.

Post-genomic update on a classical candidate gene for coronary artery disease: ESR1.

BACKGROUND: After age, sex is the most important risk factor for coronary artery disease (CAD). The mechanism through which women are protected from CAD is still largely unknown, but the observed sex difference suggests the involvement of the reproductive steroid hormone signaling system. Genetic association studies of the gene-encoding Estrogen Receptor α (ESR1) have shown conflicting results, although only a limited range of variation in the gene has been investigated. METHODS AND RESULTS: We exploited information made available by advanced new methods and resources in complex disease genetics to revisit the question of ESR1's role in risk of CAD. We performed a meta-analysis of 14 genome-wide association studies (CARDIoGRAM discovery analysis, N=≈87,000) to search for population-wide and sex-specific associations between CAD risk and common genetic variants throughout the coding, noncoding, and flanking regions of ESR1. In addition to samples from the MIGen (N=≈6000), WTCCC (N=≈7400), and Framingham (N=≈3700) studies, we extended this search to a larger number of common and uncommon variants by imputation into a panel of haplotypes constructed using data from the 1000 Genomes Project. Despite the widespread expression of ERα in vascular tissues, we found no evidence for involvement of common or low-frequency genetic variation throughout the ESR1 gene in modifying risk of CAD, either in the general population or as a function of sex. CONCLUSIONS: We suggest that future research on the genetic basis of sex-related differences in CAD risk should initially prioritize other genes in the reproductive steroid hormone biosynthesis system.

Isolation of a fragment homologous to the rp49 constitutive gene of Drosophila in the Neotropical malaria vector Anopheles aquasalis (Diptera: Culicidae)

The constitutive ribosomal gene rp49 is frequently used as an endogenous control in Drosophila gene expression experiments. Using the degenerate primer PCR technique we have cloned a fragment homologous to this gene in Anopheles aquasalis Curry, a Neotropical vector of malaria. In addition, based on this first sequence, a new primer was designed, which allowed the isolation of fragments of rp49 in two other species, Aedes aegypti (Linnaeus) and Culex quinquefasciatus Say, suggesting that it could be used to clone fragments of this gene in a number of other mosquito species. Primers were also designed to specifically amplify rp49 cDNA fragments in An. aquasalis and Ae. aegypti, showing that rp49 could be used as a good constitutive control in gene expression studies of these and other vectorially important mosquito species.

The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.

Development of a real time polymerase chain reaction for quantitation of Schistosoma mansoni DNA

This report describes the development of a SYBR Green I based real time polymerase chain reaction (PCR) protocol for detection on the ABI Prism 7000 instrument. Primers targeting the gene encoding the SSU rRNA were designed to amplify with high specificity DNA from Schistosoma mansoni, in a real time quantitative PCR system. The limit of detection of parasite DNA for the system was 10 fg of purified genomic DNA, that means less than the equivalent to one parasite cell (genome ~580 fg DNA). The efficiency was 0.99 and the correlation coefficient (R²) was 0.97. When different copy numbers of the target amplicon were used as standards, the assay could detect at least 10 copies of the specific target. The primers used were designed to amplify a 106 bp DNA fragment (Tm 83ºC). The assay was highly specific for S. mansoni, and did not recognize DNA from closely related non-schistosome trematodes. The real time PCR allowed for accurate quantification of S. mansoni DNA and no time-consuming post-PCR detection of amplification products by gel electrophoresis was required. The assay is potentially able to quantify S. mansoni DNA (and indirectly parasite burden) in a number of samples, such as snail tissue, serum and feces from patients, and cercaria infested water. Thus, these PCR protocols have potential to be used as tools for monitoring of schistosome transmission and quantitative diagnosis of human infection.

Gene finding in the chicken genome.

BACKGROUND: Despite the continuous production of genome sequence for a number of organisms, reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularly true for genomes for which there is not a large collection of known gene sequences, such as the recently published chicken genome. We used the chicken sequence to test comparative and homology-based gene-finding methods followed by experimental validation as an effective genome annotation method. RESULTS: We performed experimental evaluation by RT-PCR of three different computational gene finders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram was computed and each component of it was evaluated. The results showed that de novo comparative methods can identify up to about 700 chicken genes with no previous evidence of expression, and can correctly extend about 40% of homology-based predictions at the 5' end. CONCLUSIONS: De novo comparative gene prediction followed by experimental verification is effective at enhancing the annotation of the newly sequenced genomes provided by standard homology-based methods.

Single-stranded oligonucleotide-mediated in vivo gene repair in the rd1 retina.

PURPOSE: The aim of this study was to test whether oligonucleotide-targeted gene repair can correct the point mutation in genomic DNA of PDE6b(rd1) (rd1) mouse retinas in vivo. METHODS: Oligonucleotides (ODNs) of 25 nucleotide length and complementary to genomic sequence subsuming the rd1 point mutation in the gene encoding the beta-subunit of rod photoreceptor cGMP-phosphodiesterase (beta-PDE), were synthesized with a wild type nucleotide base at the rd1 point mutation position. Control ODNs contained the same nucleotide bases as the wild type ODNs but with varying degrees of sequence mismatch. We previously developed a repeatable and relatively non-invasive technique to enhance ODN delivery to photoreceptor nuclei using transpalpebral iontophoresis prior to intravitreal ODN injection. Three such treatments were performed on C3H/henJ (rd1) mouse pups before postnatal day (PN) 9. Treatment outcomes were evaluated at PN28 or PN33, when retinal degeneration was nearly complete in the untreated rd1 mice. The effect of treatment on photoreceptor survival was evaluated by counting the number of nuclei of photoreceptor cells and by assessing rhodopsin immunohistochemistry on flat-mount retinas and sections. Gene repair in the retina was quantified by allele-specific real time PCR and by detection of beta-PDE-immunoreactive photoreceptors. Confirmatory experiments were conducted using independent rd1 colonies in separate laboratories. These experiments had an additional negative control ODN that contained the rd1 mutant nucleotide base at the rd1 point mutation site such that the sole difference between treatment with wild type and control ODN was the single base at the rd1 point mutation site. RESULTS: Iontophoresis enhanced the penetration of intravitreally injected ODNs in all retinal layers. Using this delivery technique, significant survival of photoreceptors was observed in retinas from eyes treated with wild type ODNs but not control ODNs as demonstrated by cell counting and rhodopsin immunoreactivity at PN28. Beta-PDE immunoreactivity was present in retinas from eyes treated with wild type ODN but not from those treated with control ODNs. Gene correction demonstrated by allele-specific real time PCR and by counts of beta-PDE-immunoreactive cells was estimated at 0.2%. Independent confirmatory experiments showed that retinas from eyes treated with wild type ODN contained many more rhodopsin immunoreactive cells compared to retinas treated with control (rd1 sequence) ODN, even when harvested at PN33. CONCLUSIONS: Short ODNs can be delivered with repeatable efficiency to mouse photoreceptor cells in vivo using a combination of intravitreal injection and iontophoresis. Delivery of therapeutic ODNs to rd1 mouse eyes resulted in genomic DNA conversion from mutant to wild type sequence, low but observable beta-PDE immunoreactivity, and preservation of rhodopsin immunopositive cells in the outer nuclear layer, suggesting that ODN-directed gene repair occurred and preserved rod photoreceptor cells. Effects were not seen in eyes treated with buffer or with ODNs having the rd1 mutant sequence, a definitive control for this therapeutic approach. Importantly, critical experiments were confirmed in two laboratories by several different researchers using independent mouse colonies and ODN preparations from separate sources. These findings suggest that targeted gene repair can be achieved in the retina following enhanced ODN delivery.

Mice Transgenic For Human Dominant Mutations Of The GUCY2D Gene

Purpose: Animal models are essential to study pathological mechanisms and to test new therapeutic strategies. Many mouse models mimic human rod loss but only a limited number simulate cone dystrophies. The importance of cone function for human vision highlights the need to engineer a model for cone degeneration. An approach of lentiviral-directed transgenesis was tested in mice to express a dominant mutant gene described in a human cone dystrophy.Methods: Lentiviral vectors (LV) encoding either hrGFPII or the human double mutant GUCY2DE837D/R838S cDNA under the control of a region of the pig arrestin-3 promoter (Arr3) were produced and used for lentiviral-derived transgenesis. PCR-genotyping determined the transgenic mouse ratio. The expression of GFP was then analyzed both in vivo and by immunohistochemistry in Arr3-GFPII mice. Functional analysis was performed by ERG at 5, 9, 16 and 24 weeks for Arr3-GUCY2DE837D/R838S mice. Mice were sacrificed at 10 months of age for both histological analysis and RNA extraction.Results: While all the newborns from the transgenesis using the LV-Arr3-GFPII were transgenic, one third of the newborns from the LV-Arr3-GUCY2DE837D/R838S transgenesis were positive. Expression of GFPII was demonstrated by in vivo imaging, while expression of the mutant GUCY2D transcript was detetected using RT-PCR. No severe alteration of the functional response was observed up to 24 weeks of age in the transgenic mice. No obvious modification of the retinal morphology was identified either.Conclusions: Lentiviral-directed transgenesis is a rapid and straightforward method to engineer transgenic mice. Protein expression can be specifically targeted to the retina and thus could help to study the effect of expression of dominant mutant proteins. In our case, Arr3-GUCY2DE837D/R838S mice have a less severe phenotype than that described for human patients. Further analyses are required to understand this difference but several modifications of the expression cassette might also help to increase the expression of the mutant protein and reinforce the phenotype. Interestingly, the same construct is less effective in mouse versus pig retina (see Arsenijevic et al. ARVO 2011 abstract).

Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications.

BACKGROUND: The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes. RESULTS: We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons. CONCLUSION: There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular - but possibly clusters of genes more generally - might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

Genetic clusters and sex-biased gene flow in a unicolonial Formica ant.

BACKGROUND: Animal societies are diverse, ranging from small family-based groups to extraordinarily large social networks in which many unrelated individuals interact. At the extreme of this continuum, some ant species form unicolonial populations in which workers and queens can move among multiple interconnected nests without eliciting aggression. Although unicoloniality has been mostly studied in invasive ants, it also occurs in some native non-invasive species. Unicoloniality is commonly associated with very high queen number, which may result in levels of relatedness among nestmates being so low as to raise the question of the maintenance of altruism by kin selection in such systems. However, the actual relatedness among cooperating individuals critically depends on effective dispersal and the ensuing pattern of genetic structuring. In order to better understand the evolution of unicoloniality in native non-invasive ants, we investigated the fine-scale population genetic structure and gene flow in three unicolonial populations of the wood ant F. paralugubris. RESULTS: The analysis of geo-referenced microsatellite genotypes and mitochondrial haplotypes revealed the presence of cryptic clusters of genetically-differentiated nests in the three populations of F. paralugubris. Because of this spatial genetic heterogeneity, members of the same clusters were moderately but significantly related. The comparison of nuclear (microsatellite) and mitochondrial differentiation indicated that effective gene flow was male-biased in all populations. CONCLUSION: The three unicolonial populations exhibited male-biased and mostly local gene flow. The high number of queens per nest, exchanges among neighbouring nests and restricted long-distance gene flow resulted in large clusters of genetically similar nests. The positive relatedness among clustermates suggests that kin selection may still contribute to the maintenance of altruism in unicolonial populations if competition occurs among clusters.

Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.