Genetically-regulated transcriptomics & copy number variation of proctitis points to altered mitochondrial and DNA repair mechanisms in individuals of European ancestry.

BACKGROUND: Proctitis is an inflammation of the rectum and may be induced by radiation treatment for cancer. The genetic heritability of developing radiotoxicity and prior role of genetic variants as being associated with side-effects of radiotherapy necessitates further investigation for underlying molecular mechanisms. In this study, we investigated gene expression regulated by genetic variants, and copy number variation in prostate cancer survivors with radiotoxicity. METHODS: We investigated proctitis as a radiotoxic endpoint in prostate cancer patients who received radiotherapy (n = 222). We analyzed the copy number variation and genetically regulated gene expression profiles of whole-blood and prostate tissue associated with proctitis. The SNP and copy number data were genotyped on Affymetrix® Genome-wide Human SNP Array 6.0. Following QC measures, the genotypes were used to obtain gene expression by leveraging GTEx, a reference dataset for gene expression association based on genotype and RNA-seq information for prostate (n = 132) and whole-blood tissue (n = 369). RESULTS: In prostate tissue, 62 genes were significantly associated with proctitis, and 98 genes in whole-blood tissue. Six genes - CABLES2, ATP6AP1L, IFIT5, ATRIP, TELO2, and PARD6G were common to both tissues. The copy number analysis identified seven regions associated with proctitis, one of which (ALG1L2) was also associated with proctitis based on transcriptomic profiles in the whole-blood tissue. The genes identified via transcriptomics and copy number variation association were further investigated for enriched pathways and gene ontology. Some of the enriched processes were DNA repair, mitochondrial apoptosis regulation, cell-to-cell signaling interaction processes for renal and urological system, and organismal injury. CONCLUSIONS: We report gene expression changes based on genetic polymorphisms. Integrating gene-network information identified these genes to relate to canonical DNA repair genes and processes. This investigation highlights genes involved in DNA repair processes and mitochondrial malfunction possibly via inflammation. Therefore, it is suggested that larger studies will provide more power to infer the extent of underlying genetic contribution for an individual's susceptibility to developing radiotoxicity.

Paternity calculations in a di-spermy case.

In a criminal paternity case, which involved analysis of the product of conception, a rare circumstance was observed. The product of conception was triploidy, apparently due to an egg fertilized by two sperm. Since there is little guidance on how to calculate the probability of the DNA evidence given some basic hypotheses, the formulae were derived and are presented herein. These approaches could provide guidance for similar situations if they arise.

Resolving the etiology of atopic disorders by using genetic analysis of racial ancestry.

Atopic dermatitis (AD), food allergy, allergic rhinitis, and asthma are common atopic disorders of complex etiology. The frequently observed atopic march from early AD to asthma, allergic rhinitis, or both later in life and the extensive comorbidity of atopic disorders suggest common causal mechanisms in addition to distinct ones. Indeed, both disease-specific and shared genomic regions exist for atopic disorders. Their prevalence also varies among races; for example, AD and asthma have a higher prevalence in African Americans when compared with European Americans. Whether this disparity stems from true genetic or race-specific environmental risk factors or both is unknown. Thus far, the majority of the genetic studies on atopic diseases have used populations of European ancestry, limiting their generalizability. Large-cohort initiatives and new analytic methods, such as admixture mapping, are currently being used to address this knowledge gap. Here we discuss the unique and shared genetic risk factors for atopic disorders in the context of ancestry variations and the promise of high-throughput "-omics"-based systems biology approach in providing greater insight to deconstruct their genetic and nongenetic etiologies. Future research will also focus on deep phenotyping and genotyping of diverse racial ancestry, gene-environment, and gene-gene interactions.

Selection pressure on human STR loci and its relevance in repeat expansion disease.

Short Tandem Repeats (STRs) comprise repeats of one to several base pairs. Because of the high mutability due to strand slippage during DNA synthesis, rapid evolutionary change in the number of repeating units directly shapes the range of repeat-number variation according to selection pressure. However, the remaining questions include: Why are STRs causing repeat expansion diseases maintained in the human population; and why are these limited to neurodegenerative diseases? By evaluating the genome-wide selection pressure on STRs using the database we constructed, we identified two different patterns of relationship in repeat-number polymorphisms between DNA and amino-acid sequences, although both patterns are evolutionary consequences of avoiding the formation of harmful long STRs. First, a mixture of degenerate codons is represented in poly-proline (poly-P) repeats. Second, long poly-glutamine (poly-Q) repeats are favored at the protein level; however, at the DNA level, STRs encoding long poly-Qs are frequently divided by synonymous SNPs. Furthermore, significant enrichments of apoptosis and neurodevelopment were biological processes found specifically in genes encoding poly-Qs with repeat polymorphism. This suggests the existence of a specific molecular function for polymorphic and/or long poly-Q stretches. Given that the poly-Qs causing expansion diseases were longer than other poly-Qs, even in healthy subjects, our results indicate that the evolutionary benefits of long and/or polymorphic poly-Q stretches outweigh the risks of long CAG repeats predisposing to pathological hyper-expansions. Molecular pathways in neurodevelopment requiring long and polymorphic poly-Q stretches may provide a clue to understanding why poly-Q expansion diseases are limited to neurodegenerative diseases.

Effects of the Ion PGM™ Hi-Q™ sequencing chemistry on sequence data quality.

Massively parallel sequencing (MPS) offers substantial improvements over current forensic DNA typing methodologies such as increased resolution, scalability, and throughput. The Ion PGM™ is a promising MPS platform for analysis of forensic biological evidence. The system employs a sequencing-by-synthesis chemistry on a semiconductor chip that measures a pH change due to the release of hydrogen ions as nucleotides are incorporated into the growing DNA strands. However, implementation of MPS into forensic laboratories requires a robust chemistry. Ion Torrent's Hi-Q™ Sequencing Chemistry was evaluated to determine if it could improve on the quality of the generated sequence data in association with selected genetic marker targets. The whole mitochondrial genome and the HID-Ion STR 10-plex panel were sequenced on the Ion PGM™ system with the Ion PGM™ Sequencing 400 Kit and the Ion PGM™ Hi-Q™ Sequencing Kit. Concordance, coverage, strand balance, noise, and deletion ratios were assessed in evaluating the performance of the Ion PGM™ Hi-Q™ Sequencing Kit. The results indicate that reliable, accurate data are generated and that sequencing through homopolymeric regions can be improved with the use of Ion Torrent's Hi-Q™ Sequencing Chemistry. Overall, the quality of the generated sequencing data supports the potential for use of the Ion PGM™ in forensic genetic laboratories.

Empirical testing of a 23-AIMs panel of SNPs for ancestry evaluations in four major US populations.

Ancestry informative markers (AIMs) can be used to determine population affiliation of the donors of forensic samples. In order to examine ancestry evaluations of the four major populations in the USA, 23 highly informative AIMs were identified from the International HapMap project. However, the efficacy of these 23 AIMs could not be fully evaluated in silico. In this study, these 23 SNPs were multiplexed to test their actual performance in ancestry evaluations. Genotype data were obtained from 189 individuals collected from four American populations. One SNP (rs12149261) on chromosome 16 was removed from this panel because it was duplicated on chromosome 1. The resultant 22-AIMs panel was able to empirically resolve the four major populations as in the in silico study. Eight individuals were assigned to a different group than indicated on their samples. The assignments of the 22 AIMs for these samples were consistent with AIMs results from the ForenSeq(TM) panel. No departures from Hardy-Weinberg equilibrium (HWE) and linkage disequilibrium (LD) were detected for all 22 SNPs in four US populations (after removing the eight problematic samples). The principal component analysis (PCA) results indicated that 181 individuals from these populations were assigned to the expected groups. These 22 SNPs can contribute to the candidate AIMs pool for potential forensic identification purposes in major US populations.

Selection of highly informative SNP markers for population affiliation of major US populations.

Ancestry informative markers (AIMs) can be used to detect and adjust for population stratification and predict the ancestry of the source of an evidence sample. Autosomal single nucleotide polymorphisms (SNPs) are the best candidates for AIMs. It is essential to identify the most informative AIM SNPs across relevant populations. Several informativeness measures for ancestry estimation have been used for AIMs selection: absolute allele frequency differences (δ), F statistics (F ST), and informativeness for assignment measure (In). However, their efficacy has not been compared objectively, particularly for determining affiliations of major US populations. In this study, these three measures were directly compared for AIMs selection among four major US populations, i.e., African American, Caucasian, East Asian, and Hispanic American. The results showed that the F ST panel performed slightly better for population resolution based on principal component analysis (PCA) clustering than did the δ panel and both performed better than the In panel. Therefore, the 23 AIMs selected by the F ST measure were used to characterize the four major American populations. Genotype data of nine sample populations were used to evaluate the efficiency of the 23-AIMs panel. The results indicated that individuals could be correctly assigned to the major population categories. Our AIMs panel could contribute to the candidate pool of AIMs for potential forensic identification purposes.

Genetic and linguistic correlation of the Kra-Dai-speaking groups in Thailand.

The Kra-Dai linguistic family includes Thai and Lao as well as a great number of languages spoken by ethnic minorities in Southeast Asia. In Thailand, a dozen of other Kra-Dai languages are spoken in addition to Thai, the national language. The genetic structure of the Kra-Dai-speaking populations in Thailand has been studied extensively using uniparentally inherited markers. To extend this line of genetic investigation, this study used 15 autosomal microsatellites of 500 individuals from 11 populations, belonging to nine Kra-Dai ethnicities, namely, the Kaleung, Phu Thai, Saek, Nyo, Lao Isan, Yuan, Black Tai, Phuan and Lue. These ethnolinguistic groups are dispersed in three different geographic regions of Thailand, that is, Northern, Northeastern and Central. The results show a very low average of pairwised F(st) (0.0099), as well as no population substructure based on STRUCTURE analysis, indicating genetic homogeneity within the Kra-Dai-speaking group, possibly owing to shared linguistic ancestry. The Mantel test, an analysis of molecular variance, and the approximate Bayesian computation procedure employed to evaluate potential factors for driving genetic diversity revealed that language is the predominant factor affecting genetic variations, whereas geography is not. The result of distance-based clustering analyses and spatial analysis of molecular variance revealed genetic distinctions of some populations, reflecting the effects of genetic drift and gene flow on allele frequency within populations, in concordance with the result of R-matrix regression. The genetic and linguistic affiliations of the contemporary Kra-Dai-speaking groups are consistent with each other despite certain deviation due to various evolutionary factors that may have occurred during their migrations and resettlements.

Role of genomics in eliminating health disparities.

The Texas Center for Health Disparities, a National Institute on Minority Health and Health Disparities Center of Excellence, presents an annual conference to discuss prevention, awareness education, and ongoing research about health disparities both in Texas and among the national population. The 2014 Annual Texas Conference on Health Disparities brought together experts in research, patient care, and community outreach on the "Role of Genomics in Eliminating Health Disparities." Rapid advances in genomics and pharmacogenomics are leading the field of medicine to use genetics and genetic risk to build personalized or individualized medicine strategies. We are at a critical juncture of ensuring such rapid advances benefit diverse populations. Relatively few forums have been organized around the theme of the role of genomics in eliminating health disparities. The conference consisted of three sessions addressing "Gene-Environment Interactions and Health Disparities," "Personalized Medicine and Elimination of Health Disparities," and "Ethics and Public Policy in the Genomic Era." This article summarizes the basic science, clinical correlates, and public health data presented by the speakers.

Rank-based genome-wide analysis reveals the association of ryanodine receptor-2 gene variants with childhood asthma among human populations.

BACKGROUND: The standard approach to determine unique or shared genetic factors across populations is to identify risk alleles in one population and investigate replication in others. However, since populations differ in DNA sequence information, allele frequencies, effect sizes, and linkage disequilibrium patterns, SNP association using a uniform stringent threshold on p values may not be reproducible across populations. Here, we developed rank-based methods to investigate shared or population-specific loci and pathways for childhood asthma across individuals of diverse ancestry. We performed genome-wide association studies on 859,790 SNPs genotyped in 527 affected offspring trios of European, African, and Hispanic ancestry using publically available asthma database in the Genotypes and Phenotypes database. RESULTS: Rank-based analyses showed that there are shared genetic factors for asthma across populations, more at the gene and pathway levels than at the SNP level. Although the top 1,000 SNPs were not shared, 11 genes (RYR2, PDE4D, CSMD1, CDH13, ROBO2, RBFOX1, PTPRD, NPAS3, PDE1C, SEMA5A, and CTNNA2) mapped by these SNPs were shared across populations. Ryanodine receptor 2 (RYR2, a statin response-related gene) showed the strongest association in European (p value=2.55×10(-7)) and was replicated in African (2.57×10(-4)) and Hispanic (1.18 × 10(-3)) Americans. Imputation analyses based on the 1000 Genomes Project uncovered additional RYR2 variants associated with asthma. Network and functional ontology analyses revealed that RYR2 is an integral part of dermatological or allergic disorder biological networks, specifically in the functional classes involving inflammatory, eosinophilic, and respiratory diseases. CONCLUSION: Our rank-based genome-wide analysis revealed for the first time an association of RYR2 variants with asthma and replicated previously discovered PDE4D asthma gene across human populations. The replication of top-ranked asthma genes across populations suggests that such loci are less likely to be false positives and could indicate true associations. Variants that are associated with asthma across populations could be used to identify individuals who are at high risk for asthma regardless of genetic ancestry.

A prioritization analysis of disease association by data-mining of functional annotation of human genes.

Complex diseases result from contributions of multiple genes that act in concert through pathways. Here we present a method to prioritize novel candidates of disease-susceptibility genes depending on the biological similarities to the known disease-related genes. The extent of disease-susceptibility of a gene is prioritized by analyzing seven features of human genes captured in H-InvDB. Taking rheumatoid arthritis (RA) and prostate cancer (PC) as two examples, we evaluated the efficiency of our method. Highly scored genes obtained included TNFSF12 and OSM as candidate disease genes for RA and PC, respectively. Subsequent characterization of these genes based upon an extensive literature survey reinforced the validity of these highly scored genes as possible disease-susceptibility genes. Our approach, Prioritization ANalysis of Disease Association (PANDA), is an efficient and cost-effective method to narrow down a large set of genes into smaller subsets that are most likely to be involved in the disease pathogenesis.

Uncertainties in estimating health risks associated with exposure to ionising radiation.

The information for the present discussion on the uncertainties associated with estimation of radiation risks and probability of disease causation was assembled for the recently published NCRP Report No. 171 on this topic. This memorandum provides a timely overview of the topic, given that quantitative uncertainty analysis is the state of the art in health risk assessment and given its potential importance to developments in radiation protection. Over the past decade the increasing volume of epidemiology data and the supporting radiobiology findings have aided in the reduction of uncertainty in the risk estimates derived. However, it is equally apparent that there remain significant uncertainties related to dose assessment, low dose and low dose-rate extrapolation approaches (e.g. the selection of an appropriate dose and dose-rate effectiveness factor), the biological effectiveness where considerations of the health effects of high-LET and lower-energy low-LET radiations are required and the transfer of risks from a population for which health effects data are available to one for which such data are not available. The impact of radiation on human health has focused in recent years on cancer, although there has been a decided increase in the data for noncancer effects together with more reliable estimates of the risk following radiation exposure, even at relatively low doses (notably for cataracts and cardiovascular disease). New approaches for the estimation of hereditary risk have been developed with the use of human data whenever feasible, although the current estimates of heritable radiation effects still are based on mouse data because of an absence of effects in human studies. Uncertainties associated with estimation of these different types of health effects are discussed in a qualitative and semi-quantitative manner as appropriate. The way forward would seem to require additional epidemiological studies, especially studies of low dose and low dose-rate occupational and perhaps environmental exposures and for exposures to x rays and high-LET radiations used in medicine. The development of models for more reliably combining the epidemiology data with experimental laboratory animal and cellular data can enhance the overall risk assessment approach by providing biologically refined data to strengthen the estimation of effects at low doses as opposed to the sole use of mathematical models of epidemiological data that are primarily driven by medium/high doses. NASA's approach to radiation protection for astronauts, although a unique occupational group, indicates the possible applicability of estimates of risk and their uncertainty in a broader context for developing recommendations on: (1) dose limits for occupational exposure and exposure of members of the public; (2) criteria to limit exposures of workers and members of the public to radon and its short-lived decay products; and (3) the dosimetric quantity (effective dose) used in radiation protection.

Genome-wide association of serum uric acid concentration: replication of sequence variants in an island population of the Adriatic coast of Croatia.

A genome-wide association study of serum uric acid (SUA) laevels was performed in a relatively isolated population of European descent from an island of the Adriatic coast of Croatia. The study sample included 532 unrelated and 768 related individuals from 235 pedigrees. Inflation due to relatedness was controlled by using genomic control. Genetic association was assessed with 2,241,249 single nucleotide polymorphisms (SNPs) in 1300 samples after adjusting for age and gender. Our study replicated four previously reported SUA loci (SLC2A9, ABCG2, RREB1, and SLC22A12). The strongest association was found with a SNP in SLC2A9 (rs13129697, P=2.33×10(-19)), which exhibited significant gender-specific effects, 35.76 µmol/L (P=2.11×10(-19)) in females and 19.58 µmol/L (P=5.40×10(-5)) in males. Within this region of high linkage disequilibrium, we also detected a strong association with a nonsynonymous SNP, rs16890979 (P=2.24×10(-17)), a putative causal variant for SUA variation. In addition, we identified several novel loci suggestive of association with uric acid levels (SEMA5A, TMEM18, SLC28A2, and ODZ2), although the P-values (P<5×10(-6)) did not reach the threshold of genome-wide significance. Together, these findings provide further confirmation of previously reported uric-acid-related genetic variants and highlight suggestive new loci for additional investigation.

Pedigree likelihood ratio for lineage markers.

Lineage-based haplotype markers (e.g., Y chromosome STRs and mitochondrial DNA sequences) are important adjunct tools to the autosomal markers for kinship analysis and for specialized kinship applications such as database searching. Traditionally, the prosecution or kinship hypothesis considers the haplotypes in the same lineage and the probability of genotype data given the lineage hypothesis is simply set at 1 if the number of mismatched loci or nucleotides between the questioned person and the references is less than a predefined threshold. In this study, a kinship hypothesis based on a fixed relationship of the questioned person in the reference family is introduced. A graphical model is proposed to calculate the probability of the genotype data given the kinship hypothesis, which is the product of haplotype frequency of the founder in the pedigree and the transmission probability from the founder to all descendants. Proper mutation models are suggested for Y chromosome STRs and mitochondrial DNA sequence variants (i.e., SNPs) to calculate the transmission probability. The methods to infer the genotypes of the untyped individuals in the pedigree and the computational complexity of handling these untyped individuals are also addressed. Lastly, numerical examples of the applications are given to demonstrate the kinship hypothesis and the algorithms.

Population genetic analyses of the NGM STR loci.

The AmpFlSTR® NGM™ PCR Amplification Kit enables amplification of 15 autosomal short tandem repeat (STR) loci. The loci are the ten STRs in the SGM Plus® Kit plus the EDNAP and ENSFI recommended STRs D10S1248, D22S1045, D2S441, D1S1656, and D12S391. Allele frequency and other forensically relevant statistics data were generated for the NGM loci in three US population groups (African Americans, Caucasians, and Hispanics). The analyses support that the NGM multiplex is one of the most informative STR multiplex kits available to the forensic science community. At the population level, there are no more detectable departures from expectations of the independence of alleles within as well as between loci than would be expected due to chance, even for the two syntenic loci vWA and D12S391; however, linkage analysis in three large pedigree families shows close linkage between these two loci with a recombination fraction of 0.108. Therefore, in contrast to the practices in calculating the rarity of a DNA profile, for kinship analyses independence between the loci, vWA and D12S391 cannot be assumed.

Y-STR loci diversity in native Alaskan populations.

Y chromosome short tandem repeat (Y-STR) loci are important genetic markers for forensic biological evidence analyses. However, paternal inheritance, reduced effective population size, and lack of independence between loci can reduce Y-STR diversity and may yield greater population substructure effects on a locus-by-locus basis compared with the autosomal STR loci. Population studies are necessary to assess the genetic variation of forensically relevant markers so that proper inferences can be made about the rarity of DNA profiles. This study examined 16 Y-STRs in three sampled populations of Native Americans from Alaska: Inupiat, Yupik, and Athabaskan. Population genetic and statistical issues addressed were: (1) the degree of diversity at locus and haplotype levels, (2) determination of the loci that contribute more so to haplotype diversity, and (3) the effects of population substructure on forensic statistical calculations of the rarity of a Y-STR profile. All three population samples were highly polymorphic at the haplotype level for the 16 Y-STR markers; however, the Native Americans demonstrated reduced genetic diversity compared with major US populations. The degree of substructure indicated that the three populations were related and admixed in terms of paternal lineage. The examination of more polymorphic loci may be needed to increase the power of discrimination of Y-STR systems in these populations.

Analysis of human CYP1A1 and CYP1A2 genes and their shared bidirectional promoter in eight world populations.

The human CYP1A1_CYP1A2 locus comprises the CYP1A1 (5,988 bp) and CYP1A2 (7,759 bp) transcribed regions, oriented head-to-head, sharing a bidirectional promoter of 23,306 bp. The older CYP1A1 gene appears more conserved and responsible for critical life function(s), whereas the younger CYP1A2 gene might have evolved more rapidly due to environmental (dietary) pressures. A population genetics study might confirm this premise. We combined 60 CYP1A1_CYP1A2 SNPs found in the present study (eight New Guinea Highlanders, eight Samoans, four Dogrib, four Teribe, four Pehuenche, and one Caucasian) with those found in a previous study (six West Africans, four Han Chinese, six Germans, four Samoans, and four Dogrib), yielding a total of 106 SNPs in 106 chromosomes. Resequencing of Oceanians plus Amerindians in the present study yielded 21 New World SNPs ( approximately 20%), of which 17 are not previously reported in any SNP database. Various tests revealed selective pressures for both genes and both haploblocks; unfortunately, differences in rates of evolution between the two genes were undetectable. Fay & Wu's H test revealed a "hitchhiking event" centered around four SNPs in the CYP1A1 3'-UTR; a study in silico identified different microRNA-binding patterns in the hitchhiked region, when the mutations were present compared with the mutations absent.

Haplotype block: a new type of forensic DNA markers.

Forensic DNA analysis is currently performed using highly discriminating short tandem repeat (STR) markers. SNPs are being investigated as adjunct tools for human identity testing because of their abundance in the human genome, utility for genotyping degraded DNA samples, and amenability to automation. While SNPs can provide an alternative approach, on a per locus basis they have a lower power of discrimination (PD) than STRs. With the discovery of block structures in the human genome, a novel set of SNP markers are available for further exploration of forensic utility. Several neighboring, tightly linked SNPs are inherited together and form a haplotype block, which as a haploblock has a higher discrimination power than the individual SNPs within the block. Candidate haplotype blocks were selected from three major populations (Caucasian, East Asian, and African) using the following parameters: maximum match probability reduction = 0.85, linkage disequilibrium (LD) r(2) ≥ 0.7, maximum F(st) = 0.06, minimum number of SNPs = 3, minimum heterozygosity = 0.2, and minimum number of haplotypes = 3. From the HapMap Phase II data, 253 haploblocks were identified on the 22 autosomal chromosomes. After removing haploblocks deviating from the Hardy-Weinberg equilibrium (HWE) or in LD with other haploblocks, 24 haploblocks remained as candidates for forensic consideration. The cumulative PD of these blocks can reach 10(-12) in the populations studied. The data support within and between haplotype independence even when they are syntenic. We propose guidelines for evidence interpretation that address the application of haplotype blocks for transfer evidence, mixture, and kinship analyses.

Genetic evidence supports demic diffusion of Han culture.

The spread of culture and language in human populations is explained by two alternative models: the demic diffusion model, which involves mass movement of people; and the cultural diffusion model, which refers to cultural impact between populations and involves limited genetic exchange between them. The mechanism of the peopling of Europe has long been debated, a key issue being whether the diffusion of agriculture and language from the Near East was concomitant with a large movement of farmers. Here we show, by systematically analysing Y-chromosome and mitochondrial DNA variation in Han populations, that the pattern of the southward expansion of Han culture is consistent with the demic diffusion model, and that males played a larger role than females in this expansion. The Han people, who all share the same culture and language, exceed 1.16 billion (2000 census), and are by far the largest ethnic group in the world. The expansion process of Han culture is thus of great interest to researchers in many fields.