Independent polled mutations leading to complex gene expression differences in cattle

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.

Clinical importance of risk variants in the dihydropyrimidine dehydrogenase gene for the prediction of early-onset fluoropyrimidine toxicity.

We investigated the clinical relevance of dihydropyrimidine dehydrogenase gene (DPYD) variants to predict severe early-onset fluoropyrimidine (FP) toxicity, in particular of a recently discovered haplotype hapB3 and a linked deep intronic splice site mutation c.1129-5923C>G. Selected regions of DPYD were sequenced in prospectively collected germline DNA of 500 patients receiving FP-based chemotherapy. Associations of DPYD variants and haplotypes with hematologic, gastrointestinal, infectious, and dermatologic toxicity in therapy cycles 1-2 and resulting FP-dose interventions (dose reduction, therapy delay or cessation) were analyzed accounting for clinical and demographic covariates. Fifteen additional cases with toxicity-related therapy delay or cessation were retrospectively examined for risk variants. The association of c.1129-5923C>G/hapB3 (4.6% carrier frequency) with severe toxicity was replicated in an independent prospective cohort. Overall, c.1129-5923G/hapB3 carriers showed a relative risk of 3.74 (RR, 95% CI = 2.30-6.09, p = 2 × 10(-5)) for severe toxicity (grades 3-5). Of 31 risk variant carriers (c.1129-5923C>G/hapB3, c.1679T>G, c.1905+1G>A or c.2846A>T), 11 (all with c.1129-5923C>G/hapB3) experienced severe toxicity (15% of 72 cases, RR = 2.73, 95% CI = 1.61-4.63, p = 5 × 10(-6)), and 16 carriers (55%) required FP-dose interventions. Seven of the 15 (47%) retrospective cases carried a risk variant. The c.1129-5923C>G/hapB3 variant is a major contributor to severe early-onset FP toxicity in Caucasian patients. This variant may substantially improve the identification of patients at risk of FP toxicity compared to established DPYD risk variants (c.1905+1G>A, c.1679T>G and c.2846A>T). Pre-therapeutic DPYD testing may prevent 20-30% of life-threatening or lethal episodes of FP toxicity in Caucasian patients.

Identification of genetic variation influencing apolipoprotein E levels: Follow-up of genome-wide linkage scans in the GENOA study

Apolipoprotein E (ApoE) plays a major role in the metabolism of high density and low density lipoproteins (HDL and LDL). Its common protein isoforms (E2, E3, E4) are risk factors for coronary artery disease (CAD) and explain between 16 to 23% of the inter-individual variation in plasma apoE levels. Linkage analysis has been completed for plasma apoE levels in the GENOA study (Genetic Epidemiology Network of Atherosclerosis). After stratification of the population by lipoprotein levels and body mass index (BMI) to create more homogeneity with regard to biological context for apoE levels, Hispanic families showed significant linkage on chromosome 17q for two strata (LOD=2.93 at 104 cM for a low cholesterol group, LOD=3.04 at 111 cM for a low cholesterol, high HDLC group). Replication of 17q linkage was observed for apoB and apoE levels in the unstratified Hispanic and African-American populations, and for apoE levels in African-American families. Replication of this 17q linkage in different populations and strata provides strong support for the presence of gene(s) in this region with significant roles in the determination of inter-individual variation in plasma apoE levels. Through a positional and functional candidate gene approach, ten genes were identified in the 17q linked region, and 62 polymorphisms in these genes were genotyped in the GENOA families. Association analysis was performed with FBAT, GEE, and variance-component based tests followed by conditional linkage analysis. Association studies with partial coverage of TagSNPs in the gene coding for apolipoprotein H (APOH) were performed, and significant results were found for 2 SNPs (APOH_20951 and APOH_05407) in the Hispanic low cholesterol strata accounting for 3.49% of the inter-individual variation in plasma apoE levels. Among the other candidate genes, we identified a haplotype block in the ACE1 gene that contains two major haplotypes associated with apoE levels as well as total cholesterol, apoB and LDLC levels in the unstratified Hispanic population. Identifying genes responsible for the remaining 60% of inter-individual variation in plasma apoE level, will yield new insights into the understanding of genetic interactions involved in the lipid metabolism, and a more precise understanding of the risk factors leading to CAD. ^

Characterization of glucocorticoid receptor (NR3C1) gene polymorphisms and a family-based association study with hypertension

Hypertension is usually defined as having values of systolic blood pressure ≥140 mmHg, diastolic blood pressure ≥90 mmHg. Hypertension is one of the main adverse effects of glucocorticoid on the cardiovascular system. Glucocorticoids are essential hormones, secreted from adrenal glands in circadian fashion. Glucocorticoid's effect on blood pressure is conveyed by the glucocorticoid receptor (NR3C1), an omnipresent nuclear transcription factor. Although polymorphisms in this gene have long been implicated to be a causal factor for cardiovascular diseases such as hypertension, no study has yet thoroughly interrogated the gene's polymorphisms for their effect on blood pressure levels. Therefore, I have first resequenced ∼30 kb of the gene, encompassing all exons, promoter regions, 5'/3' UTRs as well as at least 1.5 kb of the gene's flanking regions from 114 chromosome 5 monosomic cell lines, comprised of three major American ethnic groups—European American, African American and Mexican American. I observed 115 polymorphisms and 14 common molecularly phased haplotypes. A subset of markers was chosen for genotyping study populations of GENOA (Genetic Epidemiology Network of Atherosclerosis; 1022 non-Hispanic whites, 1228 African Americans and 954 Mexican Americans). Since these study populations include sibships, the family-based association test was performed on 4 blood pressure-related quantitative variables—pulse, systolic blood pressure, diastolic blood pressure and mean arterial pressure. Using these analyses, multiple correlated SNPs are significantly protective against high systolic blood pressure in non-Hispanic whites, which includes rsb198, a SNP formerly associated with beneficial body compositions. Haplotype association analysis also supports this finding and all p-values remained significant after permutation tests. I therefore conclude that multiple correlated SNPs on the gene may confer protection against high blood pressure in non-Hispanic whites. ^

Genetic variations in the PI3K-AKT-mTOR pathway and bladder cancer susceptibility and clinical outcomes

Bladder cancer is the fourth most common cancer in men in the United States. There is compelling evidence supporting that genetic variations contribute to the risk and outcomes of bladder cancer. The PI3K-AKT-mTOR pathway is a major cellular pathway involved in proliferation, invasion, inflammation, tumorigenesis, and drug response. Somatic aberrations of PI3K-AKT-mTOR pathway are frequent events in several cancers including bladder cancer; however, no studies have investigated the role of germline genetic variations in this pathway in bladder cancer. In this project, we used a large case control study to evaluate the associations of a comprehensive catalogue of SNPs in this pathway with bladder cancer risk and outcomes. Three SNPs in RAPTOR were significantly associated with susceptibility: rs11653499 (OR: 1.79, 95%CI: 1.24–2.60), rs7211818 (OR: 2.13, 95%CI: 1.35–3.36), and rs7212142 (OR: 1.57, 95%CI: 1.19–2.07). Two haplotypes constructed from these 3 SNPs were also associated with bladder cancer risk. In combined analysis, a significant trend was observed for increased risk with an increase in the number of unfavorable genotypes (P for trend<0.001). Classification and regression tree analysis identified potential gene-environment interactions between RPS6KA5 rs11653499 and smoking. In superficial bladder cancer, we found that PTEN rs1234219 and rs11202600, TSC1 rs7040593, RAPTOR rs901065, and PIK3R1 rs251404 were significantly associated with recurrence in patients receiving BCG. In muscle invasive and metastatic bladder cancer, AKT2 rs3730050, PIK3R1 rs10515074, and RAPTOR rs9906827 were associated with survival. Survival tree analysis revealed potential gene-gene interactions: patients carrying the unfavorable genotypes of PTEN rs1234219 and TSC1 rs704059 exhibited a 5.24-fold (95% CI: 2.44–11.24) increased risk of recurrence. In combined analysis, with the increasing number of unfavorable genotypes, there was a significant trend of higher risk of recurrence and death (P for trend<0.001) in Cox proportional hazard regression analysis, and shorter event (recurrence and death) free survival in Kaplan-Meier estimates (P log rank<0.001). This study strongly suggests that genetic variations in PI3K-AKT-mTOR pathway play an important role in bladder cancer development. The identified SNPs, if validated in further studies, may become valuable biomarkers in assessing an individual's cancer risk, predicting prognosis and treatment response, and facilitating physicians to make individualized treatment decisions. ^

Genetic variability and the hypothalamic control of energy balance

Background. Obesity is a major health problem throughout the industrialized world. Despite numerous attempts to curtail the rapid growth of obesity, its incidence continues to rise. Therefore, it is crucial to better understand the etiology of obesity beyond the concept of energy balance.^ Aims. The first aim of this study was to first investigate the relationship between eating behaviors and body size. The second goal was to identify genetic variation associated with eating behaviors. Thirdly, this study aimed to examine the joint relationships between eating behavior, body size and genetic variation.^ Methods. This study utilized baseline data ascertained in young adults from the Training Interventions and Genetics of Exercise (TIGER) Study. Variables assessed included eating behavior (Emotional Eating Scale, Eating Attitudes Test-26, and the Block98 Food Frequency Questionnaire), body size (body mass index, waist and hip circumference, waist/hip ratio, and percent body fat), genetic variation in genes implicated related to the hypothalamic control of energy balance, and appropriate covariates (age, gender, race/ethnicity, smoking status, and physical activity. For the genetic association analyses, genotypes were collapsed by minor allele frequency, and haplotypes were estimated for each gene. Additionally, Bayesian networks were constructed in order to determine the relationships between genetic variation, eating behavior and body size.^ Results. We report that the EAT-26 score, Caloric intake, percent fat, fiber intake, HEAT index, and daily servings of vegetables, meats, grains, and fats were significantly associated with at least one body size measure. Multiple SNPs in 17 genes and haplotypes from 12 genes were tested for their association with body size. Variation within both DRD4 and HTR2A was found to be associated with EAT-26 score. In addition, variation in the ghrelin gene (GHRL) was significantly associated with daily Caloric intake. A significant interaction between daily servings of grains and the HEAT index and variation within the leptin receptor gene (LEPR) was shown to influence body size.^ Conclusion. This study has shown that there is a substantial genetic component to eating behavior and that genetic variation interacts with eating behavior to influence body size.^

Haplotype and nucleotide characteristics of photobiont and mycobiont gene sequences in the lichen Cetraria aculeata

Lichens, symbiotic associations of fungi (mycobionts) and green algae or cyanobacteria (photobionts), are poikilohydric organisms that are particularly well adapted to withstand adverse environmental conditions. Terrestrial ecosystems of the Antarctic are therefore largely dominated by lichens. The effects of global climate change are especially pronounced in the maritime Antarctic and it may be assumed that the lichen vegetation will profoundly change in the future. The genetic diversity of populations is closely correlated to their ability to adapt to changing environmental conditions and to their future evolutionary potential. In this study, we present evidence for low genetic diversity in Antarctic mycobiont and photobiont populations of the widespread lichen Cetraria aculeata. We compared between 110 and 219 DNA sequences from each of three gene loci for each symbiont. A total of 222 individuals from three Antarctic and nine antiboreal, temperate and Arctic populations were investigated. The mycobiont diversity is highest in Arctic populations, while the photobionts are most diverse in temperate regions. Photobiont diversity decreases significantly towards the Antarctic but less markedly towards the Arctic, indicating that ecological factors play a minor role in determining the diversity of Antarctic photobiont populations. Richness estimators calculated for the four geographical regions suggest that the low genetic diversity of Antarctic populations is not a sampling artefact. Cetraria aculeata appears to have diversified in the Arctic and subsequently expanded its range into the Southern Hemisphere. The reduced genetic diversity in the Antarctic is most likely due to founder effects during long-distance colonization.

Biometry, isotopic composition and diet of Arctic charr (Salvelinus alpinus) from Lake Aigneau, Northern Quebec

Two modal size groups of sexually mature Arctic charr (Salvelinus alpinus) differing in shape and found at different depths in Lake Aigneau in the Canadian sub-Arctic are described and tested for genetic and ecological differentiation. Forms consisted of a small littoral resident, mean size 21.7 cm, and a large profundal resident, mean size 53.9 cm. Mitochondrial DNA analysis indicated that seven of eight haplotypes were diagnostic for either the littoral or profundal fish, with 66.6% of the variation being found within form groupings. Pairwise tests of microsatellite data indicated significant differences in nine of 12 loci and a significant difference between the forms across all tested loci. Molecular variation was partitioned to 84.1% within and 15.9% between forms and suggestive of either restricted interbreeding over time or different allopatric origins. Stable isotope signatures were also significantly different, with the profundal fish having higher d13C and d15N values than the littoral fish. Overlap and separation, respectively, in the range of form d13C and d15N signatures indicated that carbon was obtained from similar sources, but that forms fed at different trophic levels. Littoral fish relied on aquatic insects, predominantly chironomids. Profundal fish were largely piscivorous, including cannibalism. Predominantly empty stomachs and low per cent nitrogen muscle-tissue composition among profundal fish further indicated that the feeding activity was limited to the winter when ice-cover increases the density of available prey at depth. Results provide evidence of significant differences between the modal groups, with origins in both genetics and ecology.

Mitochondrial DNA reveals hidden diversity for tardigrades from across the Antarctic realm

Antarctica contains some of the most challenging environmental conditions on the planet due to freezing temperatures, prolonged winters and lack of liquid water. Whereas 99.7% of Antarctica is permanently covered by ice and snow, some coastal areas and mountain ridges have remained ice-free and are able to sustain populations of microinvertebrates. Tardigrades are one of the more dominant groups of microfauna in soil and limno-terrestrial habitats, but little is known of their diversity and distribution across Antarctica. Here, we examine tardigrades sampled from across an extensive region of continental Antarctica, and analyse and compare their partial mitochondrial cytochrome c oxidase subunit 1 (COI) gene sequences with those from the Antarctic Peninsula, maritime and sub-Antarctica, Tierra del Fuego and other worldwide locations in order to recognise operational taxonomic units (OTUs). From 439 new tardigrade COI sequences, we identified 98 unique haplotypes (85 from Antarctica) belonging to Acutuncus, Diphascon, Echiniscus, Macrobiotus, Milnesium and unidentified Parachela. Operational taxonomic units were delimited by Poisson tree processes and general mixed Yule coalescent methods, resulting in 58 and 55 putative species, respectively. Most tardigrades appear to be locally endemic (i.e. restricted to a single geographic region), but some (e.g. Acutuncus antarcticus (Richters, 1904)) are widespread across continental Antarctica. Our molecular results reveal: (i) greater diversity than has previously been appreciated with distinct OTUs that potentially represent undescribed species, and (ii) a lack of connectivity between most OTUs from continental Antarctica and those from other Antarctic geographical zones.

Cladistic association analysis of Y chromosome effects on alcohol dependence and related personality traits

Association between Y chromosome haplotype variation and alcohol dependence and related personality traits was investigated in a large sample of psychiatrically diagnosed Finnish males. Haplotypes were constructed for 359 individuals using alleles at eight loci (seven microsatellite loci and a nucleotide substitution in the DYZ3 alphoid satellite locus). A cladogram linking the 102 observed haplotype configurations was constructed by using parsimony with a single-step mutation model. Then, a series of contingency tables nested according to the cladogram hierarchy were used to test for association between Y haplotype and alcohol dependence. Finally, using only alcohol-dependent subjects, we tested for association between Y haplotype and personality variables postulated to define subtypes of alcoholism—antisocial personality disorder, novelty seeking, harm avoidance, and reward dependence. Significant association with alcohol dependence was observed at three Y haplotype clades, with significance levels of P = 0.002, P = 0.020, and P = 0.010. Within alcohol-dependent subjects, no relationship was revealed between Y haplotype and antisocial personality disorder, novelty seeking, harm avoidance, or reward dependence. These results demonstrate, by using a fully objective association design, that differences among Y chromosomes contribute to variation in vulnerability to alcohol dependence. However, they do not demonstrate an association between Y haplotype and the personality variables thought to underlie the subtypes of alcoholism.

Differential expression of major histocompatibility complex class II genes on murine macrophages associated with T cell cytokine profile and protective/suppressive effects

Protective/suppressive major histocompatibility complex (MHC) class II alleles have been identified in humans and mice where they exert a disease-protective and immunosuppressive effect. Various modes of action have been proposed, among them differential expression of MHC class II genes in different types of antigen-presenting cells impacting on the T helper type 1 (Th1)–Th2 balance. To test this possibility, the expression of H-2 molecules from the four haplotypes H-2b, H-2d, H-2k, and H-2q was determined on bone marrow-derived macrophages (BMDMs) and splenic B cells. The I-Ab and I-Ek molecules, both well characterized as protective/suppressive, are expressed at a high level on almost all CD11b+ BMDMs for 5–8 days, after which expression slowly declines. In contrast, I-Ad, I-Ak, and I-Aq expression is lower, peaks over a shorter period, and declines more rapidly. No differential expression could be detected on B cells. In addition, the differential MHC class II expression found on macrophages skews the cytokine response of T cells as shown by an in vitro restimulation assay with BMDMs as antigen-presenting cells. The results indicate that macrophages of the protective/suppressive haplotypes express MHC class II molecules at a high level and exert Th1 bias, whereas low-level expression favors a Th2 response. We suggest that the extent of expression of the class II gene gates the back signal from T cells and in this way controls the activity of macrophages. This effect mediated by polymorphic nonexon segments of MHC class II genes may play a role in determining disease susceptibility in humans and mice.

Philopatry of male marine turtles inferred from mitochondrial DNA markers

Recent studies of mitochondrial DNA (mtDNA) variation among marine turtle populations are consistent with the hypothesis that females return to beaches in their natal region to nest as adults. In contrast, less is known about breeding migrations of male marine turtles and whether they too are philopatric to natal regions. Studies of geographic structuring of restriction fragment and microsatellite polymorphisms at anonymous nuclear loci in green turtle (Chelonia mydas) populations indicate that nuclear gene flow is higher than estimates from mtDNA analyses. Regional populations from the northern and southern Great Barrier Reef were distinct for mtDNA but indistinguishable at nuclear loci, whereas the Gulf of Carpentaria (northern Australia) population was distinct for both types of marker. To assess whether this result was due to reduced philopatry of males across the Great Barrier Reef, we determined the mtDNA haplotypes of breeding males at courtship areas for comparison with breeding females from the same three locations. We used a PCR-restriction fragment length polymorphism approach to determine control region haplotypes and designed mismatch primers for the identification of specific haplotypes. The mtDNA haplotype frequencies were not significantly different between males and females at any of the three areas and estimates of Fst among the regions were similar for males and females (Fst = 0.78 and 0.73, respectively). We conclude that breeding males, like females, are philopatric to courtship areas within their natal region. Nuclear gene flow between populations is most likely occurring through matings during migrations of both males and females through nonnatal courtship areas.

Phylogenetic assessment of length variation at a microsatellite locus

Sixty-six haplotypes at a locus containing a simple dinucleotide (CA)n microsatellite repeat were isolated by PCR–single-strand conformational polymorphism from populations of the horseshoe crab Limulus polyphemus. These haplotypes were sequenced to assess nucleotide variation directly. Thirty-four distinct sequences (alleles) were identified in a region 570 bp long that included the microsatellite motif. In the repeat region itself, CA-number varied in integer values from 5 to 11 across alleles, except that a (CA)8 class was not observed. Differences among alleles were due also to polymorphisms at 22 sites in regions immediately flanking the microsatellite repeats. Nucleotide substitutions in these regions were used to estimate phylogenetic relationships among alleles, and the gene phylogeny was used to trace the evolution of length variation and CA repeat numbers. A low correlation between size variation and genealogical relationships among alleles suggests that absolute fragment size (as normally scored in microsatellite assays) is an unreliable indicator of historical affinities among alleles. This finding on the molecular fine structure of microsatellite variation suggests the need for caution in the use of repeat counts at microsatellite loci as secure indicators of allelic relationships.