Evolution and development of male-specific leg brushes in Drosophilidae.

The origin, diversification, and secondary loss of sexually dimorphic characters are common in animal evolution. In some cases, structurally and functionally similar traits have evolved independently in multiple lineages. Prominent examples of such traits include the male-specific grasping structures that develop on the front legs of many dipteran insects. In this report, we describe the evolution and development of one of these structures, the male-specific "sex brush." The sex brush is composed of densely packed, irregularly arranged modified bristles and is found in several distantly related lineages in the family Drosophilidae. Phylogenetic analysis using 250 genes from over 200 species provides modest support for a single origin of the sex brush followed by many secondary losses; however, independent origins of the sex brush cannot be ruled out completely. We show that sex brushes develop in very similar ways in all brush-bearing lineages. The dense packing of brush hairs is explained by the specification of bristle precursor cells at a near-maximum density permitted by the lateral inhibition mechanism, as well as by the reduced size of the surrounding epithelial cells. In contrast to the female and the ancestral male condition, where bristles are arranged in stereotypical, precisely spaced rows, cell migration does not contribute appreciably to the formation of the sex brush. The complex phylogenetic history of the sex brush can make it a valuable model for investigating coevolution of sex-specific morphology and mating behavior.

Distinct genetic architectures underlie divergent thorax, leg, and wing pigmentation between Drosophila elegans and D. gunungcola.

Pigmentation divergence between Drosophila species has emerged as a model trait for studying the genetic basis of phenotypic evolution, with genetic changes contributing to pigmentation differences often mapping to genes in the pigment synthesis pathway and their regulators. These studies of Drosophila pigmentation have tended to focus on pigmentation changes in one body part for a particular pair of species, but changes in pigmentation are often observed in multiple body parts between the same pair of species. The similarities and differences of genetic changes responsible for divergent pigmentation in different body parts of the same species thus remain largely unknown. Here we compare the genetic basis of pigmentation divergence between Drosophila elegans and D. gunungcola in the wing, legs, and thorax. Prior work has shown that regions of the genome containing the pigmentation genes yellow and ebony influence the size of divergent male-specific wing spots between these two species. We find that these same two regions of the genome underlie differences in leg and thorax pigmentation; however, divergent alleles in these regions show differences in allelic dominance and epistasis among the three body parts. These complex patterns of inheritance can be explained by a model of evolution involving tissue-specific changes in the expression of Yellow and Ebony between D. elegans and D. gunungcola.

Author Correction: Genotypic variability-based genome-wide association study identifies non-additive loci HLA-C and IL12B for psoriasis.

This article was originally published under a CC BY-NC-SA License, but has now been made available under a CC BY 4.0 License.

Genome-wide association study of response to tumour necrosis factor inhibitor therapy in rheumatoid arthritis.

Rheumatoid arthritis (RA) is characterised by chronic synovial joint inflammation. Treatment has been revolutionised by tumour necrosis factor alpha inhibitors (TNFi) but each available drug shows a significant non-response rate. We conducted a genome-wide association study of 1752 UK RA TNFi-treated patients to identify predictors of change in the Disease Activity Score 28 (DAS28) and subcomponents over 3-6 months. The rs7195994 variant at the FTO gene locus was associated with infliximab response when looking at a change in the swollen joint count (SJC28) subcomponent (p = 9.74 × 10-9). Capture Hi-C data show chromatin interactions in GM12878 cells between rs2540767, in high linkage disequilibrium with rs7195994 (R2 = 0.9) and IRX3, a neighbouring gene of FTO. IRX3 encodes a transcription factor involved in adipocyte remodelling and is regarded as the obesity gene at the FTO locus. Importantly, the rs7195994 association remained significantly associated following adjustment for BMI. In addition, using capture Hi-C data we showed interactions between TNFi-response associated variants and 16 RA susceptibility variants.

Genome-wide association study of response to methotrexate in early rheumatoid arthritis patients.

Methotrexate (MTX) monotherapy is a common first treatment for rheumatoid arthritis (RA), but many patients do not respond adequately. In order to identify genetic predictors of response, we have combined data from two consortia to carry out a genome-wide study of response to MTX in 1424 early RA patients of European ancestry. Clinical endpoints were change from baseline to 6 months after starting treatment in swollen 28-joint count, tender 28-joint count, C-reactive protein and the overall 3-component disease activity score (DAS28). No single nucleotide polymorphism (SNP) reached genome-wide statistical significance for any outcome measure. The strongest evidence for association was with rs168201 in NRG3 (p = 10-7 for change in DAS28). Some support was also seen for association with ZMIZ1, previously highlighted in a study of response to MTX in juvenile idiopathic arthritis. Follow-up in two smaller cohorts of 429 and 177 RA patients did not support these findings, although these cohorts were more heterogeneous.

Genotypic variability-based genome-wide association study identifies non-additive loci HLA-C and IL12B for psoriasis.

Genome-wide association studies (GWASs) have identified a number of loci for psoriasis but largely ignored non-additive effects. We report a genotypic variability-based GWAS (vGWAS) that can prioritize non-additive loci without requiring prior knowledge of interaction types or interacting factors in two steps, using a mixed model to partition dichotomous phenotypes into an additive component and non-additive environmental residuals on the liability scale and then the Levene's (Brown-Forsythe) test to assess equality of the residual variances across genotype groups genome widely. The vGWAS identified two genome-wide significant (P < 5.0e-08) non-additive loci HLA-C and IL12B that were also genome-wide significant in an accompanying GWAS in the discovery cohort. Both loci were statistically replicated in vGWAS of an independent cohort with a small sample size. HLA-C and IL12B were reported in moderate gene-gene and/or gene-environment interactions in several occasions. We found a moderate interaction with age-of-onset of psoriasis, which was replicated indirectly. The vGWAS also revealed five suggestive loci (P < 6.76e-05) including FUT2 that was associated with psoriasis with environmental aspects triggered by virus infection and/or metabolic factors. Replication and functional investigation are needed to validate the suggestive vGWAS loci.

Analysis of SNP-SNP interactions and bone quantitative ultrasound parameter in early adulthood.

BACKGROUND: Osteoporosis individual susceptibility is determined by the interaction of multiple genetic variants and environmental factors. The aim of this study was to conduct SNP-SNP interaction analyses in candidate genes influencing heel quantitative ultrasound (QUS) parameter in early adulthood to identify novel insights into the mechanism of disease. METHODS: The study population included 575 healthy subjects (mean age 20.41; SD 2.36). To assess bone mass QUS was performed to determine Broadband ultrasound attenuation (BUA, dB/MHz). A total of 32 SNPs mapping to loci that have been characterized as genetic markers for QUS and/or BMD parameters were selected as genetic markers in this study. The association of all possible SNP pairs with QUS was assessed by linear regression and a SNP-SNP interaction was defined as a significant departure from additive effects. RESULTS: The pairwise SNP-SNP analysis showed multiple interactions. The interaction comprising SNPs rs9340799 and rs3736228 that map in the ESR1 and LRP5 genes respectively, revealed the lowest p value after adjusting for confounding factors (p-value = 0.001, ß (95% CI) = 14.289 (5.548, 23.029). In addition, our model reported others such as TMEM135-WNT16 (p = 0.007, ß(95%CI) = 9.101 (2.498, 15.704), ESR1-DKK1 (p = 0.012, ß(95%CI) = 13.641 (2.959, 24.322) or OPG-LRP5 (p = 0.012, ß(95%CI) = 8.724 (1.936, 15.512). However, none of the detected interactions remain significant considering the Bonferroni significance threshold for multiple testing (p<0.0001). CONCLUSION: Our analysis of SNP-SNP interaction in candidate genes of QUS in Caucasian young adults reveal several interactions, especially between ESR1 and LRP5 genes, that did not reach statistical significance. Although our results do not support a relevant genetic contribution of SNP-SNP epistatic interactions to QUS in young adults, further studies in larger independent populations would be necessary to support these preliminary findings.

Socially responsive effects of brain oxidative metabolism on aggression.

Despite ongoing high energetic demands, brains do not always use glucose and oxygen in a ratio that produces maximal ATP through oxidative phosphorylation. In some cases glucose consumption exceeds oxygen use despite adequate oxygen availability, a phenomenon known as aerobic glycolysis. Although metabolic plasticity seems essential for normal cognition, studying its functional significance has been challenging because few experimental systems link brain metabolic patterns to distinct behavioral states. Our recent transcriptomic analysis established a correlation between aggression and decreased whole-brain oxidative phosphorylation activity in the honey bee (Apis mellifera), suggesting that brain metabolic plasticity may modulate this naturally occurring behavior. Here we demonstrate that the relationship between brain metabolism and aggression is causal, conserved over evolutionary time, cell type-specific, and modulated by the social environment. Pharmacologically treating honey bees to inhibit complexes I or V in the oxidative phosphorylation pathway resulted in increased aggression. In addition, transgenic RNAi lines and genetic manipulation to knock down gene expression in complex I in fruit fly (Drosophila melanogaster) neurons resulted in increased aggression, but knockdown in glia had no effect. Finally, honey bee colony-level social manipulations that decrease individual aggression attenuated the effects of oxidative phosphorylation inhibition on aggression, demonstrating a specific effect of the social environment on brain function. Because decreased neuronal oxidative phosphorylation is usually associated with brain disease, these findings provide a powerful context for understanding brain metabolic plasticity and naturally occurring behavioral plasticity.

A role for HLA-DRB1*1101 and DRB1*0801 in cognitive ability and its decline with age.

Cognitive abilities (memory, processing speed, vocabulary, and fluid intelligence) are correlated with educational attainment and occupational status, as well as physical and mental health. The variation in cognitive abilities observed within a population has a substantial genetic contribution (heritability ∼50%) and yet the identification of genetic polymorphisms from both genome-wide association and candidate studies have to date only uncovered a limited number of genetic variants that exert small genetic effects. Here we impute human leukocyte antigens (HLA) using existing genome-wide association data from 1,559 non-pathological elderly volunteers who have been followed for changes in cognitive functioning between a 12- and 18-year period. Specifically, we investigate DRB1*05 (*11/*12) and DRB1*01, which have previously been associated with cognitive ability. We also analyze DRB1*0801, which shares close sequence homology with DRB1*1101. Together with DRB1*1101, DRB1*0801 has been associated with several diseases including multiple sclerosis and primary biliary cirrhosis, which themselves are associated with cognitive impairment. We observed that both DRB1*0801 and DRB1*1101 were significantly associated with vocabulary ability (cross-sectional and longitudinal scores) and that the effects were in opposite directions with DRB1*0801 associated with lower score and faster decline. This opposing affect is similar to that reported by other groups in systemic lupus erythematosus, type 1 diabetes, and primary biliary cirrhosis. DRB1*0801 was also significantly associated with reduced memory ability. We observed no associations between cognitive abilities and DRB1*01 or DRB1*12.

Clinical features of idiopathic inflammatory polymyopathy in the Hungarian Vizsla.

BACKGROUND: A retrospective study of the clinicopathological features of presumed and confirmed cases of idiopathic inflammatory polymyopathy in the Hungarian Vizsla dog and guidelines for breeding. RESULTS: 369 medical records were reviewed (1992-2013) and 77 Hungarian Vizslas were identified with a case history consistent with idiopathic inflammatory polymyopathy. Inclusion criteria were: group 1 (confirmed diagnosis); histopathology and clinical findings compatible with an inflammatory polymyopathy and group 2 (probable diagnosis); clinical findings compatible with a polymyopathy including dysphagia, sialorrhea, temporal muscle atrophy, elevated serum creatine kinase (CK) activity, and sufficient clinical history to suggest that other neuromuscular disorders could be ruled out. Some group 2 dogs had muscle biopsy, which suggested muscle disease but did not reveal an inflammatory process. The mean age of onset was 2.4 years; male dogs were slightly overrepresented. Common presenting signs were dysphagia, sialorrhea, masticatory muscle atrophy, and regurgitation. Common muscle histopathological findings included degenerative and regenerative changes, with multifocal mononuclear cell infiltration with lymphoplasmacytic myositis of variable severity. A positive response to immunosuppressive treatment supported an immune-mediated aetiology. The mean age at death and survival time were 6.4 and 3.9 years, respectively. Recurrence of clinical signs and aspiration pneumonia were common reasons for euthanasia. CONCLUSIONS: Diagnosis of Vizsla idiopathic inflammatory polymyopathy can be challenging due to lack of specific tests, however the presence of dysphagia, regurgitation and masticatory muscle atrophy in this breed with negative serological tests for masticatory muscle myositis and myasthenia gravis, along with muscle biopsies suggesting an inflammatory process, support the diagnosis. However, there is an urgent need for a more specific diagnostic test. The average of inbreeding coefficient (CoI) of 16.3% suggests an increased expression of a Dog Leukocyte Antigen Class II haplotype, leading to an increased disease risk. The prognosis remains guarded, as treatment can only manage the disease. Recurrence of clinical signs and perceived poor quality of life are the most common reasons for humane euthanasia.

Genetics of canine anal furunculosis in the German shepherd dog.

Canine anal furunculosis (AF) is characterised by ulceration and fistulation of perianal tissue and is a disease that particularly affects German shepherd dogs (GSDs). There are some similarities between AF and perianal Crohn's disease (CD) in man. An immune-mediated aetiopathogenesis for AF has been suggested due to tissue pathology, a major histocompatibility complex (MHC) association and clinical response to ciclosporin. Genome-wide association studies (GWAS) can be conducted in dogs with fewer markers and individuals than would be required in a human study. A discovery GWAS was performed on 21 affected and 25 control GSDs from the UK. No SNPs reached genome-wide significance levels at this stage. However, 127 nominally associated SNPs were genotyped in further 76 cases and 191 controls from the UK and Finland. Sequencing of these regions was undertaken to discover novel genetic variation. Association testing of these variants in the UK and Finnish cohorts revealed nine significantly associated SNPs, six of which cause non-synonymous changes in protein sequence. The ADAMTS16 and CTNND2 gene regions were most significantly associated with disease. Members of the butyrophilin protein family, important in intestinal inflammatory regulation, were also associated with disease, but their independence from the MHC region remains to be established. The CTNND2 gene region is also interesting as this locus was implicated in human ulcerative colitis and CD, albeit at a different candidate gene: DAP. We suggest that this represents a common association between inflammatory bowel disease-related conditions in both species and believe that future studies will strengthen this link.

MHC class II association study in eight breeds of dog with hypoadrenocorticism.

Canine hypoadrenocorticism is an endocrine disorder characterised by inadequate secretion of steroid hormones from the adrenal glands. Pathology results from immune-mediated destruction of the adrenal cortex, which is similar to that seen in the human Addison's disease. Both the canine and human diseases have similar clinical presentation, with the diagnosis based on performing a dynamic adrenocorticotropic hormone stimulation test. MHC class II has previously been associated with the human and canine diseases. In the current study, we conducted an MHC class II association study in eight breeds of dog with diagnoses of hypoadrenocorticism. We demonstrated significant differences in dog leukocyte antigen (DLA) haplotype frequencies in six of these breeds: Cocker spaniel, Springer spaniel, Labrador, West Highland white terrier (WHWT), Bearded collie, and Standard poodle. In the Springer spaniel, the DLA-DRB1*015:01--DQA1*006:01--DQB1*023:01 haplotype was significantly associated with disease risk (p = 0.014, odds ratio (OR) = 5.14) and showed a similar trend in the Cocker spaniel. This haplotype is related to one associated with hypoadrenocorticism in the Nova Scotia duck tolling retriever. Similar haplotypes shared between breeds were demonstrated, with DLA-DRB1*001:01--DQA1*001:01--DQB1*002:01 more prevalent in both affected Labrador (p = 0.0002, OR = 3.06) and WHWT (p = 0.01, OR = 2.11). Other haplotypes that have not previously been associated with the disease were identified. The inter-breed differences in DLA haplotypes associated with susceptibility to canine hypoadrenocorticism could represent divergent aetiologies. This could have implications for clinical diagnosis and future comparative studies. Alternatively, it may suggest that the gene of interest is closely linked to the MHC.

A candidate gene analysis of canine hypoadrenocorticism in 3 dog breeds.

Canine hypoadrenocorticism is believed to be an immune-related condition. It is rare in the overall dog population but shows a breed-related predisposition with Standard poodles and Portuguese water dogs having a greater prevalence of the condition. It shares many similarities with human primary adrenal insufficiency and is believed to be a naturally occurring, spontaneous model for the human condition. Short haplotype blocks and low levels of linkage disequilibrium in the human genome mean that the identification of genetic contributors to the condition requires large sample numbers. Pedigree dogs have high linkage disequilibrium and long haplotypes within a breed, increasing the potential of identifying novel genes that contribute to canine genetic disease. We investigated 222 SNPs from 42 genes that have been associated or may be implicated in human Addison's disease. We conducted case-control analyses in 3 pedigree dog breeds (Labrador retriever: affected n = 30, unaffected = 76; Cocker Spaniel: affected n = 19, unaffected = 53; Springer spaniel: affected n = 26, unaffected = 46) and identified 8 associated alleles in genes COL4A4, OSBPL9, CTLA4, PTPN22, and STXBP5 in 3 pedigree breeds. Association with immune response genes PTPN22 and CTLA4 in certain breeds suggests an underlying immunopathogenesis of the disease. These results suggest that canine hypoadrenocorticism could be a useful model for studying comparative genetics relevant to human Addison's disease.

Diversity and function of arthropod endosymbiont toxins.

Bacterial endosymbionts induce dramatic phenotypes in their arthropod hosts, including cytoplasmic incompatibility, feminization, parthenogenesis, male killing, parasitoid defense, and pathogen blocking. The molecular mechanisms underlying these effects remain largely unknown but recent evidence suggests that protein toxins secreted by the endosymbionts play a role. Here, we describe the diversity and function of endosymbiont proteins with homology to known bacterial toxins. We focus on maternally transmitted endosymbionts belonging to the Wolbachia, Rickettsia, Arsenophonus, Hamiltonella, Spiroplasma, and Cardinium genera because of their ability to induce the above phenotypes. We identify at least 16 distinct toxin families with diverse enzymatic activities, including AMPylases, nucleases, proteases, and glycosyltransferases. Notably, several annotated toxins contain domains with homology to eukaryotic proteins, suggesting that arthropod endosymbionts mimic host biochemistry to manipulate host physiology, similar to bacterial pathogens.

Co-evolving wing spots and mating displays are genetically separable traits in Drosophila.

The evolution of sexual traits often involves correlated changes in morphology and behavior. For example, in Drosophila, divergent mating displays are often accompanied by divergent pigment patterns. To better understand how such traits co-evolve, we investigated the genetic basis of correlated divergence in wing pigmentation and mating display between the sibling species Drosophila elegans and Drosophila gunungcola. Drosophila elegans males have an area of black pigment on their wings known as a wing spot and appear to display this spot to females by extending their wings laterally during courtship. By contrast, D. gunungcola lost both of these traits. Using Multiplexed Shotgun Genotyping (MSG), we identified a ∼440 kb region on the X chromosome that behaves like a genetic switch controlling the presence or absence of male-specific wing spots. This region includes the candidate gene optomotor-blind (omb), which plays a critical role in patterning the Drosophila wing. The genetic basis of divergent wing display is more complex, with at least two loci on the X chromosome and two loci on autosomes contributing to its evolution. Introgressing the X-linked region affecting wing spot development from D. gunungcola into D. elegans reduced pigmentation in the wing spots but did not affect the wing display, indicating that these are genetically separable traits. Consistent with this observation, broader sampling of wild D. gunungcola populations confirmed that the wing spot and wing display are evolving independently: some D. gunungcola males performed wing displays similar to D. elegans despite lacking wing spots. These data suggest that correlated selection pressures rather than physical linkage or pleiotropy are responsible for the coevolution of these morphological and behavioral traits. They also suggest that the change in morphology evolved prior to the change in behavior.

The yellow gene influences Drosophila male mating success through sex comb melanization.

Drosophila melanogaster males perform a series of courtship behaviors that, when successful, result in copulation with a female. For over a century, mutations in the yellow gene, named for its effects on pigmentation, have been known to reduce male mating success. Prior work has suggested that yellow influences mating behavior through effects on wing extension, song, and/or courtship vigor. Here, we rule out these explanations, as well as effects on the nervous system more generally, and find instead that the effects of yellow on male mating success are mediated by its effects on pigmentation of male-specific leg structures called sex combs. Loss of yellow expression in these modified bristles reduces their melanization, which changes their structure and causes difficulty grasping females prior to copulation. These data illustrate why the mechanical properties of anatomy, not just neural circuitry, must be considered to fully understand the development and evolution of behavior.

Pleiotropic Effects of ebony and tan on Pigmentation and Cuticular Hydrocarbon Composition in Drosophila melanogaster.

Pleiotropic genes are genes that affect more than one trait. For example, many genes required for pigmentation in the fruit fly Drosophila melanogaster also affect traits such as circadian rhythms, vision, and mating behavior. Here, we present evidence that two pigmentation genes, ebony and tan, which encode enzymes catalyzing reciprocal reactions in the melanin biosynthesis pathway, also affect cuticular hydrocarbon (CHC) composition in D. melanogaster females. More specifically, we report that ebony loss-of-function mutants have a CHC profile that is biased toward long (>25C) chain CHCs, whereas tan loss-of-function mutants have a CHC profile that is biased toward short (<25C) chain CHCs. Moreover, pharmacological inhibition of dopamine synthesis, a key step in the melanin synthesis pathway, reversed the changes in CHC composition seen in ebony mutants, making the CHC profiles similar to those seen in tan mutants. These observations suggest that genetic variation affecting ebony and/or tan activity might cause correlated changes in pigmentation and CHC composition in natural populations. We tested this possibility using the Drosophila Genetic Reference Panel (DGRP) and found that CHC composition covaried with pigmentation as well as levels of ebony and tan expression in newly eclosed adults in a manner consistent with the ebony and tan mutant phenotypes. These data suggest that the pleiotropic effects of ebony and tan might contribute to covariation of pigmentation and CHC profiles in Drosophila.

Evaluation of metrics for benchmarking antimicrobial use in the UK dairy industry.

The issue of antimicrobial resistance is of global concern across human and animal health. In 2016, the UK government committed to new targets for reducing antimicrobial use (AMU) in livestock. Although a number of metrics for quantifying AMU are defined in the literature, all give slightly different interpretations. This paper evaluates a selection of metrics for AMU in the dairy industry: total mg, total mg/kg, daily dose and daily course metrics. Although the focus is on their application to the dairy industry, the metrics and issues discussed are relevant across livestock sectors. In order to be used widely, a metric should be understandable and relevant to the veterinarians and farmers who are prescribing and using antimicrobials. This means that clear methods, assumptions (and possible biases), standardised values and exceptions should be published for all metrics. Particularly relevant are assumptions around the number and weight of cattle at risk of treatment and definitions of dose rates and course lengths; incorrect assumptions can mean metrics over-represent or under-represent AMU. The authors recommend that the UK dairy industry work towards the UK-specific metrics using the UK-specific medicine dose and course regimens as well as cattle weights in order to monitor trends nationally.