First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.     

Infectious diseases as main causes of mortality to beavers Castor fiber after translocation to the Netherlands

Between 1988 and 1994, 58 beavers were translocated from the Elbe region (Germany) to the Netherlands. In 43 animals, radio-transmitters were implanted with a pulse interval which was dependent on body temperature; subsequently, 22 of the released animals were found dead and the cause of death was investigated. Infectious diseases were the most important cause of death (50%), in particular yersiniosis and leptospirosis. Yersiniosis is relatively rare in the source population along the Elbe, and leptospirosis has not been reported before in beavers. We hypothesize that the immune system of the translocated animals was weakened by stress at a time of exposure to new pathogens. In future translocations of beavers, we recommend to reduce the stress, improve the hygiene, and to vaccinate against Yersinia pseudotuberculosis and leptospirosis prior to the release.     

Methicillin-resistant Staphylococcus aureus Bacterial Nitric-oxide Synthase Affects Antibiotic Sensitivity and Skin Abscess Development

Staphylococcus aureus infections present an enormous global health concern complicated by an alarming increase in antibiotic resistance. S. aureus is among the few bacterial species that express nitric-oxide synthase (bNOS) and thus can catalyze NO production from l-arginine. Here we generate an isogenic bNOS-deficient mutant in the epidemic community-acquired methicillin-resistant S. aureus (MRSA) USA300 clone to study its contribution to virulence and antibiotic susceptibility. Loss of bNOS increased MRSA susceptibility to reactive oxygen species and host cathelicidin antimicrobial peptides, which correlated with increased MRSA killing by human neutrophils and within neutrophil extracellular traps. bNOS also promoted resistance to the pharmaceutical antibiotics that act on the cell envelope such as vancomycin and daptomycin. Surprisingly, bNOS-deficient strains gained resistance to aminoglycosides, suggesting that the role of bNOS in antibiotic susceptibility is more complex than previously observed in Bacillus species. Finally, the MRSA bNOS mutant showed reduced virulence with decreased survival and smaller abscess generation in a mouse subcutaneous infection model. Together, these data indicate that bNOS contributes to MRSA innate immune and antibiotic resistance phenotypes. Future development of specific bNOS inhibitors could be an attractive option to simultaneously reduce MRSA pathology and enhance its susceptibility to commonly used antibiotics.     

The Himalayas: Barrier and conduit for gene flow

The Himalayan mountain range is strategically located at the crossroads of the major cultural centers in Asia, the Middle East and Europe. Although previous Y‐chromosome studies indicate that the Himalayas served as a natural barrier for gene flow from the south to the Tibetan plateau, this region is believed to have played an important role as a corridor for human migrations between East and West Eurasia along the ancient Silk Road. To evaluate the effects of the Himalayan mountain range in shaping the maternal lineages of populations residing on either side of the cordillera, we analyzed mitochondrial DNA variation in 344 samples from three Nepalese collections (Newar, Kathmandu and Tamang) and a general population of Tibet. Our results revealed a predominantly East Asian‐specific component in Tibet and Tamang, whereas Newar and Kathmandu are both characterized by a combination of East and South Central Asian lineages. Interestingly, Newar and Kathmandu harbor several deep‐rooted Indian lineages, including M2, R5, and U2, whose coalescent times from this study (U2, >40 kya) and previous reports (M2 and R5, >50 kya) suggest that Nepal was inhabited during the initial peopling of South Central Asia. Comparisons with our previous Y‐chromosome data indicate sex‐biased migrations in Tamang and a founder effect and/or genetic drift in Tamang and Newar. Altogether, our results confirm that while the Himalayas acted as a geographic barrier for human movement from the Indian subcontinent to the Tibetan highland, it also served as a conduit for gene flow between Central and East Asia. Am J Phys Anthropol 151:169–182, 2013. © 2013 Wiley Periodicals, Inc.     

Dealing with paralogy in RADseq data: in silico detection and single nucleotide polymorphism validation in Robinia pseudoacacia L.

The RADseq technology allows researchers to efficiently develop thousands of polymorphic loci across multiple individuals with little or no prior information on the genome. However, many questions remain about the biases inherent to this technology. Notably, sequence misalignments arising from paralogy may affect the development of single nucleotide polymorphism (SNP) markers and the estimation of genetic diversity. We evaluated the impact of putative paralog loci on genetic diversity estimation during the development of SNPs from a RADseq dataset for the nonmodel tree species Robinia pseudoacacia L. We sequenced nine genotypes and analyzed the frequency of putative paralogous RAD loci as a function of both the depth of coverage and the mismatch threshold allowed between loci. Putative paralogy was detected in a very variable number of loci, from 1% to more than 20%, with the depth of coverage having a major influence on the result. Putative paralogy artificially increased the observed degree of polymorphism and resulting estimates of diversity. The choice of the depth of coverage also affected diversity estimation and SNP validation: A low threshold decreased the chances of detecting minor alleles while a high threshold increased allelic dropout. SNP validation was better for the low threshold (4×) than for the high threshold (18×) we tested. Using the strategy developed here, we were able to validate more than 80% of the SNPs tested by means of individual genotyping, resulting in a readily usable set of 330 SNPs, suitable for use in population genetics applications.     

Decorin reverses the repressive effect of autocrine-produced TGF-beta on mouse macrophage activation

Several cytokines or growth factors induce macrophages to proliferate, become activated, differentiate, or die through apoptosis. Like the major macrophage activator IFN-gamma, the extracellular matrix protein decorin inhibits proliferation and protects macrophages from the induction of apoptosis. Decorin enhances the IFN-gamma-induced expression of the IAalpha and IAbeta MHC class II genes. Moreover, it increases the IFN-gamma- or LPS-induced expression of inducible NO synthase, TNF-alpha, IL-1beta, and IL-6 genes and the secretion of these cytokines. Using a number of extracellular matrix proteins, we found a negative correlation between adhesion and proliferation. However, the effects of decorin on macrophage activation do not seem to be mediated through its effect on adhesion or proliferation. Instead, this proteoglycan abolishes the binding of TGF-beta to macrophages, as shown by Scatchard analysis of (125)I-labeled TGF-beta, which, in the absence of decorin, showed a K(d) of 0.11 +/- 0.03 nM and approximately 5000 receptors/cell. This was confirmed when we treated macrophages with Abs to block the endogenously produced TGF-beta, which enhanced macrophage activation in a way similar to decorin. The increase in activation mediated by decorin demonstrates that macrophages are under negative regulation that can be reversed by proteins of the extracellular matrix.     

Structure of oxalate decarboxylase from Bacillus subtilis at 1.75 A resolution

Oxalate decarboxylase is a manganese-dependent enzyme that catalyzes the conversion of oxalate to formate and carbon dioxide. We have determined the structure of oxalate decarboxylase from Bacillus subtilis at 1.75 A resolution in the presence of formate. The structure reveals a hexamer with 32-point symmetry in which each monomer belongs to the cupin family of proteins. Oxalate decarboxylase is further classified as a bicupin because it contains two cupin folds, possibly resulting from gene duplication. Each oxalate decarboxylase cupin domain contains one manganese binding site. Each of the oxalate decarboxylase domains is structurally similar to oxalate oxidase, which catalyzes the manganese-dependent oxidative decarboxylation of oxalate to carbon dioxide and hydrogen peroxide. Amino acid side chains in the two metal binding sites of oxalate decarboxylase and the metal binding site of oxalate oxidase are very similar. Four manganese binding residues (three histidines and one glutamate) are conserved as well as a number of hydrophobic residues. The most notable difference is the presence of Glu333 in the metal binding site of the second cupin domain of oxalate decarboxylase. We postulate that this domain is responsible for the decarboxylase activity and that Glu333 serves as a proton donor in the production of formate. Mutation of Glu333 to alanine reduces the catalytic activity by a factor of 25. The function of the other domain in oxalate decarboxylase is not yet known.     

Target selection and determination of function in structural genomics

The first crucial step in any structural genomics project is the selection and prioritization of target proteins for structure determination. There may be a number of selection criteria to be satisfied, including that the proteins have novel folds, that they be representatives of large families for which no structure is known, and so on. The better the selection at this stage, the greater is the value of the structures obtained at the end of the experimental process. This value can be further enhanced once the protein structures have been solved if the functions of the given proteins can also be determined. Here we describe the methods used at either end of the experimental process: firstly, sensitive sequence comparison techniques for selecting a high-quality list of target proteins, and secondly the various computational methods that can be applied to the eventual 3D structures to determine the most likely biochemical function of the proteins in question.