In insects, airborne chemical signals are mainly detected by two receptor families, odorant receptors (ORs) and ionotropic receptors (IRs). Functions of ORs have been intensively investigated in Diptera and Lepidoptera, while the functions and evolution of the more ancient IR family remain largely unexplored beyond Diptera. Here, we identified a repertoire of 26 IRs from transcriptomes of female and male antennae, and ovipositors in the moth Agrotis segetum. We observed that a large clade formed by IR75p and IR75q expansions is closely related to the acid-sensing IRs identified in Diptera. We functionally assayed each of the five AsegIRs from this clade using Xenopus oocytes and found that two receptors responded to the tested ligands. AsegIR75p.1 responded to several compounds but hexanoic acid was revealed to be the primary ligand, and AsegIR75q.1 responded primarily to octanoic acid, and less so to nonanoic acid. It has been reported that the C6-C10 medium-chain fatty acids repel various insects including many drosophilids and mosquitos. We show that the C6-C10 medium-chain fatty acids elicited antennal responses of both sexes of A. segetum, while only octanoic acid had repellent effect to the moths in a behavioral assay. In addition, using fluorescence in situ hybridization, we demonstrated that the five IRs and their co-receptor AsegIR8a are not located in coeloconic sensilla as found in Drosophila, but in basiconic or trichoid sensilla. Our results significantly expand the current knowledge of the insect IR family. Based on the functional data in combination with phylogenetic analysis, we propose that subfunctionalization after gene duplication plays an important role in the evolution of ligand specificities of the acid-sensing IRs in Lepidoptera.

Leopards are the only big cats still widely distributed across the continents of Africa and Asia. They occur in a wide range of habitats and are often found in close proximity to humans. But despite their ubiquity, leopard phylogeography and population history have not yet been studied with genomic tools. Here, we present population-genomic data from 26 modern and historical samples encompassing the vast geographical distribution of this species. We find that Asian leopards are broadly monophyletic with respect to African leopards across almost their entire nuclear genomes. This profound genetic pattern persists despite the animals' high potential mobility, and despite evidence of transfer of African alleles into Middle Eastern and Central Asian leopard populations within the last 100,000 years. Our results further suggest that Asian leopards originated from a single out-of-Africa dispersal event 500-600 thousand years ago and are characterized by higher population structuring, stronger isolation by distance, and lower heterozygosity than African leopards. Taxonomic categories do not take into account the variability in depth of divergence among subspecies. The deep divergence between the African subspecies and Asian populations contrasts with the much shallower divergence among putative Asian subspecies. Reconciling genomic variation and taxonomy is likely to be a growing challenge in the genomics era.

Outbreaks of the zoonotic disease tularemia occurred in north-east Bulgaria in the 1960s. Then came 30 years of epidemiological silence until new outbreaks occurred in west Bulgaria in the 1990s. To investigate how bacterial strains of Francisella tularensis causing tularemia in wildlife and humans in the 1960s and the 1990s were related, we explored their genetic diversity. Ten F. tularensis genomes from the 1960s (n=3) and the 1990s (n=7) were sequenced, assigned to canonical single-nucleotide polymorphism (canSNP) clades, and compared to reference genomes. We developed four new canSNP polymerase chain reaction (PCR) assays based on the genome sequence information. The genetic analysis showed that the outbreaks in the 1960s as well as in the 1990s involved multiple clones and new genetic diversity. The smallest genetic difference found between any of the Bulgarian strains was five SNPs between the strains L2 and 81 isolated 43 years apart, indicating that F. tularensis may persist locally over long time periods without causing outbreaks. The existence of genetically highly similar strain-pairs isolated the same year in the same area from different hosts supports a hypothesis of local expansion of clones during outbreaks. Close relationship (two SNPs) was found between one strain isolated 1961 in northeast Bulgaria and one strain isolated 5 years before in USSR. Historical data coinciding with the actual time point describe the introduction of water rats from USSR into the Bulgarian outbreak area, which may explain the close genetic relationship and the origin of the outbreak. Genome analysis of strains from two outbreaks in the 1960s and the 1990s provided valuable information on the genetic diversity and persistence of F. tularensis in Bulgaria.

Infections by HIV increase the risk of acquiring secondary viral and bacterial infections and methods are needed to determine the spectrum of co-infections for proper treatment. We used rolling circle amplification (RCA) and Ion Proton sequencing to investigate the vaginal microbiome of 20 HIV positive women from South Africa. A total of 46 different human papillomavirus (HPV) types were found, many of which are not detected by existing genotyping assays. Moreover, the complete genomes of two novel HPV types were determined. Abundance of HPV infections was highly correlated with real-time PCR estimates, indicating that the RCA-Proton method can be used for quantification of individual pathogens. We also identified a large number of other viral, bacterial and parasitic co-infections and the spectrum of these co-infections varied widely between individuals. Our method provides rapid detection of a broad range of pathogens and the ability to reconstruct complete genomes of novel infectious agents.

Methylmercury is a potent human neurotoxin which biomagnifies in aquatic food webs. Although anaerobic microorganisms containing the hgcA gene potentially mediate the formation of methylmercury in natural environments, the diversity of these mercury-methylating microbial communities remains largely unexplored. Previous studies have implicated sulfate-reducing bacteria as the main mercury methylators in aquatic ecosystems. In the present study, we characterized the diversity of mercury-methylating microbial communities of boreal lake sediments using high-throughput sequencing of 16S rRNA and hgcA genes. Our results show that in the lake sediments, Methanomicrobiales and Geobacteraceae also represent abundant members of the mercury-methylating communities. In fact, incubation experiments with a mercury isotopic tracer and molybdate revealed that only between 38% and 45% of mercury methylation was attributed to sulfate reduction. These results suggest that methanogens and iron-reducing bacteria may contribute to more than half of the mercury methylation in boreal lakes.IMPORTANCE Despite the global awareness that mercury, and methylmercury in particular, is a neurotoxin to which millions of people continue to be exposed, there are sizable gaps in the understanding of the processes and organisms involved in methylmercury formation in aquatic ecosystems. In the present study, we shed light on the diversity of the microorganisms responsible for methylmercury formation in boreal lake sediments. All the microorganisms identified are associated with the processing of organic matter in aquatic systems. Moreover, our results show that the well-known mercury-methylating sulfate-reducing bacteria constituted only a minor portion of the potential mercury methylators. In contrast, methanogens and iron-reducing bacteria were important contributors to methylmercury formation, highlighting their role in mercury cycling in the environment.

Chemotherapy-induced myelosuppression, including thrombocytopenia, is a recurrent problem during cancer treatments that may require dose alterations or cessations that could affect the antitumor effect of the treatment. To identify genetic markers associated with treatment-induced thrombocytopenia, we whole-exome sequenced 215 non-small cell lung cancer patients homogeneously treated with gemcitabine/carboplatin. The decrease in platelets (defined as nadir/baseline) was used to assess treatment-induced thrombocytopenia. Association between germline genetic variants and thrombocytopenia was analyzed at single-nucleotide variant (SNV) (based on the optimal false discovery rate, the severity of predicted consequence, and effect), gene, and pathway levels. These analyses identified 130 SNVs/INDELs and 25 genes associated with thrombocytopenia (P-value < 0.002). Twenty-three SNVs were validated in an independent genome-wide association study (GWAS). The top associations include rs34491125 in JMJD1C (P-value = 9.07 × 10-5), the validated variants rs10491684 in DOCK8 (P-value = 1.95 × 10-4), rs6118 in SERPINA5 (P-value = 5.83 × 10-4), and rs5877 in SERPINC1 (P-value = 1.07 × 10-3), and the genes CAPZA2 (P-value = 4.03 × 10-4) and SERPINC1 (P-value = 1.55 × 10-3). The SNVs in the top-scoring pathway "Factors involved in megakaryocyte development and platelet production" (P-value = 3.34 × 10-4) were used to construct weighted genetic risk score (wGRS) and logistic regression models that predict thrombocytopenia. The wGRS predict which patients are at high or low toxicity risk levels, for CTCAE (odds ratio (OR) = 22.35, P-value = 1.55 × 10-8), and decrease (OR = 66.82, P-value = 5.92 × 10-9). The logistic regression models predict CTCAE grades 3-4 (receiver operator characteristics (ROC) area under the curve (AUC) = 0.79), and large decrease (ROC AUC = 0.86). We identified and validated genetic variations within hematopoiesis-related pathways that provide a solid foundation for future studies using genetic markers for predicting chemotherapy-induced thrombocytopenia and personalizing treatments.

A meta-analysis of genome-wide association studies (GWAS) identified eight loci that are associated with heart rate variability (HRV), but candidate genes in these loci remain uncharacterized. We developed an image- and CRISPR/Cas9-based pipeline to systematically characterize candidate genes for HRV in live zebrafish embryos. Nine zebrafish orthologues of six human candidate genes were targeted simultaneously in eggs from fish that transgenically express GFP on smooth muscle cells (Tg[acta2:GFP]), to visualize the beating heart. An automated analysis of repeated 30 s recordings of beating atria in 381 live, intact zebrafish embryos at 2 and 5 days post-fertilization highlighted genes that influence HRV (hcn4 and si:dkey-65j6.2 [KIAA1755]); heart rate (rgs6 and hcn4); and the risk of sinoatrial pauses and arrests (hcn4). Exposure to 10 or 25 µM ivabradine-an open channel blocker of HCNs-for 24 h resulted in a dose-dependent higher HRV and lower heart rate at 5 days post-fertilization. Hence, our screen confirmed the role of established genes for heart rate and rhythm (RGS6 and HCN4); showed that ivabradine reduces heart rate and increases HRV in zebrafish embryos, as it does in humans; and highlighted a novel gene that plays a role in HRV (KIAA1755).