BMC Genetics - Latest Articles

Hypermethylation of the DLC1 CpG island does not alter gene expression in canine lymphoma

Jeffrey Bryan — 2009-11-13T00:00:00Z

Background: This study is a comparative epigenetic evaluation of the methylation status of the DLC1 tumor suppressor gene in naturally-occurring canine lymphoma. Canine non-Hodgkin's lymphoma (NHL) has been proposed to be a relevant preclinical model that occurs spontaneously and may share causative factors with human NHL due to a shared home environment. The canine DLC1 mRNA sequence was derived from normal tissue. Using lymphoid samples from 21 dogs with NHL and 7 normal dogs, the methylation status of the promoter CpG island of the gene was defined for each sample using combined bisulfite restriction analysis (COBRA), methylation-specific PCR (MSP), and bisulfite sequencing methods. Relative gene expression was determined using real-time PCR. Results: The mRNA sequence of canine DLC1 is highly similar to the human orthologue and contains all protein functional groups, with 97% or greater similarity in functional regions. Hypermethylation of the 5' and 3' flanking regions of the promoter was statistically significantly associated with the NHL phenotype, but was not associated with silencing of expression or differences in survival. Conclusions: The canine DLC1 is constructed highly similarly to the human gene, which has been shown to be an important tumor suppressor in many forms of cancer. As in human NHL, the promoter CpG island of DLC1 in canine NHL samples is abnormally hypermethylated, relative to normal lymphoid tissue. This study confirms that hypermethylation occurs in canine cancers, further supporting the use of companion dogs as comparative models of disease for evaluation of carcinogenesis, biomarker diagnosis, and therapy.

Composite likelihood estimation of demographic parameters

Daniel Garrigan — 2009-11-12T00:00:00Z

Background: Most existing likelihood-based methods for fitting historical demographic models to DNA sequence polymorphism data to do not scale feasibly up to the level of whole-genome data sets. Computational economies can be achieved by incorporating two forms of pseudo-likelihood: composite and approximatelikelihood methods. Composite likelihood enables scaling up to large data sets because it takes the product of marginal likelihoods as an estimator of the likelihood of the complete data set. This approach is especially useful when a large number of genomic regions constitutes the data set. Additionally, approximate likelihood methods can reduce the dimensionality of the data by summarizing the information in the original data by either a sufficient statistic, or a set of statistics. Both composite and approximate likelihood methods hold promise for analyzing large data sets or for use in situations where the underlying demographic model is complex and has many parameters. This paper considers a simple demographic model of allopatric divergence between two populations, in which one of the population is hypothesized to have experienced a founder event, or population bottleneck. A large resequencing data set from human populations is summarized by the joint frequency spectrum, which is a matrix of the genomic frequency spectrum of derivedbase frequencies in two populations. A Bayesian Metropolis-coupled Markov chain Monte Carlo (MCMCMC) method for parameterestimation is developed that uses both composite and likelihood methods and is applied to the three different pairwise combinations of the human population resequence data. The accuracy of the method is also tested on data sets sampled from a simulated population model with known parameters. Results: The Bayesian MCMCMC method also estimates the ratio of effective population size for the X chromosome versus that of the autosomes. The method is shown to estimate, with reasonable accuracy, demographic parameters from 12 simulated data sets that vary in the magnitude of a founder event and a skew in the effective population size of the X chromosome relative to the autosomes. The behavior of the Markov chain is also examined and shown to convergence to its stationary distribution, while also showing high levels of parameter mixing. The analysis of three pairwise comparisons of sub-Saharan African human populations with non-African human populations do not provide unequivocal support for a strong non-African founder event from these nuclear data. The estimates do however suggest a skew in the ratio of X chromosome to autosome effective population size that is greater than one. However in all three cases, the 95% highest posterior density interval for this ratio does include three-fourths, the value expected under an equal breeding sex ratio. Conclusions: The implementation of composite and approximate likelihood methods in a framework that includes MCMCMC demographic parameter estimation shows great promise for being flexible and computationally efficient enough to scale up to the level of whole-genome polymorphism and divergence analysis. Further work must be done to characterize the effects of the assumption of linkage equilibrium among genomic regions that is crucial to the validity of applying the composite likelihood method.

Variation in genetic admixture and population structure among Latinos - The Los Angeles Latino Eye Study (LALES)

Corina Shtir — 2009-11-10T00:00:00Z

Background: Population structure and admixture have strong confounding effects on genetic association studies. Discordant frequencies for age-related macular degeneration (AMD) risk alleles and for AMD incidence and prevalence rates are reported across different ethnic groups. We examined the genomic ancestry characterizing 538 Latinos drawn from the Los Angeles Latino Eye Study [LALES] as part of an ongoing AMD-association study. To help assess the degree of Native American ancestry inherited by Latino populations we sampled 25 Mayans and 5 Mexican Indians collected through Coriell's Institute. Levels of European, Asian, and African descent in Latinos were inferred through the USC Multiethnic Panel (USC MEP), formed from a sample from the Multiethnic Cohort (MEC) study, the Yoruba African samples from HapMap II, the Singapore Chinese Health Study, and a prospective cohort from Shanghai, China. A total of 233 ancestry informative markers were genotyped for 538 LALES Latinos, 30 Native Americans, and 355 USC MEP individuals (African Americans, Japanese, Chinese, European Americans, Latinos, and Native Hawaiians). Sensitivity of ancestry estimates to relative sample size was considered. Results: We detected strong evidence for recent population admixture in LALES Latinos. Gradients of increasing Native American background and of correspondingly decreasing European ancestry were observed as a function of birth origin from North to South. The strongest excess of homozygosity, a reflection of recent population admixture, was observed in non-US born Latinos that recently populated the US. A set of 42 SNPs especially informative for distinguishing between Native Americans and Europeans were identified. Conclusion: These findings reflect the historic migration patterns of Native Americans and suggest that while the 'Latino' label is used to categorize the entire population, there exists a strong degree of heterogeneity within that population, and that it will be important to assess this heterogeneity within future association studies on Latino populations. Our study raises awareness of the diversity within "Latinos" and the necessity to assess appropriate risk and treatment management.

Optimization of selection contribution and mate allocations in monoecious tree breeding populations

Jon Hallander — 2009-11-06T00:00:00Z

Background: The combination of optimized contribution dynamic selection and various mating schemes was investigated over seven generations for a typical tree breeding scenario. The allocation of mates was optimized using a simulated annealing algorithm for various object functions including random mating (RM), positive assortative mating (PAM) and minimization of pair-wise coancestry between mates (MCM) all combined with minimization of variance in family size and coancestry. The present study considered two levels of heritability (0.05 and 0.25), two restrictions on relatedness (group coancestry; 1 and 2%) and two maximum permissible numbers of crosses in each generation (100 and 400). The infinitesimal genetic model was used to simulate the genetic architecture of the trait that was the subject of selection. A framework of the long term genetic contribution of ancestors was used to examine the impacts of the mating schemes on population parameters. Results: MCM schemes produced on average, an increased rate of genetic gain in the breeding population, although the difference between schemes was small but significant after seven generations (up to 7.1% more than obtained with RM). In addition, MCM reduced the level of inbreeding by as much as 37% compared with RM, although the rate of inbreeding was similar after three generations of selection. PAM schemes yielded levels of genetic gain similar to those produced by RM, but the increase in the level of inbreeding was substantial (up to 43%). Conclusion: The main reason why MCM schemes yielded higher genetic gains was the improvement in managing the long term genetic contribution of founders in the population; this was achieved by connecting unrelated families. In addition, the accumulation of inbreeding was reduced by MCM schemes since the variance in long term genetic contributions of founders was smaller than in the other schemes. Consequently, by combining an MCM scheme with an algorithm that optimizes contributions of the selected individuals, a higher long term response is obtained while reducing the risk within the breeding program.

Developing a set of ancestry-sensitive DNA markers reflecting continental origins of humans.

Paula Kersbergen — 2009-10-27T00:00:00Z

Background: The identification and use of Ancestry-Sensitive Markers (ASMs), i.e. genetic polymorphisms facilitating the genetic reconstruction of geographical origins of individuals, is far from straightforward. Results: Here we describe the ascertainment and application of five different sets of 47 single nucleotide polymorphisms (SNPs) allowing the inference of major human groups of different continental origin. For this, we first used 74 cell lines, representing human males from six different geographical areas and screened them with the Affymetrix Mapping 10K assay. In addition to using summary statistics estimating the genetic diversity among multiple groups of individuals defined by geography or language, we also used the program STRUCTURE to detect genetically distinct subgroups. Subsequently, we used a pairwise FST ranking procedure among all pairs of genetic subgroups in order to identify a single best performing set of ASMs. Our initial results were independently confirmed by genotyping this set of ASMs in 22 individuals from Somalia, Afghanistan and Sudan and in 919 samples from the CEPH Human Genome Diversity Panel (HGDP-CEPH) Conclusion: By means of our pairwise population FST ranking approach we identified a set of 47 SNPs that could serve as a panel of ASMs at a continental level.

Microsatellite markers of water buffalo, Bubalus bubalis - development, characterisation and linkage disequilibrium studies

Muniyandi Nagarajan — 2009-10-21T00:00:00Z

Background: Microsatellite markers are highly polymorphic and widely used in genome mapping and population genetic studies in livestock species. River buffalo, Bubalus bubalis is an economically important livestock species, though only a limited number of microsatellite markers have been reported thus far in this species. Results: In the present study, using two different approaches 571 microsatellite markers have been characterized for water buffalo. Of the 571 microsatellite markers, 498 were polymorphic with average heterozygosity of 0.51 on a panel of 24 unrelated buffalo. Fisher exact test was used to detect LD between the marker pairs. Among the 137550 pairs of marker combination, 14.58% pairs showed significant LD (P < 0.05). Further to check the suitability of these microsatellite markers to map these on a radiation hybrid map of buffalo genome, the markers were tested on Chinese hamster genomic DNA for amplification. Only seven of these markers showed amplification in Chinese hamster, and thus 564, of these can be added to the radiation hybrid map of this species. Conclusion: The high conservation of cattle microsatellite loci in water buffalo promises the usefulness of the cattle microsatellites markers on buffalo. The polymorphic markers characterised in this study will contribute to genetic linkage and radiation hybrid mapping of water buffalo and population genetic studies.

Large introns in relation to alternative splicing and gene evolution: a case study of Drosophila bruno-3

Nikolai Kandul — 2009-10-19T00:00:00Z

Background: Alternative splicing (AS) of maturing mRNA can generate structurally and functionally distinct transcripts from the same gene. Recent bioinformatic analyses of available genome databases inferred a positive correlation between intron length and AS. To study the interplay between intron length and AS empirically and in more detail, we analyzed the diversity of alternatively spliced transcripts (ASTs) in the Drosophila RNA-binding Bruno-3 (Bru-3) gene. This gene was known to encode thirteen exons separated by introns of diverse sizes, ranging from 71 to 41,973 nucleotides in D. melanogaster. Although Bru-3's structure is expected to be conducive to AS, only two ASTs of this gene were previously described. Results: Cloning of RT-PCR products of the entire ORF from four species representing three diverged Drosophila lineages provided an evolutionary perspective, high sensitivity, and long-range contiguity of splice choices currently unattainable by high-throughput methods. Consequently, we identified three new exons, a new exon fragment and thirty-three previously unknown ASTs of Bru-3. All exon-skipping events in the gene were mapped to the exons surrounded by introns of at least 800 nucleotides, whereas exons split by introns of less than 250 nucleotides were always spliced contiguously in mRNA. Cases of exon loss and creation during Bru-3 evolution in Drosophila were also localized within large introns. Notably, we identified a true de novo exon gain: exon 8 was created along the lineage of the obscura group from intronic sequence between cryptic splice sites conserved among all Drosophila species surveyed. Exon 8 was included in mature mRNA by the species representing all the major branches of the obscura group. To our knowledge, the origin of exon 8 is the first documented case of exonization of intronic sequence outside vertebrates. Conclusion: We found that large introns can promote AS via exon-skipping and exon turnover during evolution likely due to frequent errors in their removal from maturing mRNA. Large introns could be a reservoir of genetic diversity, because they have a greater number of mutable sites than short introns. Taken together, gene structure can constrain and/or promote gene evolution.

Automated DNA mutation detection using universal conditions direct sequencing: application to ten muscular dystrophy genes

Richard Bennett — 2009-10-18T00:00:00Z

Background: One of the most common and efficient methods for detecting mutations in genes is PCR amplification followed by direct sequencing. Until recently, the process of designing PCR assays has been to focus on individual assay parameters rather than concentrating on matching conditions for a set of assays. Primers for each individual assay were selected based on location and sequence concerns. The two primer sequences were then iteratively adjusted to make the individual assays work properly. This generally resulted in groups of assays with different annealing temperatures that required the use of multiple thermal cyclers or multiple passes in a single thermal cycler making diagnostic testing time-consuming, laborious and expensive.These factors have severely hampered diagnostic testing services, leaving many families without an answer for the exact cause of a familial genetic disease. A search of GeneTests for sequencing analysis of the entire coding sequence for genes that are known to cause muscular dystrophies returns only a small list of laboratories that perform comprehensive gene panels.The hypothesis for the study was that a complete set of universal assays can be designed to amplify and sequence any gene or family of genes using computer aided design tools. If true, this would allow automation and optimization of the mutation detection process resulting in reduced cost and increased throughput. Results: An automated process has been developed for the detection of deletions, duplications/insertions and point mutations in any gene or family of genes and has been applied to ten genes known to bear mutations that cause muscular dystrophy: DMD; CAV3; CAPN3; FKRP; TRIM32; LMNA; SGCA; SGCB; SGCG; SGCD. Using this process, mutations have been found in five DMD patients and four LGMD patients (one in the FKRP gene, one in the CAV3 gene, and two likely causative heterozygous pairs of variations in the CAPN3 gene of two other patients). Methods and assay sequences are reported in this paper. Conclusion: This automated process allows laboratories to discover DNA variations in a short time and at low cost.

Population genetics of the understory fishtail palm Chamaedorea ernesti-augusti in Belize: high genetic connectivity with local differentiation

Angelica Cibrian-Jaramillo — 2009-10-09T00:00:00Z

Background: Developing a greater understanding of population genetic structure in lowland tropical plant species is highly relevant to our knowledge of increasingly fragmented forests and to the conservation of threatened species. Specific studies are particularly needed for taxa whose population dynamics are further impacted by human harvesting practices. One such case is the fishtail or xaté palm (Chamaedorea ernesti-augusti) of Central America, whose wild-collected leaves are becoming progressively more important to the global ornamental industry. We use microsatellite markers to describe the population genetics of this species in Belize and test the effects of climate change and deforestation on its recent and historical effective population size. Results: We found high levels of inbreeding coupled with moderate or high allelic diversity within populations. Overall high gene flow was observed, with a north and south gradient and ongoing differentiation at smaller spatial scales. Immigration rates among populations were more difficult to discern, with minimal evidence for isolation by distance. We infer a tenfold reduction in effective population size ca. 10,000 years ago, but fail to detect changes attributable to Mayan or contemporary deforestation. Conclusion: Populations of C. ernesti-augusti are genetically heterogeneous demes at a local spatial scale, but are widely connected at a regional level in Belize. We suggest that the inferred patterns in population genetic structure are the result of the colonization of this species into Belize following expansion of humid forests in combination with demographic and mating patterns. Within populations, we hypothesize that low aggregated population density over large areas, short distance pollen dispersal via thrips, low adult survival, and low fruiting combined with early flowering may contribute towards local inbreeding via genetic drift. Relatively high levels of regional connectivity are likely the result of animal-mediated long-distance seed dispersal. The greatest present threat to the species is the potential onset of inbreeding depression as the result of increased human harvesting activities. Future genetic studies in understory palms should focus on both fine-scale and landscape-level genetic structure.

Molecular cloning and characterization of the porcine prostaglandin transporter (SLCO2A1): evaluation of its role in F4 mediated neonatal diarrhoea

Mario Van Poucke — 2009-10-06T00:00:00Z

Background: Because prostaglandins are involved in many (patho)physiological processes, SLCO2A1 was already characterized in several species in an attempt to unravel specific processes/deficiencies. Here, we describe the molecular cloning and characterization of the porcine ortholog in order to evaluate its possible involvement in F4 enterotoxigenic E. coli mediated neonatal diarrhoea, based on a positional candidate gene approach study. Results: Porcine SLCO2A1 is organized in 14 exons, containing an open reading frame of 1935 bp, encoding a 12-transmembrane organic anion cell surface transporter of 644 aa. The -388 to -5 upstream region comprises a (CpG)48 island containing a number of conserved promoter elements, including a TATA box. A potential alternative promoter region was found in the conserved -973 to -700 upstream region. No consensus polyadenylation signal was discovered in the 3' UTR. Repeat sequences were found in 15% of all the non coding sequences.As expected for a multifunctional protein, a wide tissue distribution was observed. mRNA expression was found in the adrenal gland, bladder, caecum, colon (centripetal coil/centrifugal coil), diaphragm, duodenum, gallbladder, heart, ileum, jejunum, kidney, liver, longissimus dorsi muscle, lung, lymph node, mesenterium, rectum, spleen, stomach, tongue and ureter, but not in the aorta, oesophagus and pancreas.The promoter region and the exons (including the splice sites) of SLCO2A1 were resequenced in 5 F4ab/ac receptor positive and 5 F4ab/ac receptor negative pigs. Two silent and 2 missense (both S → L at position 360 and 633) mutations were found, but none was associated with the F4ab/ac receptor phenotype. In addition, no phenotype associated differential mRNA expression or alternative/abberant splicing/polyadenylation was found in the jejunum. Conclusion: The molecular cloning and characterization of porcine SLCO2A1 not only contributes to the already existing knowledge about the transporter in general, but enables studies on porcine prostaglandin related processes/deficiencies as patient and/or model. Here we examined its possible involvement as receptor in F4 enterotoxigenic E. coli mediated neonatal diarrhoea. Because no phenotype associated differences could be found in the gene sequence nor in its jejunal transcription profile of F4ab/ac receptor positive/negative pigs, SLCO2A1 can most likely be excluded as receptor for F4 bacteria.