DnaSP v5: A software for comprehensive analysis of DNA polymorphism data

DnaSP is a software package for a comprehensive analysis of DNA polymorphism data. Version 5 implements a number of new features and analytical methods allowing extensive DNA polymorphism analyses on large datasets. Among other features, the newly implemented methods allow for: (i) analyses on multiple data files; (ii) haplotype phasing; (iii) analyses on insertion/deletion polymorphism data; (iv) visualizing sliding window results integrated with available genome annotations in the UCSC browser.

DnaSP Version 5. DNA Sequence Polymorphism

DnaSP, DNA Sequence Polymorphism, is a software package for the analysis of nucleotide polymorphism from aligned DNA sequence data. DnaSP can estimate several measures of DNA sequence variation within and between populations (in noncoding, synonymous or nonsynonymous sites, or in various sorts of codon positions), as well as linkage disequilibrium, recombination, gene flow and gene conversion parameters. DnaSP can also carry out several tests of neutrality: Hudson, Kreitman and Aguadé (1987), Tajima (1989), McDonald and Kreitman (1991), Fu and Li (1993), and Fu (1997) tests. Additionally, DnaSP can estimate the confidence intervals of some test-statistics by the coalescent. The results of the analyses are displayed on tabular and graphic form.

VariScan Software

VariScan is a software package for the analysis of DNA sequence polymorphisms at the whole genome scale. Among other features, the software:(1) can conduct many population genetic analyses; (2) incorporates a multiresolution wavelet transform-based method that allows capturing relevant information from DNA polymorphism data; and (3) it facilitates the visualization of the results in the most commonly used genome browsers.

Vanishing chromosomal inversion clines in Drosophila subobscura from Chile: Is behavioral thermoregulation to blame?

Chromosomal inversion clines paralleling the long-standing ones in native Palearctic populations of Drosophila subobscura evolved swiftly after this species invaded the Americas in the late 1970s and early 1980s. However, the new clines did not consistently continue to converge on the Old World baseline. Our recent survey of Chilean populations of D. subobscura shows that inversion clines have faded or even changed sign with latitude. Here, we investigate the hypothesis that this fading of inversion clines might be due to the Bogert effect, namely, that flies' thermoregulatory behavior has eventually compensated for environmental variation in temperature, thus buffering selection on thermal-related traits. We show that latitudinal divergence in thermal preference (T-p) has evolved in Chile for females, with higher-latitude flies having a lower mean T-p. Plastic responses in T-p also lessen latitudinal thermal variation because flies developed at colder temperatures prefer warmer microclimates. Our results are consistent with the idea that active behavioral thermoregulation might buffer environmental variation and reduce the potential effect of thermal selection on other traits as chromosomal arrangements.

Scrfl restriction fragment length polymorphism at the D7S23 locus (probe pKM.19), closely linked to cystic fibrosis

Source/Description: pKM.19 is a 1.0 kb EcoRI genomic fragment in pUC13 (ref. 1,2). pPl was isolated independently but contains the same fragment as pKM.19 (ref. 3)...

Inversion polymorphism in two Serbian natural populations of Drosophila subobscura: Analysis of long-term changes

To study whether inversions (or arrangements) by themselves or karyotypes are the main global warming adaptive target of natural selection, two Drosophila subobscura Serbian populations (Apatin and Petnica) were re analyzed using different statistical approaches. Both populations were sampled in an approximately 15 years period: Apatin in 1994 and 2008 + 2009 and Petnica in 1995 and 2010. For all chromosomes, the four collections studied were in Hardy-Weinberg equilibrium. Thus, it seems that inversions (or arrangements) combined at random to constitute populations" karyotypes. However, there were differences in karyotypic fre quencies along the years, although they were significant only for Apatin population. It is possible to conclude that inversions (or arrangements) are likely the target of natural selection, because they presented long term changes, but combine at random to generate the corresponding karyotypic combinations. As a consequence, the frequencies of karyotypes also change along time.

A polymorphic DNA probe from chromosome 7 (7q22)

Source/Description: p6a-l is a O.9 kb HindIII/BamHl genomic fragment subclone or cosmic cNX.6a in pUC13. cNX.6a was isolated from a non-methylated enriched library from the CMGT cell line Cll (1,2).

Chromosomal polymorphism of D. subobscura: no differences between wild males and sons of wild females

When analyzing the chromosomal polymorphism of D. subobscura natural populations it is assumed that the information provided by wild males and sons of wild females is equivalent. Thus, using both in the analysis it is possible to increase the sample size. However, it is important to verify whether there are significant differences between both groups or not. The aim of this research has been to statistically compare the results of chromosomal polymorphism of both groups. We have used data from Avala Mountain (Serbia) where D. subobscura flies were collected from the 30th May to the 5th June 2011. Avala is located 18 km south of Belgrade and the trapping place is a forest with polydominant communities of Fagetum submontanum Table 1. Number and percentage of adult flies collected in Font Groga (Barcelona, Spain) on 9th October 2013. Males and sons of wild females were crossed with virgin females of the Küsnacht strain. Third instar larvae from F1 were dissected to obtain the salivary glands and the polytene chromosomes were stained and squashed in aceto-orcein solution. No significant differences were observed for any chromosome of the karyotype: A (p-value = 0.485), J (p-value = 0.230), U (p-value =0.572), E (p-value = 0.536), and O (p-value = 0.338). Thus, it seems that the two groups can be grouped together to obtain the chromosomal polymorphism of the population.

Study of linkage between miniature and singed genes in Drosophila melanogaster

We have developed a practical exercise for undergraduate students whose main aim is to identify, using genetic crosses, a pair of D. melanogaster mutations (miniature and singed). Each student receives a vial with the problem strain containing two unknown mutations. The first step is to observe and describe both mutations. Then, the students carry out genetic crosses between mutant and normal strains: (P) ♀ mutant strain × ♂ normal strain (P) ♀ normal strain × ♂ mutant strain A different offspring is expected in these crosses: in the first one we will obtain normal females and m sn males, whereas in the second all individuals will present normal phenotype. It is possible to deduce that both are sex linked mutations. With this information and to simplify the amount of work, only F1 individuals from the first cross will be used (m+sn+ / m sn × m sn / Y chrom.) to obtain the F2 generation. By counting the number of miniature (recombinant type), singed (recombinant type), miniature-singed (parental type) and normal (parental type) flies it is possible to estimate the recombination frequency between both genes. Knowing the phenotype, their chromosomal location (X chromosome) and the genetic distance between both mutations, it is possible to identify them by finding all this information in a Drosophila melanogaster genetic map. Additionally, a statistical analysis can be carried out to compare the number of expected F2 individuals with those observed in the experiment. As the distance between both genes is 15.1 m.u., then the expected percentages for each phenotype would be: normal (42.45%), miniature-signed (42.45%), miniature (7.55%) and singed (7.55%). Multiplying the frequency of each class by the total number of individuals obtained in the F2 it is possible to estimate the expected number of flies for each class. Finally, a χ2 test can be computed to ascertain whether there are significant differences between expected and observed number of individuals.

Genetic origin, admixture, and asymmetry in maternal and paternal human lineages in Cuba

Background: Before the arrival of Europeans to Cuba, the island was inhabited by two Native American groups, the Tainos and the Ciboneys. Most of the present archaeological, linguistic and ancient DNA evidence indicates a South American origin for these populations. In colonial times, Cuban Native American people were replaced by European settlers and slaves from Africa. It is still unknown however, to what extent their genetic pool intermingled with and was 'diluted' by the arrival of newcomers. In order to investigate the demographic processes that gave rise to the current Cuban population, we analyzed the hypervariable region I (HVS-I) and five single nucleotide polymorphisms (SNPs) in the mitochondrial DNA (mtDNA) coding region in 245 individuals, and 40 Y-chromosome SNPs in 132 male individuals. Results: The Native American contribution to present-day Cubans accounted for 33% of the maternal lineages, whereas Africa and Eurasia contributed 45% and 22% of the lineages, respectively. This Native American substrate in Cuba cannot be traced back to a single origin within the American continent, as previously suggested by ancient DNA analyses. Strikingly, no Native American lineages were found for the Y-chromosome, for which the Eurasian and African contributions were around 80% and 20%, respectively. Conclusion: While the ancestral Native American substrate is still appreciable in the maternal lineages, the extensive process of population admixture in Cuba has left no trace of the paternal Native American lineages, mirroring the strong sexual bias in the admixture processes taking place during colonial times.

Satisfaction survey with DNA cards method to collect genetic samples for pharmacogenetics studies

Background: Pharmacogenetic studies are essential in understanding the interindividual variability of drug responses. DNA sample collection for genotyping is a critical step in genetic studies. A method using dried blood samples from finger-puncture, collected on DNA-cards, has been described as an alternative to the usual venepuncture technique. The purpose of this study is to evaluate the implementation of the DNA cards method in a multicentre clinical trial, and to assess the degree of investigators' satisfaction and the acceptance of the patients perceived by the investigators.Methods: Blood samples were collected on DNA-cards. The quality and quantity of DNA recovered were analyzed. Investigators were questioned regarding their general interest, previous experience, safety issues, preferences and perceived patient satisfaction. Results: 151 patients' blood samples were collected. Genotyping of GST polymorphisms was achieved in all samples (100%). 28 investigators completed the survey. Investigators perceived patient satisfaction as very good (60.7%) or good (39.3%), without reluctance to finger puncture. Investigators preferred this method, which was considered safer and better than the usual methods. All investigators would recommend using it in future genetic studies. Conclusion: Within the clinical trial setting, the DNA-cards method was very well accepted by investigators and patients (in perception of investigators), and was preferred to conventional methods due to its ease of use and safety.

Association of MMP-2 polymorphisms with severe and very severe COPD: a case control study of MMPs-1, 9 and 12 in a European population.

BACKGROUND: Genetic factors play a role in chronic obstructive pulmonary disease (COPD) but are poorly understood. A number of candidate genes have been proposed on the basis of the pathogenesis of COPD. These include the matrix metalloproteinase (MMP) genes which play a role in tissue remodelling and fit in with the protease--antiprotease imbalance theory for the cause of COPD. Previous genetic studies of MMPs in COPD have had inadequate coverage of the genes, and have reported conflicting associations of both single nucleotide polymorphisms (SNPs) and SNP haplotypes, plausibly due to under-powered studies. METHODS: To address these issues we genotyped 26 SNPs, providing comprehensive coverage of reported SNP variation, in MMPs- 1, 9 and 12 from 977 COPD patients and 876 non-diseased smokers of European descent and evaluated their association with disease singly and in haplotype combinations. We used logistic regression to adjust for age, gender, centre and smoking history. RESULTS: Haplotypes of two SNPs in MMP-12 (rs652438 and rs2276109), showed an association with severe/very severe disease, corresponding to GOLD Stages III and IV. CONCLUSIONS: Those with the common A-A haplotype for these two SNPs were at greater risk of developing severe/very severe disease (p = 0.0039) while possession of the minor G variants at either SNP locus had a protective effect (adjusted odds ratio of 0.76; 95% CI 0.61 - 0.94). The A-A haplotype was also associated with significantly lower predicted FEV1 (42.62% versus 44.79%; p = 0.0129). This implicates haplotypes of MMP-12 as modifiers of disease severity.

Segregation of regulatory polymorphisms with effects on the gluteus medius transcriptome in a purebred pig population

Background: The main goal of the present study was to analyse the genetic architecture of mRNA expression in muscle, a tissue with an outmost economic importance for pig breeders. Previous studies have used F2 crosses to detect porcine expression QTL (eQTL), so they contributed with data that mostly represents the between-breed component of eQTL variation. Herewith, we have analysed eQTL segregation in an outbred Duroc population using two groups of animals with divergent fatness profiles. This approach is particularly suitable to analyse the within-breed component of eQTL variation, with a special emphasis on loci involved in lipid metabolism. Methodology/Principal Findings: GeneChip Porcine Genome arrays (Affymetrix) were used to determine the mRNA expression levels of gluteus medius samples from 105 Duroc barrows. A whole-genome eQTL scan was carried out with a panel of 116 microsatellites. Results allowed us to detect 613 genome-wide significant eQTL unevenly distributed across the pig genome. A clear predominance of trans- over cis-eQTL, was observed. Moreover, 11 trans-regulatory hotspots affecting the expression levels of four to 16 genes were identified. A Gene Ontology study showed that regulatory polymorphisms affected the expression of muscle development and lipid metabolism genes. A number of positional concordances between eQTL and lipid trait QTL were also found, whereas limited evidence of a linear relationship between muscle fat deposition and mRNA levels of eQTL regulated genes was obtained. Conclusions/Significance: Our data provide substantial evidence that there is a remarkable amount of within-breed genetic variation affecting muscle mRNA expression. Most of this variation acts in trans and influences biological processes related with muscle development, lipid deposition and energy balance. The identification of the underlying causal mutations and the ascertainment of their effects on phenotypes would allow gaining a fundamental perspective about how complex traits are built at the molecular level.

Gene Set Analysis for improving genetic association studies

Introduction. Genetic epidemiology is focused on the study of the genetic causes that determine health and diseases in populations. To achieve this goal a common strategy is to explore differences in genetic variability between diseased and nondiseased individuals. Usual markers of genetic variability are single nucleotide polymorphisms (SNPs) which are changes in just one base in the genome. The usual statistical approach in genetic epidemiology study is a marginal analysis, where each SNP is analyzed separately for association with the phenotype. Motivation. It has been observed, that for common diseases the single-SNP analysis is not very powerful for detecting genetic causing variants. In this work, we consider Gene Set Analysis (GSA) as an alternative to standard marginal association approaches. GSA aims to assess the overall association of a set of genetic variants with a phenotype and has the potential to detect subtle effects of variants in a gene or a pathway that might be missed when assessed individually. Objective. We present a new optimized implementation of a pair of gene set analysis methodologies for analyze the individual evidence of SNPs in biological pathways. We perform a simulation study for exploring the power of the proposed methodologies in a set of scenarios with different number of causal SNPs under different effect sizes. In addition, we compare the results with the usual single-SNP analysis method. Moreover, we show the advantage of using the proposed gene set approaches in the context of an Alzheimer disease case-control study where we explore the Reelin signal pathway.