Genetic Variation In Individuals Research Articles

De novo inference of stratification and local admixture in sequencing studies

Analysis of population structures and genome local ancestry has become increasingly important in population and disease genetics. With the advance of next generation sequencing technologies, complete genetic variants in individuals' genomes are quickly generated, providing unprecedented opportunities for learning population evolution histories and identifying local genetic signatures at the SNP resolution. The successes of those studies critically rely on accurate and powerful computational tools that can fully utilize the sequencing information. Although many algorithms have been developed for population structure inference and admixture mapping, many of them only work for independent SNPs in genotype or haplotype format, and require a large panel of reference individuals. In this paper, we propose a novel probabilistic method for detecting population structure and local admixture. The method takes input of sequencing data, genotype data and haplotype data. The method characterizes the dependence of genetic variants via haplotype segmentation, such that all variants detected in a sequencing study can be fully utilized for inference. The method further utilizes a infinite-state Bayesian Markov model to perform de novo stratification and admixture inference. Using simulated datasets from HapMapII and 1000Genomes, we show that our method performs superior than several existing algorithms, particularly when limited or no reference individuals are available. Our method is applicable to not only human studies but also studies of other species of interests, for which little reference information is available.Software Availability: http://stat.psu.edu/~yuzhang/software/dbm.tar

Genes Contributing to Pain Sensitivity in the Normal Population: An Exome Sequencing Study

Sensitivity to pain varies considerably between individuals and is known to be heritable. Increased sensitivity to experimental pain is a risk factor for developing chronic pain, a common and debilitating but poorly understood symptom. To understand mechanisms underlying pain sensitivity and to search for rare gene variants (MAF<5%) influencing pain sensitivity, we explored the genetic variation in individuals' responses to experimental pain. Quantitative sensory testing to heat pain was performed in 2,500 volunteers from TwinsUK (TUK): exome sequencing to a depth of 70× was carried out on DNA from singletons at the high and low ends of the heat pain sensitivity distribution in two separate subsamples. Thus in TUK1, 101 pain-sensitive and 102 pain-insensitive were examined, while in TUK2 there were 114 and 96 individuals respectively. A combination of methods was used to test the association between rare variants and pain sensitivity, and the function of the genes identified was explored using network analysis. Using causal reasoning analysis on the genes with different patterns of SNVs by pain sensitivity status, we observed a significant enrichment of variants in genes of the angiotensin pathway (Bonferroni corrected p = 3.8×10−4). This pathway is already implicated in animal models and human studies of pain, supporting the notion that it may provide fruitful new targets in pain management. The approach of sequencing extreme exome variation in normal individuals has provided important insights into gene networks mediating pain sensitivity in humans and will be applicable to other common complex traits.

Bias, accuracy, and impact of indirect genetic effects in infectious diseases.

Selection for improved host response to infectious disease offers a desirable alternative to chemical treatment but has proven difficult in practice, due to low heritability estimates of disease traits. Disease data from field studies is often binary, indicating whether an individual has become infected or not following exposure to an infectious disease. Numerous studies have shown that from this data one can infer genetic variation in individuals’ underlying susceptibility. In a previous study, we showed that with an indirect genetic effect (IGE) model it is possible to capture some genetic variation in infectivity, if present, as well as in susceptibility. Infectivity is the propensity of transmitting infection upon contact with a susceptible individual. It is an important factor determining the severity of an epidemic. However, there are severe shortcomings with the Standard IGE models as they do not accommodate the dynamic nature of disease data. Here we adjust the Standard IGE model to (1) make expression of infectivity dependent on the individuals’ disease status (Case Model) and (2) to include timing of infection (Case-ordered Model). The models are evaluated by comparing impact of selection, bias, and accuracy of each model using simulated binary disease data. These were generated for populations with known variation in susceptibility and infectivity thus allowing comparisons between estimated and true breeding values. Overall the Case Model provided better estimates for host genetic susceptibility and infectivity compared to the Standard Model in terms of bias, impact, and accuracy. Furthermore, these estimates were strongly influenced by epidemiological characteristics. However, surprisingly, the Case-Ordered model performed considerably worse than the Standard and the Case Models, pointing toward limitations in incorporating disease dynamics into conventional variance component estimation methodology and software used in animal breeding.

Population structure analysis using rare and common functional variants

Next-generation sequencing technologies now make it possible to genotype and measure hundreds of thousands of rare genetic variations in individuals across the genome. Characterization of high-density genetic variation facilitates control of population genetic structure on a finer scale before large-scale genotyping in disease genetics studies. Population structure is a well-known, prevalent, and important factor in common variant genetic studies, but its relevance in rare variants is unclear. We perform an extensive population structure analysis using common and rare functional variants from the Genetic Analysis Workshop 17 mini-exome sequence. The analysis based on common functional variants required 388 principal components to account for 90% of the variation in population structure. However, an analysis based on rare variants required 532 significant principal components to account for similar levels of variation. Using rare variants, we detected fine-scale substructure beyond the population structure identified using common functional variants. Our results show that the level of population structure embedded in rare variant data is different from the level embedded in common variant data and that correcting for population structure is only as good as the level one wishes to correct.

GENOME-WIDE INTERACTION STUDY OF CUMULATIVE PARTICULATE MATTER EXPOSURE ON AGE-RELATED LUNG FUNCTION DECLINE IN ADULTS, THE SAPALDIA COHORT.

Background and Aims: Air pollution exposure accelerates age-related lung function decline. There is however evidence that i) genetic variation in individuals determines the susceptibility to ambien...

Inter-chromosomal variation in the pattern of human population genetic structure

Emerging technologies now make it possible to genotype hundreds of thousands of genetic variations in individuals, across the genome. The study of loci at finer scales will facilitate the understanding of genetic variation at genomic and geographic levels. We examined global and chromosomal variations across HapMap populations using 3.7 million single nucleotide polymorphisms to search for the most stratified genomic regions of human populations and linked these regions to ontological annotation and functional network analysis. To achieve this, we used five complementary statistical and genetic network procedures: principal component (PC), cluster, discriminant, fixation index (FST) and network/pathway analyses. At the global level, the first two PC scores were sufficient to account for major population structure; however, chromosomal level analysis detected subtle forms of population structure within continental populations, and as many as 31 PCs were required to classify individuals into homogeneous groups. Using recommended population ancestry differentiation measures, a total of 126 regions of the genome were catalogued. Gene ontology and networks analyses revealed that these regions included the genes encoding oculocutaneous albinism II (OCA2), hect domain and RLD 2 (HERC2), ectodysplasin A receptor (EDAR) and solute carrier family 45, member 2 (SLC45A2). These genes are associated with melanin production, which is involved in the development of skin and hair colour, skin cancer and eye pigmentation. We also identified the genes encoding interferon-γ (IFNG) and death-associated protein kinase 1 (DAPK1), which are associated with cell death, inflammatory and immunological diseases. An in-depth understanding of these genomic regions may help to explain variations in adaptation to different environments. Our approach offers a comprehensive strategy for analysing chromosome-based population structure and differentiation, and demonstrates the application of complementary statistical and functional network analysis in human genetic variation studies.

Population-based study of genetic variation in individuals with autism spectrum disorders from Croatia

BackgroundGenome-wide studies on autism spectrum disorders (ASDs) have mostly focused on large-scale population samples, but examination of rare variations in isolated populations may provide additional insights into the disease pathogenesis.MethodsAs a first step in the genetic analysis of ASD in Croatia, we characterized genetic variation in a sample of 103 subjects with ASD and 203 control individuals, who were genotyped using the Illumina HumanHap550 BeadChip. We analyzed the genetic diversity of the Croatian population and its relationship to other populations, the degree of relatedness via Runs of Homozygosity (ROHs), and the distribution of large (>500 Kb) copy number variations.ResultsCombining the Croatian cohort with several previously published populations in the FastME analysis (an alternative to Neighbor Joining) revealed that Croatian subjects cluster, as expected, with Southern Europeans; in addition, individuals from the same geographic region within Europe cluster together. Whereas Croatian subjects could be separated from a sample of healthy control subjects of European origin from North America, Croatian ASD cases and controls are well mixed. A comparison of runs of homozygosity indicated that the number and the median length of regions of homozygosity are higher for ASD subjects than for controls (p = 6 × 10-3). Furthermore, analysis of copy number variants found a higher frequency of large chromosomal rearrangements (>2 Mb) in ASD cases (5/103) than in ethnically matched control subjects (1/197, p = 0.019).ConclusionsOur findings illustrate the remarkable utility of high-density genotype data for subjects from a limited geographic area in dissecting genetic heterogeneity with respect to population and disease related variation.

Epigenetic differentiation and relationship to adaptive genetic divergence in discrete populations of the violet Viola cazorlensis

*In plants, epigenetic variations based on DNA methylation are often heritable and could influence the course of evolution. Before this hypothesis can be assessed, fundamental questions about epigenetic variation remain to be addressed in a real-world context, including its magnitude, structuring within and among natural populations, and autonomy in relation to the genetic context. *Extent and patterns of cytosine methylation, and the relationship to adaptive genetic divergence between populations, were investigated for wild populations of the southern Spanish violet Viola cazorlensis (Violaceae) using the methylation-sensitive amplified polymorphism (MSAP) technique, a modification of the amplified fragment length polymorphism method (AFLP) based on the differential sensitivity of isoschizomeric restriction enzymes to site-specific cytosine methylation. *The genome of V. cazorlensis plants exhibited extensive levels of methylation, and methylation-based epigenetic variation was structured into distinct between- and within- population components. Epigenetic differentiation of populations was correlated with adaptive genetic divergence revealed by a Bayesian population-genomic analysis of AFLP data. Significant associations existed at the individual genome level between adaptive AFLP loci and the methylation state of methylation-susceptible MSAP loci. *Population-specific, divergent patterns of correlated selection on epigenetic and genetic individual variation could account for the coordinated epigenetic-genetic adaptive population differentiation revealed by this study.

Genetic variations and response to antiepileptic drug

Despite the state-of-the-art medical treatment with antiepileptic drugs (AED), at least one-third of newly diagnosed epilepsy patients may show poor response to AED. Although the biological mechanism of diverse responses to AED is poorly understood, it is likely that multifactorial factors could be responsible for the different responses. Clinical factors such as etiology of epilepsy and pre-treatment seizure frequency can predict the different responses to AED, whereas the role of genetic variations in individuals contributing to variation in response to AED is still conflicting. The inconsistency in proving the genetic roles in different responses to AED may be partly explained by (1) diversity of the definition of response to AED, (2) heterogeneity of the subjects, and (3) lack of multigenetic approach. While we hypothesize that the genetic variations in individuals may independently contribute to different responses to AED, multigenetic interactions including both the genetic variations of pharmacokinetics and pharmacodynamics may take place.

The potential for automated question answering in the context of genomic medicine: an assessment of existing resources and properties of answers

Knowledge gained in studies of genetic disorders is reported in a growing body of biomedical literature containing reports of genetic variation in individuals that map to medical conditions and/or response to therapy. These scientific discoveries need to be translated into practical applications to optimize patient care. Translating research into practice can be facilitated by supplying clinicians with research evidence. We assessed the role of existing tools in extracting answers to translational research questions in the area of genomic medicine. We: evaluate the coverage of translational research terms in the Unified Medical Language Systems (UMLS) Metathesaurus; determine where answers are most often found in full-text articles; and determine common answer patterns. Findings suggest that we will be able to leverage the UMLS in development of natural language processing algorithms for automated extraction of answers to translational research questions from biomedical text in the area of genomic medicine.

Beginning to Unravel the Common, Low-Penetrance Genetic Components Behind Colorectal Cancer by Whole-Genome Association Scans

Beginning to Unravel the Common, Low-Penetrance Genetic Components Behind Colorectal Cancer by Whole-Genome Association Scans

Inbreeding and disease resistance in a social insect: effects of heterozygosity on immunocompetence in the termite Zootermopsis angusticollis

Recent research has shown that low genetic variation in individuals can increase susceptibility to infection and group living may exacerbate pathogen transmission. In the eusocial diploid termites, cycles of outbreeding and inbreeding characterizing basal species can reduce genetic variation within nestmates during the life of a colony, but the relationship of genetic heterogeneity to disease resistance is poorly understood. Here we show that, one generation of inbreeding differentially affects the survivorship of isolated and grouped termites (Zootermopsis angusticollis) depending on the nature of immune challenge and treatment. Inbred and outbred isolated and grouped termites inoculated with a bacterial pathogen, exposed to a low dose of fungal pathogen or challenged with an implanted nylon monofilament had similar levels of immune defence. However, inbred grouped termites exposed to a relatively high concentration of fungal conidia had significantly greater mortality than outbred grouped termites. Inbred termites also had significantly higher cuticular microbial loads, presumably due to less effective grooming by nestmates. Genetic analyses showed that inbreeding significantly reduced heterozygosity and allelic diversity. Decreased heterozygosity thus appeared to increase disease susceptibility by affecting social behaviour or some other group-level process influencing infection control rather than affecting individual immune physiology.

Human leukocyte antigen polymorphisms: variable humoral immune responses to viral vaccines

Antibody formation in response to antigen stimulation remains the basis for measuring an individual’s response and protection for most viral vaccines. A significant proportion of the variation in individual humoral immune response to vaccination appears to be genetic. The collection of genes found on chromosome 6 forming the human leukocyte antigen system provides one of the greatest sources of genetic variation in individuals with respect to their immunological responses. Recent research has demonstrated significant associations between vaccine response and human leukocyte antigen alleles. These associations not only explain why vaccine-induced humoral immune responses vary among individuals and between populations, but these variations may also hold the key to the development of future generations of vaccines.

Reproductive ecology of the freshwater red alga Batrachospermum delicatulum (Batrachospermales, Rhodophyta) in three tropical streams

SUMMARY Batrachospermum delicatulum specimens from three stream segments were analyzed from a tropical region in south‐ eastern Brazil (20°18′– 20°49′S, 49°13′– 49°46′W). Physical and chemical parameters and the spatial placement of thalli were investigated along with the reproductive characteristics of the gametophytic phase. Sequence data of the cox 2‐ 3 spacer region was also utilized to evaluate genetic variation in individuals within and among stream segments. Gametophyte occurred under relatively diverse environmental conditions, whereas thalli abundance was weakly or not correlated to environmental variables within the stream segments. All specimens examined were dioecious. The ratio of male/female plants was relatively low (0.5 to 1.3) and male plants tended to occur as clumps (two or three plants together). High reproductive success was observed, as indicated by the occurrence of 100% fertilized (carposporophytic) female plants. This is similar to previous reports for this and other dioecious species, which is remarkable considering the relatively low proportion of male/female plants. Results support the two hypotheses to explain the high reproductive success in dioecious species. The occurrence of male plants in clumps was evidence for a strict spatial relationship (i.e. male plants located in upstream position of female plants in order to release spermatia, which would be carried by eddies through female plants). In contrast, the occurrence of male and female plants adjacent to each other allowed outcrossing among neighboring plants with intermingled male and female branches, which seemed more applicable to some situations (low turbulence habitats). The cox 2‐ 3 spacer region from the 18 individuals sequenced was 376 bp and the DNA sequence was identical with no base pair substitutions. Likewise, a previous study of another Batrachospermum species showed that the same haplotypes were present in all stream segments from the same drainage basin, even though the stream segments were a considerable distance apart. Short distance dispersal either by small birds or waterway connectivity might explain these findings.

Inherited colour vision deficiencies: From Dalton to molecular genetics

In recent years, great advances have been made in our understanding of the molecular basis of colour vision defects, as well as of the patterns of genetic variation in individuals with normal colour vision. Molecular genetic analyses have explained the diversity of types and degrees of severity in colour vision anomalies, their frequencies, pronounced individual variations in test results, etc. New techniques have even enabled the determination of John Dalton's real colour vision defect, 150 years after his death. Inherited colour vision deficiencies most often result from the mutations of genes that encode cone opsins. Cone opsin genes are linked to chromosomes 7 (the S or "blue" gene) and X (the L or "red" gene and the M or "green" gene). The L and M genes are located on the q arm of the X chromosome in a head-to-tail array, composed of 2 to 6 (typically 3) genes--a single L is followed by one or more M genes. Only the first two genes of the array are expressed and contribute to the colour vision phenotype. The high degree of homology (96%) between the L and M genes predisposes them to unequal recombination, leading to gene deletion or the formation of hybrid genes (comprising portions of both the L and M genes), explaining the majority of the common red-green colour vision deficiencies. The severity of any deficiency is influenced by the difference in spectral sensitivity between the opsins encoded by the first two genes of the array. A rare defect, S monochromacy, is caused either by the deletion of the regulatory region of the array or by mutations that inactivate the L and M genes. Most recent research concerns the molecular basis of complete achromatopsia, a rare disorder that involves the complete loss of all cone function. This is not caused by mutations in opsin genes, but in other genes that encode cone-specific proteins, e.g. channel proteins and transducin.

Organization of genetic variation in individuals of arbuscular mycorrhizal fungi.

Arbuscular mycorrhizal (AM) fungi (Glomeromycota) are thought to be the oldest group of asexual multicellular organisms. They colonize the roots of most land plants, where they facilitate mineral uptake from the soil in exchange for plant-assimilated carbon. Cells of AM fungi contain hundreds of nuclei. Unusual polymorphism of ribosomal DNA observed in individual spores of AM fungi inspired a hypothesis that heterokaryosis--that is, the coexistence of many dissimilar nuclei in cells--occurs throughout the AM fungal life history. Here we report a genetic approach to test the hypothesis of heterokaryosis in AM fungi. Our study of the transmission of polymorphic genetic markers in natural isolates of Glomus etunicatum, coupled with direct amplification of rDNA from microdissected nuclei by polymerase chain reaction, supports the alternative hypothesis of homokaryosis, in which nuclei populating AM fungal individuals are genetically uniform. Intrasporal rDNA polymorphism contained in each nucleus signals a relaxation of concerted evolution, a recombination-driven process that is responsible for homogenizing rDNA repeats. Polyploid organization of glomeromycotan genomes could accommodate intranuclear rDNA polymorphism and buffer these apparently asexual organisms against the effects of accumulating mutations.

Pyrosequencing AB.

Pyrosequencing AB develops, manufactures and sells research products that enable life sciences researchers to add utility to the massive amount of genomic information available, with the goal of improving human health. As more genomic information becomes available, research efforts are shifting toward implementing practical applications of this genetic information. Pyrosequencing has developed comprehensive genotyping solutions that provide researchers with rapid and accurate technology for assessing individual genetic variation. Genotypic information is anticipated to play an important role in the emerging field of applied genomics and as a result, significant opportunities may be derived from analyzing genetic variation in individuals and across populations. These data will be key to understanding how genetic differences may be used to determine an individual's predisposition to disease, detect and monitor disease progression, and select the most appropriate therapy and dose and to avoid adverse events. This is pharmacogenomics, an area in which Pyrosequencing technology is being applied and will no doubt continue to enable researchers to ultimately improve human health.

Comparison of Three Anadromous Rainbow Trout (Oncorhynchus mykiss) Populations Using DNA Fingerprinting and Mixed DNA Samples

DNA fingerprinting of 15 individuals from three populations of steelhead trout (Oncorhynchus mykiss) was compared with DNA fingerprinting of mixtures of 15 individuals from the same populations. Three oligonucleotide fingerprinting probes were used to detect genetic variation in individuals and mixtures. Scanning image analysis coupled with customized software was used to assign band identity. The degree of band sharing of individuals within and between populations was then determined. Although band-sharing estimates were similar within and between populations, bands found predominately in one population were identified. DNA fingerprinting of mixtures detected bands found only in the majority of the individuals in the mixture, allowing differences between populations to be more readily identified.

A Technique for the Investigation of Genetic Variation in Filamentous Gametophytes of Trichomanes (Hymenophyllaceae)

The gametophyte of most ferns is a small, thalloid, heart-shaped, short-lived structure. In some species of Trichomanes (Hymenophyllaceae), the gametophytes are filamentous and consist of a weft of inter-twining branched filaments. Such gametophytes are perennial and can reproduce asexually by means of gemmae. Independent colonies of gametophytes have been reported from North America (Farrar, 1967, 1992), Britain (Rumsey et al., 1990) and Central Europe (Vogel et al., 1993, Jerome and Rasbach, 1994). These colonies of filamentous gametophytes can exist as very small patches (<10mm2) or can cover huge areas (20m2). Such a colony, whether large or small, may have arisen from a single spore, from several spores or, asexually, by establishment from one or more gemmae. A colony of filaments may represent either a single individual, or consist of many individuals interwoven. This presents a problem for the study of genetic variation of individuals in these independent colonies when using conventional methods, eg. allozyme analysis (Hamrick, 1990). Since the maximum length of a continuous individual filament that can be isolated from a colony is only about 60-100 cells, conventional allozyme analysis and DNA extraction methods are inapplicable. To extract DNA from a single gametophyte filament, ultrasonication was used directly to provide a crude DNA extract for DNA amplification. To establish the extraction and amplification techniques, gametophytes of four taxa of Trichomanes were compared using a large single-copy region of chloroplast DNA. Taberlet et al. (1991) have shown that the use of their universal primers for DNA amplification of the chloroplast intergenic spacers between the trnT and trnF genes is informative for a wide range of plant species. MATERIALS & METHODS