Genetic diversity among Mycoplasma species bovine group 7: clonal isolates from an outbreak of polyarthritis, mastitis, and abortion in dairy cattle

A comprehensive genetic analysis of 60 Mycoplasma sp. bovine group 7 isolates from different geographic origins and epidemiological settings is presented. Twenty-four isolates were recovered from the joints of calves during sporadic episodes of polyarthritis in geographically distinct regions of Queensland and New South Wales, Australia, including two clones of the type strain PG5O. A further three Australian isolates were also recovered from the tympanic bulla, retropharyngeal lymph node and the lung and another three isolates had unconfirmed histories. Six isolates originated from Germany, Portugal, Nigeria, and France. Twenty-four epidemiologically related isolates of Mycoplasma sp. bovine group 7 were recovered from multiple tissue sites and body fluids of infected calves with polyarthritis, mastitic milk, and from the stomach contents, lung and liver from aborted foetuses in three large, centrally managed dairy herds in New South Wales, Australia. Restriction endonuclease analysis (REA) of genomic DNA differentiated 29 Cfol profiles among these 60 isolates and grouped all 24 epidemiologically related isolates in a defined pattern showing a clonal origin. Three isolates of this clonal cluster were recovered from mastitic milk and the synovial exudate of clinically-affected calves and appeared sporadically for periods up to 18 months after the initial outbreak of polyarthritis indicating a persistent, close association of the organism with cattle in these herds. The Cfol profile representative of the clonal cluster was distinguishable from profiles of isolates recovered from multiple, unrelated cases of polyarthritis in Queensland and New South Wales and from other countries. All 24 isolates from the clonal cluster possessed a plasmid (pBG7AU) with a molecular size of 1022 bp. DNA sequence analysis of pBG7AU identified two open reading frames sharing 81 and 99% DNA sequence similarity with hypothetical replication control proteins A and B respectively, previously described in plasmid pADB201 isolated from M. mycoides subspecies mycoides. Other isolates of bovine group 7, epidemiologically unrelated to the clonal cluster, including two clones of the type strain PG5O, possessed a similar-sized plasmid. These data confirm that Mycoplasma sp. bovine group 7 is capable of migrating to, and multiplying within, different tissue sites within a single animal and among different animals within a herd.

Copy number variation of the Beta-defensin genes in europeans: no supporting evidence for association with lung function, chronic obstructive pulmonary disease or asthma

Lung function measures are heritable, predict mortality and are relevant in diagnosis of chronic obstructive pulmonary disease (COPD). COPD and asthma are diseases of the airways with major public health impacts and each have a heritable component. Genome-wide association studies of SNPs have revealed novel genetic associations with both diseases but only account for a small proportion of the heritability. Complex copy number variation may account for some of the missing heritability. A well-characterised genomic region of complex copy number variation contains beta-defensin genes (DEFB103, DEFB104 and DEFB4), which have a role in the innate immune response. Previous studies have implicated these and related genes as being associated with asthma or COPD. We hypothesised that copy number variation of these genes may play a role in lung function in the general population and in COPD and asthma risk. We undertook copy number typing of this locus in 1149 adult and 689 children using a paralogue ratio test and investigated association with COPD, asthma and lung function. Replication of findings was assessed in a larger independent sample of COPD cases and smoking controls. We found evidence for an association of beta-defensin copy number with COPD in the adult cohort (OR = 1.4, 95%CI:1.02-1.92, P = 0.039) but this finding, and findings from a previous study, were not replicated in a larger follow-up sample(OR = 0.89, 95%CI:0.72-1.07, P = 0.217). No robust evidence of association with asthma in children was observed. We found no evidence for association between beta-defensin copy number and lung function in the general populations. Our findings suggest that previous reports of association of beta-defensin copy number with COPD should be viewed with caution. Suboptimal measurement of copy number can lead to spurious associations. Further beta-defensin copy number measurement in larger sample sizes of COPD cases and children with asthma are needed.

Genetic characterization of antimicrobial resistance in coagulase-negative staphylococci from bovine mastitis milk

Coagulase-negative staphylococci (CNS; n=417) were isolated from bovine milk and identified by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Nineteen different species were identified, and Staphylococcus xylosus, Staphylococcus chromogenes, Staphylococcus haemolyticus, and Staphylococcus sciuri were the most prevalent species. Resistance to oxacillin (47.0% of the isolates), fusidic acid (33.8%), tiamulin (31.9%), penicillin (23.3%), tetracycline (15.8%), streptomycin (9.6%), erythromycin (7.0%), sulfonamides (5%), trimethoprim (4.3%), clindamycin (3.4%), kanamycin (2.4%), and gentamicin (2.4%) was detected. Resistance to oxacillin was attributed to the mecA gene in 9.7% of the oxacillin-resistant isolates. The remaining oxacillin-resistant CNS did not contain the mecC gene or mecA1 promoter mutations. The mecA gene was detected in Staphylococcus fleurettii, Staphylococcus epidermidis, Staph. haemolyticus, and Staph. xylosus. Resistance to tetracycline was attributed to the presence of tet(K) and tet(L), penicillin resistance to blaZ, streptomycin resistance to str and ant(6)-Ia, and erythromycin resistance to erm(C), erm(B), and msr. Resistance to tiamulin and fusidic acid could not be attributed to an acquired resistance gene. In total, 15.1% of the CNS isolates were multidrug resistant (i.e., resistant to 2 or more antimicrobials). The remaining CNS isolates were susceptible to antimicrobials commonly used in mastitis treatment. Methicillin-resistant CNS isolates were diverse, as determined by mecA gene sequence analysis, staphylococcal cassette chromosome mec typing, and pulsed-field gel electrophoresis. Arginine catabolic mobile element types 1 and 3 were detected in both methicillin-resistant and methicillin-susceptible Staph. epidermidis and were associated with sequence types ST59 and ST111. Because this study revealed the presence of multidrug-resistant CNS in a heterogeneous CNS population, we recommend antibiogram analysis of CNS in persistent infections before treatment with antimicrobials.

Genome-wide analysis in German shepherd dogs reveals association of a locus on CFA 27 with atopic dermatitis

Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1 × 10(-5)) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n cases = 91, n controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0 × 10(-9)), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (praw = 3.1 × 10(-7), pgenome = 0.03). The total associated region was defined as a ~1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ~209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.

A genome-wide association study to detect QTL for commercially important traits in Swiss Large White boars

The improvement of meat quality and production traits has high priority in the pork industry. Many of these traits show a low to moderate heritability and are difficult and expensive to measure. Their improvement by targeted breeding programs is challenging and requires knowledge of the genetic and molecular background. For this study we genotyped 192 artificial insemination boars of a commercial line derived from the Swiss Large White breed using the PorcineSNP60 BeadChip with 62,163 evenly spaced SNPs across the pig genome. We obtained 26 estimated breeding values (EBVs) for various traits including exterior, meat quality, reproduction, and production. The subsequent genome-wide association analysis allowed us to identify four QTL with suggestive significance for three of these traits (p-values ranging from 4.99×10⁻⁶ to 2.73×10⁻⁵). Single QTL for the EBVs pH one hour post mortem (pH1) and carcass length were on pig chromosome (SSC) 14 and SSC 2, respectively. Two QTL for the EBV rear view hind legs were on SSC 10 and SSC 16.

Genetic diversity in the modern horse illustrated from genome-wide SNP data.

Horses were domesticated from the Eurasian steppes 5,000-6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. F(ST) calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection.

Genomewide association for a dominant pigmentation gene in sheep

Most published genomewide association studies (GWAS) in sheep have investigated recessively inherited monogenic traits. The objective here was to assess the feasibility of performing GWAS for a dominant trait for which the genetic basis was already known. A total of 42 Manchega and Rasa Aragonesa sheep that segregate solid black or white coat pigmentation were genotyped using the SNP50 BeadChip. Previous analysis in Manchegas demonstrated a complete association between the pigmentation trait and alleles of the MC1R gene, setting an a priori expectation for GWAS. Multiple methods were used to identify and quantify the strength of population substructure between black and white animals, before allelic association testing was performed for 49 034 SNPs. Following correction for substructure, GWAS identified the most strongly associated SNP (s26449) was also the closest to the MC1R gene. The finding was strongly supported by the permutation tree-based random forest (RF) analysis. Importantly, GWAS identified unlinked SNP with only slightly lower p-values than for s26449. Random forest analysis indicated these were false positives, suggesting interpretation based on both approaches was beneficial. The results indicate that a combined analytical approach can be successful in studies where a modest number of animals are available and substantial population stratification exists.

A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control

HIV-1 sequence diversity is affected by selection pressures arising from host genomic factors. Using paired human and viral data from 1071 individuals, we ran >3000 genome-wide scans, testing for associations between host DNA polymorphisms, HIV-1 sequence variation and plasma viral load (VL), while considering human and viral population structure. We observed significant human SNP associations to a total of 48 HIV-1 amino acid variants (p<2.4 × 10−12). All associated SNPs mapped to the HLA class I region. Clinical relevance of host and pathogen variation was assessed using VL results. We identified two critical advantages to the use of viral variation for identifying host factors: (1) association signals are much stronger for HIV-1 sequence variants than VL, reflecting the ‘intermediate phenotype’ nature of viral variation; (2) association testing can be run without any clinical data. The proposed genome-to-genome approach highlights sites of genomic conflict and is a strategy generally applicable to studies of host–pathogen interaction.

The Individualized Genetic Barrier Predicts Treatment Response in a Large Cohort of HIV-1 Infected Patients

The success of combination antiretroviral therapy is limited by the evolutionary escape dynamics of HIV-1. We used Isotonic Conjunctive Bayesian Networks (I-CBNs), a class of probabilistic graphical models, to describe this process. We employed partial order constraints among viral resistance mutations, which give rise to a limited set of mutational pathways, and we modeled phenotypic drug resistance as monotonically increasing along any escape pathway. Using this model, the individualized genetic barrier (IGB) to each drug is derived as the probability of the virus not acquiring additional mutations that confer resistance. Drug-specific IGBs were combined to obtain the IGB to an entire regimen, which quantifies the virus' genetic potential for developing drug resistance under combination therapy. The IGB was tested as a predictor of therapeutic outcome using between 2,185 and 2,631 treatment change episodes of subtype B infected patients from the Swiss HIV Cohort Study Database, a large observational cohort. Using logistic regression, significant univariate predictors included most of the 18 drugs and single-drug IGBs, the IGB to the entire regimen, the expert rules-based genotypic susceptibility score (GSS), several individual mutations, and the peak viral load before treatment change. In the multivariate analysis, the only genotype-derived variables that remained significantly associated with virological success were GSS and, with 10-fold stronger association, IGB to regimen. When predicting suppression of viral load below 400 cps/ml, IGB outperformed GSS and also improved GSS-containing predictors significantly, but the difference was not significant for suppression below 50 cps/ml. Thus, the IGB to regimen is a novel data-derived predictor of treatment outcome that has potential to improve the interpretation of genotypic drug resistance tests.

Contribution of Genetic Background, Traditional Risk Factors, and HIV-Related Factors to Coronary Artery Disease Events in HIV-Positive Persons

Background Persons infected with human immunodeficiency virus (HIV) have increased rates of coronary artery disease (CAD). The relative contribution of genetic background, HIV-related factors, antiretroviral medications, and traditional risk factors to CAD has not been fully evaluated in the setting of HIV infection. Methods In the general population, 23 common single-nucleotide polymorphisms (SNPs) were shown to be associated with CAD through genome-wide association analysis. Using the Metabochip, we genotyped 1875 HIV-positive, white individuals enrolled in 24 HIV observational studies, including 571 participants with a first CAD event during the 9-year study period and 1304 controls matched on sex and cohort. Results A genetic risk score built from 23 CAD-associated SNPs contributed significantly to CAD (P = 2.9×10−4). In the final multivariable model, participants with an unfavorable genetic background (top genetic score quartile) had a CAD odds ratio (OR) of 1.47 (95% confidence interval [CI], 1.05–2.04). This effect was similar to hypertension (OR = 1.36; 95% CI, 1.06–1.73), hypercholesterolemia (OR = 1.51; 95% CI, 1.16–1.96), diabetes (OR = 1.66; 95% CI, 1.10–2.49), ≥1 year lopinavir exposure (OR = 1.36; 95% CI, 1.06–1.73), and current abacavir treatment (OR = 1.56; 95% CI, 1.17–2.07). The effect of the genetic risk score was additive to the effect of nongenetic CAD risk factors, and did not change after adjustment for family history of CAD. Conclusions In the setting of HIV infection, the effect of an unfavorable genetic background was similar to traditional CAD risk factors and certain adverse antiretroviral exposures. Genetic testing may provide prognostic information complementary to family history of CAD.

GENETIC ANALYSIS OF THE FUNCTION OF THE DROSOPHILA DOUBLESEX-RELATED FACTOR dmrt93B

DMRT (Doublesex and Mab-3 related transcription factor) proteins generally associated with sexual differentiation in many organisms share a common DNA binding domain and are often expressed in reproductive tissues. Aside from doublesex, which is a central factor in the regulation of sex determination, Drosophila possesses three different dmrt genes that are of unknown function. Because the association with sexual differentiation and reproduction is not universal and some DMRT proteins have been found to play other developmental roles we chose to further characterize one of these Drosophila genes. We carried out genetic analysis of dmrt93B, which was previously found to be expressed sex-specifically in the developing somatic gonad and to affect testis morphogenesis in RNAi knockdowns. In order to disrupt this gene, the GAL4 yeast transcriptional activator followed by a polyadenylation signal was inserted after the dmrt93B start codon and introduced into the genome by homologous recombination. Analysis of the knock-in mutation as well as a small deletion removing all dmrt93B sequence demonstrate that loss of function causes partial lethality at the late pupal stage. Surprisingly, these mutations have no significant effect on gonad formation or male fertility. Analysis of GAL4-driven GFP reporter expression indicates that the dmrt93B promoter activity is highly specific to neurons in the suboesophageal and proventricular ganglion in larva and adult of both sexes suggesting a possible role in digestive tract function. Using the Capillary Feeder (CAFÉ) assay to measure daily food intake we find that reduction in this gene’s function leads to an increase in food consumption. These results suggest dmrt93 plays an important role in the formation or maintenance of neurons that affect feeding and support the idea that dmrt genes may not be restricted to roles in sexual differentiation.

Detection of nonrandom association of alleles from the distribution of the number of heterozygous loci in a sample.

The distribution of the number of heterozygous loci in two randomly chosen gametes or in a random diploid zygote provides information regarding the nonrandom association of alleles among different genetic loci. Two alternative statistics may be employed for detection of nonrandom association of genes of different loci when observations are made on these distributions: observed variance of the number of heterozygous loci (s2k) and a goodness-of-fit criterion (X2) to contrast the observed distribution with that expected under the hypothesis of random association of genes. It is shown, by simulation, that s2k is statistically more efficient than X2 to detect a given extent of nonrandom association. Asymptotic normality of s2k is justified, and X2 is shown to follow a chi-square (chi 2) distribution with partial loss of degrees of freedom arising because of estimation of parameters from the marginal gene frequency data. Whenever direct evaluations of linkage disequilibrium values are possible, tests based on maximum likelihood estimators of linkage disequilibria require a smaller sample size (number of zygotes or gametes) to detect a given level of nonrandom association in comparison with that required if such tests are conducted on the basis of s2k. Summarization of multilocus genotype (or haplotype) data, into the different number of heterozygous loci classes, thus, amounts to appreciable loss of information.

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Multiple independent genetic factors at NOS1AP modulate the QT interval in a multi-ethnic population.

Extremes of electrocardiographic QT interval are associated with increased risk for sudden cardiac death (SCD); thus, identification and characterization of genetic variants that modulate QT interval may elucidate the underlying etiology of SCD. Previous studies have revealed an association between a common genetic variant in NOS1AP and QT interval in populations of European ancestry, but this finding has not been extended to other ethnic populations. We sought to characterize the effects of NOS1AP genetic variants on QT interval in the multi-ethnic population-based Dallas Heart Study (DHS, n = 3,072). The SNP most strongly associated with QT interval in previous samples of European ancestry, rs16847548, was the most strongly associated in White (P = 0.005) and Black (P = 3.6 x 10(-5)) participants, with the same direction of effect in Hispanics (P = 0.17), and further showed a significant SNP x sex-interaction (P = 0.03). A second SNP, rs16856785, uncorrelated with rs16847548, was also associated with QT interval in Blacks (P = 0.01), with qualitatively similar results in Whites and Hispanics. In a previously genotyped cohort of 14,107 White individuals drawn from the combined Atherosclerotic Risk in Communities (ARIC) and Cardiovascular Health Study (CHS) cohorts, we validated both the second locus at rs16856785 (P = 7.63 x 10(-8)), as well as the sex-interaction with rs16847548 (P = 8.68 x 10(-6)). These data extend the association of genetic variants in NOS1AP with QT interval to a Black population, with similar trends, though not statistically significant at P<0.05, in Hispanics. In addition, we identify a strong sex-interaction and the presence of a second independent site within NOS1AP associated with the QT interval. These results highlight the consistent and complex role of NOS1AP genetic variants in modulating QT interval.

A genome-wide association study of resistance to HIV infection in highly exposed uninfected individuals with hemophilia A.

Human genetic variation contributes to differences in susceptibility to HIV-1 infection. To search for novel host resistance factors, we performed a genome-wide association study (GWAS) in hemophilia patients highly exposed to potentially contaminated factor VIII infusions. Individuals with hemophilia A and a documented history of factor VIII infusions before the introduction of viral inactivation procedures (1979-1984) were recruited from 36 hemophilia treatment centers (HTCs), and their genome-wide genetic variants were compared with those from matched HIV-infected individuals. Homozygous carriers of known CCR5 resistance mutations were excluded. Single nucleotide polymorphisms (SNPs) and inferred copy number variants (CNVs) were tested using logistic regression. In addition, we performed a pathway enrichment analysis, a heritability analysis, and a search for epistatic interactions with CCR5 Δ32 heterozygosity. A total of 560 HIV-uninfected cases were recruited: 36 (6.4%) were homozygous for CCR5 Δ32 or m303. After quality control and SNP imputation, we tested 1 081 435 SNPs and 3686 CNVs for association with HIV-1 serostatus in 431 cases and 765 HIV-infected controls. No SNP or CNV reached genome-wide significance. The additional analyses did not reveal any strong genetic effect. Highly exposed, yet uninfected hemophiliacs form an ideal study group to investigate host resistance factors. Using a genome-wide approach, we did not detect any significant associations between SNPs and HIV-1 susceptibility, indicating that common genetic variants of major effect are unlikely to explain the observed resistance phenotype in this population.