Search

The Y chromosome harbors nine multi-copy ampliconic gene families expressed exclusively in testis. The gene copies within each family are >99% identical to each other, which poses a major challenge in evaluating their copy number. Recent studies demonstrated high variation in Y ampliconic gene copy number among humans. However, how this variation affects expression levels in human testis remains understudied. Here we developed a novel computational tool Ampliconic Copy Number Estimator (AmpliCoNE) that utilizes read sequencing depth information to estimate Y ampliconic gene copy number per family. We applied this tool to whole-genome sequencing data of 149 men with matched testis expression data whose samples are part of the Genotype-Tissue Expression (GTEx) project. We found that the Y ampliconic gene families with low copy number in humans were deleted or pseudogenized in non-human great apes, suggesting relaxation of functional constraints. Among the Y ampliconic gene families, higher copy number leads to higher expression. Within the Y ampliconic gene families, copy number does not influence gene expression, rather a high tolerance for variation in gene expression was observed in testis of presumably healthy men. No differences in gene expression levels were found among major Y haplogroups. Age positively correlated with expression levels of the HSFY and PRY gene families in the African subhaplogroup E1b, but not in the European subhaplogroups R1b and I1. We also found that expression of five out of six Y ampliconic gene families is coordinated with that of their non-Y (i.e. X or autosomal) homologs. Indeed, five ampliconic gene families had consistently lower expression levels when compared to their non-Y homologs suggesting dosage regulation, while the HSFY family had higher expression levels than its X homolog and thus lacked dosage regulation. [ABSTRACT FROM AUTHOR]

Interspecific hybridization can introduce genetic variation that aids in adaptation to new or changing environments. Here, we investigate how hybrid adaptation to temperature and nutrient limitation may alter parental genome representation over time. We evolved Saccharomyces cerevisiae x Saccharomyces uvarum hybrids in nutrient-limited continuous culture at 15°C for 200 generations. In comparison to previous evolution experiments at 30°C, we identified a number of responses only observed in the colder temperature regime, including the loss of the S. cerevisiae allele in favor of the cryotolerant S. uvarum allele for several portions of the hybrid genome. In particular, we discovered a genotype by environment interaction in the form of a loss of heterozygosity event on chromosome XIII; which species’ haplotype is lost or maintained is dependent on the parental species’ temperature preference and the temperature at which the hybrid was evolved. We show that a large contribution to this directionality is due to a temperature dependent fitness benefit at a single locus, the high affinity phosphate transporter gene PHO84. This work helps shape our understanding of what forces impact genome evolution after hybridization, and how environmental conditions may promote or disfavor the persistence of hybrids over time. [ABSTRACT FROM AUTHOR]

Highly pathogenic avian influenza (HPAI) virus subtype H7N3 has been circulating in poultry in Mexico since 2012 and vaccination has been used to control the disease. In this study, eight Mexican H7N3 HPAI viruses from 2015–2017 were isolated and fully sequenced. No evidence of reassortment was detected with other avian influenza (AI) viruses, but phylogenetic analyses show divergence of all eight gene segments into three genetic clusters by 2015, with 94.94 to 98.78 percent nucleotide homology of the HA genes when compared to the index virus from 2012. The HA protein of viruses from each cluster showed a different number of basic amino acids (n = 5–7) in the cleavage site, and six different patterns at the predicted N-glycosylation sites. Comparison of the sequences of the Mexican lineage H7N3 HPAI viruses and American ancestral wild bird AI viruses to characterize the virus evolutionary dynamics showed that the nucleotide substitution rates in PB2, PB1, PA, HA, NP, and NS genes greatly increased once the virus was introduced into poultry. The global nonsynonymous and synonymous ratios imply strong purifying selection driving the evolution of the virus. Forty-nine positively selected sites out of 171 nonsynonymous mutations were identified in the Mexican H7N3 HPAI viruses, including 7 amino acid changes observed in higher proportion in North American poultry origin AI viruses isolates than in wild bird-origin viruses. Continuous monitoring and molecular characterization of the H7N3 HPAI virus is important for better understanding of the virus evolutionary dynamics and further improving control measures including vaccination. [ABSTRACT FROM AUTHOR]

Background: More than 11,000 laboratories and companies developed their own next-generation sequencing (NGS) for screening and diagnosis of various diseases including cancer. Although inconsistencies of mutation calls as high as 43% in databases such as GDSC (Genomics of Drug Sensitivity in Cancer) and CCLE (Cancer Cell Line Encyclopedia) have been reported, not many studies on the reasons for the inconsistencies have been published. Methods: Targeted-NGS analysis of 151 genes in 35 cell lines common to GDSC and CCLE was performed, and the results were compared with those from GDSC and CCLE wherein whole-exome- or highly-multiplex NGS were employed. Results: In the comparison, GDSC and CCLE had a high rate (40–45%) of false-negative (FN) errors which would lead to high rate of inconsistent mutation calls, suggesting that highly-multiplex NGS may have high rate of FN errors. We also posited the possibility that targeted NGS, especially for the detection of low-level cancer cells in cancer tissues might suffer significant FN errors. Conclusion: FN errors may be the most important errors in NGS testing for cancer; their evaluation in laboratory-developed NGS tests is needed. [ABSTRACT FROM AUTHOR]

The South American plains vizcacha, Lagostomus maximus inhabits primarily the Pampean and adjoining Espinal, Monte and Chaquenean regions of Argentina. In order to study the population genetic structure of L. maximus, a fragment of 560 bp of the mitochondrial DNA hypervariable region 1from 90 individuals collected from the 3 subspecies and 8 groups along Argentina was amplified and analyzed. We found 9 haplotypes. The haplotype network did not show an apparent phylogeographical signal. Although low levels of genetic variation were found in all the subspecies and groups analyzed, a radiation of L. maximus would have occurred from the North and Center of the Pampean region toward the rest of its geographic range in Argentina. Low levels of genetic diversity, the existence of a single genetically distinct population in Argentina and changes of its effective size indicate that metapopulation processes and changes in human population dynamics during the late-Holocene were important factors shaping the population genetic structure of L. maximus in Argentina. [ABSTRACT FROM AUTHOR]

Most current methods for detecting natural selection from DNA sequence data are limited in that they are either based on summary statistics or a composite likelihood, and as a consequence, do not make full use of the information available in DNA sequence data. We here present a new importance sampling approach for approximating the full likelihood function for the selection coefficient. Our method treats the ancestral recombination graph (ARG) as a latent variable that is integrated out using previously published Markov Chain Monte Carlo (MCMC) methods. The method can be used for detecting selection, estimating selection coefficients, testing models of changes in the strength of selection, estimating the time of the start of a selective sweep, and for inferring the allele frequency trajectory of a selected or neutral allele. We perform extensive simulations to evaluate the method and show that it uniformly improves power to detect selection compared to current popular methods such as nSL and SDS, and can provide reliable inferences of allele frequency trajectories under many conditions. We also explore the potential of our method to detect extremely recent changes in the strength of selection. We use the method to infer the past allele frequency trajectory for a lactase persistence SNP (MCM6) in Europeans. We also infer the trajectory of a SNP (EDAR) in Han Chinese, finding evidence that this allele’s age is much older than previously claimed. We also study a set of 11 pigmentation-associated variants. Several genes show evidence of strong selection particularly within the last 5,000 years, including ASIP, KITLG, and TYR. However, selection on OCA2/HERC2 seems to be much older and, in contrast to previous claims, we find no evidence of selection on TYRP1. [ABSTRACT FROM AUTHOR]

We investigated the phylogeography of the dugong (Dugong dugon) across its original range using museum material from 14 natural history museum and university collections. The mitochondrial DNA control region was successfully amplified from samples of bone or tooth powder from 162 individuals. These samples range from 1827 to 1996 and span the historical distribution range of the dugong. We were able to successfully amplify overlapping fragments of the D-loop region of the mitochondrial DNA (mtDNA) resulting in sequences of a 355 bp fragment for 162 individuals for the final analyses. This included a new sequence (189 bp) from a previously unidentified piece of skin of the extinct Steller’s sea cow (Hydrodamalis gigas), as an outgroup. The resulting dugong sequences match those from previous studies of dugongs from Australia and Indonesia, but revealed several new and divergent mtDNA lineages in the Indian Ocean. One mtDNA lineage includes most specimens from the Western Indian Ocean, with another distinct lineage isolated to nearby Madagascar and Comores. There is little geographic structuring detectable among other populations in the Western Indian Ocean and all populations from that region appear to have historically contained comparatively low levels of genetic diversity. The genetic diversity of several Indian Ocean samples collected after 1950 was lower than that of the samples collected earlier from similar locations, a result coincident with the anecdotal reductions in population size. The new lineages and potential loss of diversity highlight the particular conservation importance and vulnerability of dugong populations in the Western Indian Ocean. [ABSTRACT FROM AUTHOR]

Insects possess highly developed olfactory systems which play pivotal roles in its ecological adaptations, host plant location, and oviposition behavior. Bactrocera minax is an oligophagous tephritid insect whose host selection, and oviposition behavior largely depend on the perception of chemical cues. However, there have been very few reports on molecular components related to the olfactory system of B. minax. Therefore, the transcriptome of B. minax were sequenced in this study, with 1 candidate chemosensory protein (CSP), 21 candidate odorant binding proteins (OBPs), 53 candidate odorant receptors (ORs), 29 candidate ionotropic receptors (IRs) and 4 candidate sensory neuron membrane proteins (SNMPs) being identified. After that, we sequenced the candidate olfactory genes and performed phylogenetic analysis. qRT-PCR was used to express and characterize 9 genes in olfactory and non-olfactory tissues. Compared with GFP-injected fly (control), dsOBP21-treated B. minax and dsCSP-treated B. minax had lower electrophysiological response to D-limonene (attractant), suggesting the potential involvement of BminOBP21 and BminCSP genes in olfactory perceptions of the fly. Our study establishes the molecular basis of olfaction, tributary for further functional analyses of chemosensory processes in B. minax. [ABSTRACT FROM AUTHOR]

A comprehensive, accurate, and revisable alpha taxonomy is crucial for biodiversity studies, but is challenging when data from reference specimens are difficult to collect or observe. However, recent technological advances can overcome some of these challenges. To illustrate this, we used modern approaches to tackle a centuries-old taxonomic enigma presented by Fraser’s Clawed Frog, Xenopus fraseri, including whether X. fraseri is different from other species, and if so, where it is situated geographically and phylogenetically. To facilitate these inferences, we used high-resolution techniques to examine morphological variation, and we generated and analyzed complete mitochondrial genome sequences from all Xenopus species, including >150-year-old type specimens. Our results demonstrate that X. fraseri is indeed distinct from other species, firmly place this species within a phylogenetic context, and identify its minimal geographic distribution in northern Ghana and northern Cameroon. These data also permit novel phylogenetic resolution into this intensively studied and biomedically important group. Xenopus fraseri was formerly thought to be a rainforest endemic placed alongside species in the amieti species group; in fact this species occurs in arid habitat on the borderlands of the Sahel, and is the smallest member of the muelleri species group. This study illustrates that the taxonomic enigma of Fraser’s frog was a combined consequence of sparse collection records, interspecies conservation and intraspecific polymorphism in external anatomy, and type specimens with unusual morphology. [ABSTRACT FROM AUTHOR]

In Argentina, NDM metallo-β-lactamase was first reported in 2013. By now, it has disseminated throughout the country in diverse Gram negative bacteria. Here, we report the case of a paediatric patient that underwent a 1-year hospitalisation due to erythrodermic psoriasis in 2014 and received multiple antimicrobial treatments. During his stay, five isolates were obtained from rectal swabs (rs) or blood culture (bc) suspicious of carbapenemase production: a K. quasipneumoniae subsp. quasipneumoniae (rs), Citrobacter freundii (rs), Escherichia coli (bc), Enterobacter cloacae (rs), and a Serratia marcescens (bc). The isolates were studied with broth microdilution, biparental conjugation and plasmid and whole genome sequencing (Illumina). All isolates harboured an 138,998-bp type 1 IncC plasmid that carried blaNDM-1, bleMBL, blaCMY-6, rmtC, aac(6’)-Ib, and sul1 resistance genes. Additionally, the blaNDM-plasmids contained ISKpn8 an insertion sequence previously described as associated only to blaKPC. One isolate, a colistin-resistant E. coli, also carried a mcr-1-containing an IncI2 plasmid, which did not harbour additional resistance. The whole genome of K. quasipneumoniae subsp. quasipneumoniae isolate was fully sequenced. This isolate harboured, additionally to blaNDM, three plasmid-mediated quinolone resistance genes: qnrB4, qnrB52 and aac(6’)-Ib-cr1. The E. cloacae isolate also harboured qnrA1. These findings alert to the underestimated horizontal dissemination of multidrug-resistant plasmids limiting treatment options with last resort antimicrobials. [ABSTRACT FROM AUTHOR]

The microbiome plays a key role in the biology, ecology and evolution of arthropod vectors of human pathogens. Vector-bacterial interactions could alter disease transmission dynamics through modulating pathogen replication and/or vector fitness. Nonetheless, our understanding of the factors shaping the bacterial community in arthropod vectors is incomplete. Using large-scale 16S amplicon sequencing, we examine how habitat disturbance structures the bacterial assemblages of field-collected whole-body hematophagous arthropods that vector human pathogens including mosquitoes (Culicidae), sand flies (Psychodidae), biting midges (Ceratopogonidae) and hard ticks (Ixodidae). We found that all comparisons of the bacterial community among species yielded statistically significant differences, but a difference was not observed between adults and nymphs of the hard tick, Haemaphysalis juxtakochi. While Culicoides species had the most distinct bacterial community among dipterans, tick species were composed of entirely different bacterial OTU’s. We observed differences in the proportions of some bacterial types between pristine and disturbed habitats for Coquillettidia mosquitoes, Culex mosquitoes, and Lutzomyia sand flies, but their associations differed within and among arthropod assemblages. In contrast, habitat quality was a poor predictor of differences in bacterial classes for Culicoides biting midges and hard tick species. In general, similarities in the bacterial communities among hematophagous arthropods could be explained by their phylogenetic relatedness, although intraspecific variation seems influenced by habitat disturbance. [ABSTRACT FROM AUTHOR]

Background: Pregnane X receptor (PXR) regulates the expression of drug-metabolizing enzymes and transport enzymes. NF-κB not only plays a role in liver homeostasis and injury-healing processes by regulating inflammatory responses but may also regulate the transcription of PXR. Currently, genetic polymorphisms in PXR are associated with adverse drug effects. Because little is known about the association between NF-κB1 genetic polymorphisms and adverse drug reactions, we explored the association between PXR and NF-κB1 single nucleotide polymorphisms (SNPs) and susceptibility to anti-tuberculosis drug-induced liver injury (ATDILI). Materials and methods: A total of 746 tuberculosis patients (118 with ATDILI and 628 without ATDILI) were prospectively enrolled at West China Hospital between December 2014 and April 2018. Nine selected SNPs (rs3814055, rs13059232, rs7643645 and rs3732360 in PXR and rs78872571, rs4647992, rs60371688, rs1598861 and rs3774959 in NF-κB1) were genotyped with a custom-designed 2x48-plex SNP Scan TM Kit. The frequencies of the alleles, genotypes and genetic models of the variants were compared between patients with or without ATDILI, while joint effect analysis of the SNP-SNP interactions was performed using multiplicative and additive models. The odds ratios (ORs) and the corresponding 95% confidence intervals (CIs) were calculated. Results: The T allele of rs3814055 in PXR was associated with a decreased risk for ATDILI (OR 0.61; 95% CI: 0.42–0.89, p = 0.0098). The T alleles of rs78872571 and rs4647992 in NF-κB1 were significantly associated with an increased risk for ATDILI (OR 1.91; 95% CI: 1.06–3.43, p = 0.028 and OR 1.81; 1.06–3.10, p = 0.029, respectively). The allele, genotype and genetic model frequencies were similar in the two groups for the other six SNPs (all P>0.05). There were no multiplicative or additive interactions between the SNPs. Conclusion: Our study is the first to reveal that rs3814055 variants in PXR and rs78872571 and rs4647992 variants in NF-κB1 are associated with susceptibility to ATDILI caused by first-line anti-tuberculosis combination treatment in the Han Chinese population. [ABSTRACT FROM AUTHOR]

Despite the importance of microbial activity in mobilizing arsenic in groundwater aquifers, the phylogenetic distribution of contributing microbial metabolisms is understudied. Groundwater samples from Ohio aquifers were analyzed using metagenomic sequencing to identify functional potential that could drive arsenic cycling, and revealed mechanisms for direct (i.e., Ars system) and indirect (i.e., iron reduction) arsenic mobilization in all samples, despite differing geochemical conditions. Analyses of 194 metagenome-assembled genomes (MAGs) revealed widespread functionality related to arsenic mobilization throughout the bacterial tree of life. While arsB and arsC genes (components of an arsenic resistance system) were found in diverse lineages with no apparent phylogenetic bias, putative aioA genes (aerobic arsenite oxidase) were predominantly identified in Methylocystaceae MAGs. Both previously described and undescribed respiratory arsenate reduction potential via arrA was detected in Betaproteobacteria, Deltaproteobacteria, and Nitrospirae MAGs, whereas sulfate reduction potential was primarily limited to members of the Deltaproteobacteria and Nitrospirae. Lastly, iron reduction potential was detected in the Ignavibacteria, Deltaproteobacteria, and Nitrospirae. These results expand the phylogenetic distribution of taxa that may play roles in arsenic mobilization in subsurface systems. Specifically, the Nitrospirae are a much more functionally diverse group than previously assumed and may play key biogeochemical roles in arsenic-contaminated ecosystems. [ABSTRACT FROM AUTHOR]

Sexual development is a key evolutionary innovation of eukaryotes. In many species, mating involves interaction between compatible mating partners that can undergo cell and nuclear fusion and subsequent steps of development including meiosis. Mating compatibility in fungi is governed by the mating type (MAT) loci. In basidiomycetes, the ancestral state is hypothesized to be tetrapolar, with two genetically unlinked MAT loci containing homeodomain transcription factor genes (HD locus) and pheromone and pheromone receptor genes (P/R locus), respectively. Alleles at both loci must differ between mating partners for completion of sexual development. However, there are also basidiomycetes with bipolar mating systems, which can arise through genomic linkage of the HD and P/R loci. In the order Tremellales, bipolarity is found only in the pathogenic Cryptococcus species. Here, we describe the analysis of MAT loci from 24 species of the Trichosporonales, a sister order to the Tremellales. In all of the species analyzed, the MAT loci are fused and a single HD gene is present in each mating type, similar to the organization in the pathogenic Cryptococci. However, the HD and P/R allele combinations in the Trichosporonales are different from those in the pathogenic Cryptococci. This and the existence of tetrapolar species in the Tremellales suggest that fusion of the HD and P/R loci occurred independently in the Trichosporonales and pathogenic Cryptococci, supporting the hypothesis of convergent evolution towards fused MAT regions, similar to previous findings in other fungal groups. Unlike the fused MAT loci in several other basidiomycete lineages though, the gene content and gene order within the fused MAT loci are highly conserved in the Trichosporonales, and there is no apparent suppression of recombination extending from the MAT loci to adjacent chromosomal regions, suggesting different mechanisms for the evolution of physically linked MAT loci in these groups. [ABSTRACT FROM AUTHOR]

Background: Burkholderia pseudomallei is a soil-dwelling bacterium and the causative agent of melioidosis. The global burden and distribution of melioidosis is poorly understood, including in the Caribbean. B. pseudomallei was previously isolated from humans and soil in eastern Puerto Rico but the abundance and distribution of B. pseudomallei in Puerto Rico as a whole has not been thoroughly investigated. Methodology/Principal findings: We collected 600 environmental samples (500 soil and 100 water) from 60 sites around Puerto Rico. We identified B. pseudomallei by isolating it via culturing and/or using PCR to detect its DNA within complex DNA extracts. Only three adjacent soil samples from one site were positive for B. pseudomallei with PCR; we obtained 55 isolates from two of these samples. The 55 B. pseudomallei isolates exhibited fine-scale variation in the core genome and contained four novel genomic islands. Phylogenetic analyses grouped Puerto Rico B. pseudomallei isolates into a monophyletic clade containing other Caribbean isolates, which was nested inside a larger clade containing all isolates from Central/South America. Other Burkholderia species were commonly observed in Puerto Rico; we cultured 129 isolates from multiple soil and water samples collected at numerous sites around Puerto Rico, including representatives of B. anthina, B. cenocepacia, B. cepacia, B. contaminans, B. glumae, B. seminalis, B. stagnalis, B. ubonensis, and several unidentified novel Burkholderia spp. Conclusions/Significance: B. pseudomallei was only detected in three soil samples collected at one site in north central Puerto Rico with only two of those samples yielding isolates. All previous human and environmental B. pseudomallei isolates were obtained from eastern Puerto Rico. These findings suggest B. pseudomallei is ecologically established and widely dispersed in the environment in Puerto Rico but rare. Phylogeographic patterns suggest the source of B. pseudomallei populations in Puerto Rico and elsewhere in the Caribbean may have been Central or South America. [ABSTRACT FROM AUTHOR]

The interaction between phenotypic plasticity, e.g. learning, and evolution is an important topic both in Evolutionary Biology and Machine Learning. The evolution of learning is commonly studied in Evolutionary Biology, while the use of an evolutionary process to improve learning is of interest to the field of Machine Learning. This paper takes a different point of view by studying the effect of learning on the evolutionary process, the so-called Baldwin effect. A well-studied result in the literature about the Baldwin effect is that learning affects the speed of convergence of the evolutionary process towards some genetic configuration, which corresponds to the environment-induced plastic response. This paper demonstrates that learning can change the outcome of evolution, i.e., lead to a genetic configuration that does not correspond to the plastic response. Results are obtained both analytically and experimentally by means of an agent-based model of a foraging task, in an environment where the distribution of resources follows seasonal cycles and the foraging success on different resource types is conditioned by trade-offs that can be evolved and learned. This paper attempts to answer a question that has been overlooked: whether learning has an effect on what genotypic traits are evolved, i.e. the selection of a trait that enables a plastic response changes the selection pressure on a different trait, in what could be described as co-evolution between different traits in the same genome. [ABSTRACT FROM AUTHOR]

Bovine tuberculosis (BTB) is endemic to the African buffalo (Syncerus caffer) of Hluhluwe-iMfolozi Park (HiP) and Kruger National Park, South Africa. In HiP, the disease has been actively managed since 1999 through a test-and-cull procedure targeting BTB-positive buffalo. Prior studies in Kruger showed associations between microsatellite alleles, BTB and body condition. A sex chromosomal meiotic drive, a form of natural gene drive, was hypothesized to be ultimately responsible. These associations indicate high-frequency occurrence of two types of male-deleterious alleles (or multiple-allele haplotypes). One type negatively affects body condition and BTB resistance in both sexes. The other type has sexually antagonistic effects: negative in males but positive in females. Here, we investigate whether a similar gene drive system is present in HiP buffalo, using 17 autosomal microsatellites and microsatellite-derived Y-chromosomal haplotypes from 401 individuals, culled in 2002–2004. We show that the association between autosomal microsatellite alleles and BTB susceptibility detected in Kruger, is also present in HiP. Further, Y-haplotype frequency dynamics indicated that a sex chromosomal meiotic drive also occurred in HiP. BTB was associated with negative selection of male-deleterious alleles in HiP, unlike positive selection in Kruger. Birth sex ratios were female-biased. We attribute negative selection and female-biased sex ratios in HiP to the absence of a Y-chromosomal sex-ratio distorter. This distorter has been hypothesized to contribute to positive selection of male-deleterious alleles and male-biased birth sex ratios in Kruger. As previously shown in Kruger, microsatellite alleles were only associated with male-deleterious effects in individuals born after wet pre-birth years; a phenomenon attributed to epigenetic modification. We identified two additional allele types: male-specific deleterious and beneficial alleles, with no discernible effect on females. Finally, we discuss how our findings may be used for breeding disease-free buffalo and implementing BTB test-and-cull programs. [ABSTRACT FROM AUTHOR]

Whole mitogenome sequences (mtDNA) have been exploited for insect ecology studies, using them as molecular markers to reconstruct phylogenies, or to infer phylogeographic relationships and gene flow. Recent Anopheles phylogenomic studies have provided information regarding the time of deep lineage divergences within the genus. Here we report the complete 15,393 bp mtDNA sequences of Anopheles aquasalis, a Neotropical human malaria vector. When comparing its structure and base composition with other relevant and available anopheline mitogenomes, high similarity and conserved genomic features were observed. Furthermore, 22 mtDNA sequences comprising anopheline and Dipteran sibling species were analyzed to reconstruct phylogenies and estimate dates of divergence between taxa. Phylogenetic analysis using complete mtDNA sequences suggests that A. aquasalis diverged from the Anopheles albitarsis complex ~28 million years ago (MYA), and ~38 MYA from Anopheles darlingi. Bayesian analysis suggests that the most recent ancestor of Nyssorhynchus and Anopheles + Cellia was extant ~83 MYA, corroborating current estimates of ~79–100 MYA. Additional sampling and publication of African, Asian, and North American anopheline mitogenomes would improve the resolution of the Anopheles phylogeny and clarify early continental dispersal routes. [ABSTRACT FROM AUTHOR]

The biological significance of regional cladoceran morphotypes in the montane regions of the central Palearctic remains poorly understood. In the Holarctic Daphnia longispina complex (Cladocera: Daphniidae), several variants, lineages and species have been proposed as endemic for Southern Siberia. Daphnia turbinata Sars, for example, named after its unusual head shape, is known only from Southern Siberia. Here we sequence DNA of Daphnia from three mitochondrial genes (12S rRNA, 16S rRNA, and NADH dehydrogenase subunit 2, ND2) from 57 localities in Russia and Mongolia (the majority being from Southern Siberia) and place them in evolutionary context with existing data. Our aim was to examine regional endemism of the Daphnia longispina complex in Southern Siberian; to improve the phylogenetic understanding with improved taxonomic and regional sampling, and to better understand the influence of Pleistocene glaciation on the biogeography of these lineages. At least three lineages showed genetic evidence for endemism in Southern Siberia. There was strong support for D. turbinata as a sister lineage to to D. longispina/D. dentifera. Another endemic, Siberian D. cf. longispina, is a sister group to the longispina group in general. Within D. longispina s. str. there was an endemic Siberian clade with a western range boundary near the Yenisei River Basin. Gene flow estimates among populations (based on FST values) were very low for clades of D. longispina on a regional (the original 12S dataset), and on a pan-Eurasian (the extended 12S dataset) scale. Negative values of Fu’s FS and Tajima’s D tests prevailed for the species examined with significant values found for two D. longispina clades, D. dentifera, D. galeata and D. cristata. Our results support the notion that Southern Siberia is an important biogeographic region for cladocerans as it contained unexpected diversity of endemics (such as D. turbinata, D. cf. longispina and lineages of D. umbra and D. longsipina s.str.) and from being the geographic meeting place of expanding postglacial lineages from eastern and western refugia. [ABSTRACT FROM AUTHOR]

Leishmania (Leishmania) infantum is the zoonotic agent of visceral leishmaniasis (VL), a disease with a global distribution. The transmission scenario of VL has been undergoing changes worldwide, with the biologic cycle invading urbanized areas and dispersing the parasites into other previously free areas. The epidemiological cycle in Brazil has dispersed from the Northeast to other regions of the country. In this study, an integrative approach, including genotyping Brazilian strains of L. (L.) infantum for 14 microsatellite markers and reviewing historical records of the disease, was used to assess dispersion routes throughout central-southern Brazil. Our results support three L. (L.) infantum dispersion routes: A) dispersion from Bolivia to the states of Mato Grosso, Mato Grosso do Sul and São Paulo via the Bolivia-Brazil gas pipeline from 1998 to 2005; B) VL dispersion from Paraguay to the Brazilian side of the triple border (Foz do Iguaçu and Santa Terezinha de Itaipu) during after 2012; and C) emergence of a new L. (L.) infantum cluster in western Santa Catarina State and its dispersion to southern Paraná State (municipality of Pato Branco), after 2013. Hypotheses regarding possible entries of Leishmania (L.) infantum into the area of the triple border are presented and discussed. Understanding how VL has dispersed is vital to the development of control measures for this disease and to avoid future dispersion events. [ABSTRACT FROM AUTHOR]

Cooperation among kin is common in animal societies. Kin groups may form by individuals directly discriminating relatives based on kin recognition cues, or form passively through natal philopatry and limited dispersal. We describe the genetic landscape for a primitively eusocial wasp, Polistes dominula, and ask whether individuals choose cooperative partners that are nearby and/or that are genetic relatives. Firstly, we genotyped an entire sub-population of 1361 wasps and found genetic structuring on an extremely fine scale: the probability of finding genetic relatives decreases exponentially within just a few meters of an individual’s nest. At the same time, however, we found a lack of genetic structuring between natural nest aggregations within the population. Secondly, in a separate dataset where ~2000 wasps were genotyped, we show that wasps forced experimentally to make a new nest choice tended to choose new nests near to their original nests, and that these nests tended to contain some full sisters. However, a significant fraction of wasps chose nests that did not contain sisters, despite sisters being present in nearby nests. Although we cannot rule out a role for direct kin recognition or natal nest-mate recognition, our data suggest that kin groups may form via a philopatric rule-of-thumb, whereby wasps simply select groups and nesting sites that are nearby. The result is that most subordinate helpers obtain indirect fitness benefits by breeding cooperatively. [ABSTRACT FROM AUTHOR]

The heterotrophic flagellate Percolomonas cosmopolitus (Heterolobosea) is often observed in saline habitats worldwide, from coastal waters to saturated brines. However, only two cultures assigned to this morphospecies have been examined using molecular methods, and their 18S rRNA gene sequences are extremely different. Further the salinity tolerances of individual strains are unknown. Thus, our knowledge on the autecology and diversity in this morphospecies is deficient. Here, we report 18S rRNA gene data on seven strains similar to P. cosmopolitus from seven geographically remote locations (New Zealand, Kenya, Korea, Poland, Russia, Spain, and the USA) with sample salinities ranging from 4‰ to 280‰, and compare morphology and salinity tolerance of the nine available strains. Percolomonas cosmopolitus-like strains show few-to-no consistent morphological differences, and form six clades separated by often extremely large 18S rRNA gene divergences (up to 42.4%). Some strains grow best at salinities from 75 to 125‰ and represent halophiles. All but one of these belong to two geographically heterogeneous clusters that form a robust monophyletic group in phylogenetic trees; this likely represents an ecologically specialized subclade of halophiles. Our results suggest that P. cosmopolitus is a cluster of several cryptic species (at least), which are unlikely to be distinguished by geography. Interestingly, the 9 Percolomonas strains formed a clade in 18S rRNA gene phylogenies, unlike most previous analyses based on two sequences. [ABSTRACT FROM AUTHOR]

In response to strong selection, new mutations can arise quickly and sweep through populations, particularly, if survival and reproduction depend on certain allele copies for adaptation to rapidly changing environments, like resistance against deadly diseases or strong toxins. Since the 1950s, resistance to anticoagulant rodenticides in several rodents has emerged through single nucleotide mutations in the vitamin-K-epoxid-reductase-complex-subunit-1 (VKORC1) gene, often located in its exon 3. Detection of high prevalence and concentrations of anticoagulant rodenticides in non-target vertebrates, including carnivorous Mustelidae, let us assume that secondary exposure by feeding on poisoned prey may also cause selection along the food chain and we hypothesized that VKORC1-based resistance might also have evolved in rodents’ predators. Using newly-developed mustelid-specific primers for direct sequencing of genomic DNA, we studied VKORC1-DNA-polymorphisms in 115 mustelids of five species (Martes martes, M. foina, Mustela nivalis, M. erminea, M. putorius), obtained from northern Denmark, yielding six sites with nonsynonymous and several synonymous amino acid polymorphisms in exon 3. Comparison of these VKORC1-genotypes with hepatic rodenticide residues (obtained by HPLC combined with fluorescence or mass spectrometry) in 83 individuals (except M. martes), using generalized linear models, suggested that anticoagulant levels depended on species and specific polymorphisms. Although most VKORC-1 polymorphisms may present standing genetic variation, some are situated in resistance-mediating membrane parts of the VKORC1-encoded protein, and might be a result of selection due to exposure to anticoagulant poisons. Our new molecular markers might allow detecting indirect effects of anticoagulant rodenticides on rodent predator populations in the future. [ABSTRACT FROM AUTHOR]

Bengal tigers (Panthera tigris tigris) serve a pivotal role as an apex predator in forest ecosystems. To increase our knowledge on factors impacting the viability and health of this endangered species, we studied the gut microbiota in 32 individual Bengal tigers from three geographically separated areas (Chitwan National Park (CNP), Bardia National Park (BNP) and Suklaphanta Wildlife Reserve (SWR)) in Nepal, using noninvasive genetic sampling methods. Gut microbiota influence the immune system, impact various physiological functions, and modulates metabolic reactions, that ultimately impact the host health, behavior and development. Across the tiger populations in Nepal, we found significant differences in the composition of microbial communities based on their geographic locations. Specifically, we detected significant differences between CNP and the other two protected areas (CNP vs BNP: pseudo t = 1.944, P = 0.006; CNP vs SWR: pseudo t = 1.9942, P = 0.0071), but no differences between BNP and SWR. This mirrors what has been found for tiger gene flow in the same populations, suggesting gut microbiota composition and host gene flow may be linked. Furthermore, predictive metagenome functional content analysis (PICRUSt) revealed a higher functional enrichment and diversity for significant gut microbiota in the Chitwan tiger population and the lowest enrichment and diversity in Suklaphanta. The CNP tiger population contained higher proportions of microbiota that are associated with predicted functions relevant for metabolism of amino acid, lipid, xenobiotics biodegradation, terpenoides and polyketides than the SWR population. We conclude the tiger population structure, gut microbiota profile and associated functional metabolic categories are correlated, with geographically most separated CNP and SWR tiger population having the most distinct and different host genotype and microbiota profiles. Our work dramatically expands the understanding of tiger microbiota in wild populations and provides a valuable case study on how to investigate genetic diversity at different hierarchical levels, including hosts as well as their microbial communities. [ABSTRACT FROM AUTHOR]

The microbial communities that inhabit the distal gut of humans and other mammals exhibit large inter-individual variation. While host genetics is a known factor that influences gut microbiota composition, the mechanisms underlying this variation remain largely unknown. Bile acids (BAs) are hormones that are produced by the host and chemically modified by gut bacteria. BAs serve as environmental cues and nutrients to microbes, but they can also have antibacterial effects. We hypothesized that host genetic variation in BA metabolism and homeostasis influence gut microbiota composition. To address this, we used the Diversity Outbred (DO) stock, a population of genetically distinct mice derived from eight founder strains. We characterized the fecal microbiota composition and plasma and cecal BA profiles from 400 DO mice maintained on a high-fat high-sucrose diet for ~22 weeks. Using quantitative trait locus (QTL) analysis, we identified several genomic regions associated with variations in both bacterial and BA profiles. Notably, we found overlapping QTL for Turicibacter sp. and plasma cholic acid, which mapped to a locus containing the gene for the ileal bile acid transporter, Slc10a2. Mediation analysis and subsequent follow-up validation experiments suggest that differences in Slc10a2 gene expression associated with the different strains influences levels of both traits and revealed novel interactions between Turicibacter and BAs. This work illustrates how systems genetics can be utilized to generate testable hypotheses and provide insight into host-microbe interactions. [ABSTRACT FROM AUTHOR]

The production of anticompetitor toxins is widespread among bacteria. Because production of such toxins is costly, it is typically regulated. In particular, many toxins are produced only when the local cell density is high. It is unclear which selection pressures shaped the evolution of density-dependent regulation of toxin production. Here, we study the evolution of toxin production, resistance and the response to a cell-density cue in a model of an evolving bacterial population with spatial structure. We present results for two growth regimes: (i) an undisturbed, fixed habitat in which only small fluctuations of cell density occur, and (ii) a serial-transfer regime with large fluctuations in cell density. We find that density-dependent toxin production can evolve under both regimes. However, the selection pressures driving the evolution of regulation differ. In the fixed habitat, regulation evolves because it allows cells to produce toxin only when opportunities for reproduction are highly limited (because of a high local cell density), and the effective fitness costs of toxin production are hence low. Under serial transfers, regulation evolves because it allows cells to switch from a fast-growing non-toxic phenotype when colonising a new habitat, to a slower-growing competitive toxic phenotype when the cell density increases. Colonies of such regulating cells rapidly expand into unoccupied space because their edges consist of fast-growing, non-toxin-producing cells, but are also combative because cells at the interfaces with competing colonies do produce toxin. Because under the two growth regimes different types of regulation evolve, our results underscore the importance of growth conditions in the evolution of social behaviour in bacteria. [ABSTRACT FROM AUTHOR]

Visceral leishmaniasis (VL) is a zoonosis caused by the protozoan Leishmania infantum and in Brazil is transmitted mainly by the bite of Lutzomuyia longipalpis sand flies. Data about the presence, distribution, natural infection rate, seasonal and monthly dynamics of the vector population are important for optimizing the measures to control VL in endemic areas. This study aimed to identify sand fly fauna in an endemic area for VL to detect the prevalence of L. infantum infection in the Lu. longipalpis population and to elucidate the influence of bioclimatic factors on the monthly fluctuations of this vector. HP light traps were monthly set in the intradomicile and peridomicile of residences located in the central and beachfront areas of Camaçari, a VL endemic area. The sand fly collection was conducted in two periods: i) period 1—between December 2011 and November 2012 and ii) period 2—August 2014 and July 2015. Sand fly species were identified and detection of L. infantum infection by qPCR was performed in pools of female Lu. longipalpis. For the first time, the parasite load of positive pools was correlated with the number of Lu. longipalpis captured per month in both periods. Correlation analyses between the monthly fluctuation of the sand fly population and bioclimatic indices of the municipality in both collection periods were also performed. In both evaluated periods, more than 98% of the collected sand flies were Lu. longipalpis, confirming the predominance of this species in the region. It was captured mostly in the beachfront area in all months evaluated (99%). For the period 1, Leishmania DNA was detected in 81% of tested pools representing a minimal infection rate of 9.6%. In the period 2, 40% of the pools were positive with a minimal infection rate of 10.2%. Infected sand flies were only detected in the beachfront area in both periods. The parasite load was low and did not vary in the evaluated months despite the number of collected sand flies. No correlation was observed for climatic factors in both areas of Camaçari. These findings emphasize the high risk of Leishmania transmission in Camaçari regardless of the season and that other factors, aside from bioclimatic elements, are influencing the sand fly population monthly fluctuation. [ABSTRACT FROM AUTHOR]

Translocated saltwater crocodiles (Crocodylus porosus) in the Northern Territory (NT) of Australia often return to their original capture sites, which complicates management interventions aimed at reducing human-crocodile conflict. We examined the spatial events implicated in this homing ability, using ARGOS satellite tracking devices. Five large male C. porosus (3.03 m to 4.02 m TL) were shifted and released 100–320 km from their capture sites, and 3 additional ones (3.67 m to 4.23 m TL) were released at their site of capture as controls. Translocated crocodiles were more mobile than the controls, and moved at sea in the direction of their original capture site. However, they were unable or unwilling to swim around a geographic structure, Cobourg Peninsula, which prevented homing being achieved in all five cases. Two control crocodiles remained near their capture sites, but one, after the first year, made a 900km journey for six months, before returning to its original capture and release site. Genetic analysis of tissue samples from nests across the NT coast demonstrated significant genetic structure across the coast, and confirmed that Cobourg Peninsula contributes to genetic differentiation among populations along the NT coast. These results provide new insights into C. porosus movements, which have management significance for the maintenance of public safety. [ABSTRACT FROM AUTHOR]

Background: Genomic data have become major resources to understand complex mechanisms at fine-scale temporal and spatial resolution in functional and evolutionary genetic studies, including human diseases, such as cancers. Recently, a large number of whole genomes of evolving populations of yeast (Saccharomyces cerevisiae W303 strain) were sequenced in a time-dependent manner to identify temporal evolutionary patterns. For this type of study, a chromosome-level sequence assembly of the strain or population at time zero is required to compare with the genomes derived later. However, there is no fully automated computational approach in experimental evolution studies to establish the chromosome-level genome assembly using unique features of sequencing data. Methods and results: In this study, we developed a new software pipeline, the integrative meta-assembly pipeline (IMAP), to build chromosome-level genome sequence assemblies by generating and combining multiple initial assemblies using three de novo assemblers from short-read sequencing data. We significantly improved the continuity and accuracy of the genome assembly using a large collection of sequencing data and hybrid assembly approaches. We validated our pipeline by generating chromosome-level assemblies of yeast strains W303 and SK1, and compared our results with assemblies built using long-read sequencing and various assembly evaluation metrics. We also constructed chromosome-level sequence assemblies of S. cerevisiae strain Sigma1278b, and three commonly used fungal strains: Aspergillus nidulans A713, Neurospora crassa 73, and Thielavia terrestris CBS 492.74, for which long-read sequencing data are not yet available. Finally, we examined the effect of IMAP parameters, such as reference and resolution, on the quality of the final assembly of the yeast strains W303 and SK1. Conclusions: We developed a cost-effective pipeline to generate chromosome-level sequence assemblies using only short-read sequencing data. Our pipeline combines the strengths of reference-guided and meta-assembly approaches. Our pipeline is available online at including a Docker image, as well as a Perl script, to help users install the IMAP package, including several prerequisite programs. Users can use IMAP to easily build the chromosome-level assembly for the genome of their interest. [ABSTRACT FROM AUTHOR]

In a recent publication in PLOS ONE, Gabriele Margos and colleagues have questioned the division of the genus Borrelia into two genera on the basis that the differences in percentage of conserved proteins (POCP) between these two groups is >50%, which an earlier study has suggested as the threshold for differentiating prokaryotic genera. However, the POCP threshold is a poorly characterized and rarely used criterion for establishing distinction among prokaryotic genera. Detailed evaluation of the intergeneric POCP values for 37 genera from 3 different families (viz. Enterobacteriaceae- 24 genera, Morganellaceae-8 genera and Cystobacteraceae-5 genera) presented here shows that the POCP values for all genera within each of these families exceeded >58%. Thus, the suggested POCP threshold is not a useful criterion for delimitation of genus boundary and the objection by Margos et al. on this ground is invalid. Additionally, Margos et al. have questioned the specificities of ~15–20% of the conserved signature indels (CSIs) described in our work. However, as shown here, this concern is due to misunderstanding of the results and the CSIs in question are still highly-specific characteristics of the members of these genera and they provide important information regarding the evolutionary relationships of two new reptiles-echidna-related species viz. Borrelia turcica and Candidatus Borrelia tachyglossi to other Borrelia species. Results presented here show that both these species are deeper-branching members of the genus Borrelia and their placement within this genus is strongly supported by phylogenetic analyses and multiple uniquely shared CSIs with the other Borrelia species. Based on the large body of evidence derived from phylogenetic, genomic, molecular, phenotypic and clinical features, it is contended that the characteristics clearly distinguishing the Borrelia and Borreliella genera are far more numerous and of different kinds than those discerning most (all) other neighbouring genera of prokaryotes. Thus, the placement of these two groups of microorganisms into distinct genera, Borrelia and Borreliella, which clearly recognizes the differences among them, is highly appropriate and it should lead to a better understanding of the clinical, molecular and biological differences between these two important groups of microbes. [ABSTRACT FROM AUTHOR]

Changing environmental conditions, including those caused by human activities, reshape biological communities through both loss of native species and establishment of non-native species in the altered habitats. Dynamic interactions with the abiotic environment impact both immigration and initial establishment of non-native species into these altered habitats. The repeated emergence of disease systems in urban areas worldwide highlights the importance of understanding how dynamic migratory processes affect the current and future distribution and abundance of pathogens in urban environments. In this study, we examine the pattern of invasion of Trypanosoma cruzi—the causative agent of human Chagas disease—in the city of Arequipa, Peru. Phylogenetic analyses of 136 T. cruzi isolates from Arequipa and other South American locations suggest that only one T. cruzi lineage established a population in Arequipa as all T. cruzi isolated from vectors in Arequipa form a recent monophyletic group within the broader South American phylogeny. We discuss several hypotheses that may explain the limited number of established T. cruzi lineages despite multiple introductions of the parasite. [ABSTRACT FROM AUTHOR]

Cesky Fousek is considered to be one of the oldest pointing dog breeds in Europe and has been appreciated for its versatile working skills. Because it faced extinction in the past, the Cesky Fousek was restored from German Wirehaired and Shorthaired Pointers. Additionally, the breed was recently used in the USA with the initial intent of improvement of the Wirehaired Pointing Griffon (synonymous with Korthals Griffon) by the Bohemian Wirehaired Pointing Griffon Club of America. This study evaluates genetic diversity parameters of Cesky Fousek and compares them to the other continental pointing dogs that played a role in the formation of its gene pool. DNA from buccal swab and blood samples (n = 405) were analyzed using 18 microsatellite markers. Parameters of genetic polymorphism show that the Cesky Fousek breed has a comparable rate of variation as other hunting breeds despite the low population size and severe historical bottlenecks. Clustering analyses reveal a unique genetic status as a distinct pointing dog breed and the relatedness of the breeds is in good concordance with historical data. The present study demonstrates that despite historical admixture among lineages, separate pointing breeds constitute genetically differentiated units, mirroring unique breeding stocks and pedigree isolation among specific breed clubs, reflecting differences in breeding programs under each association. [ABSTRACT FROM AUTHOR]

In many beef and some dairy production systems, crossbreeding is used to take advantage of breed complementarity and heterosis. Admixed animals are frequently identified by their coat color and body conformation phenotypes, however, without pedigree information it is not possible to identify the expected breed composition of an admixed animal and in the presence of selection, the actual composition may differ from expectation. As the roles of DNA and genotype data become more pervasive in animal agriculture, a systematic method for estimating the breed composition (the proportions of an animal’s genome originating from ancestral pure breeds) has utility for a variety of downstream analyses including the estimation of genomic breeding values for crossbred animals, the estimation of quantitative trait locus effects, and heterosis and heterosis retention in advanced generation composite animals. Currently, there is no automated or semi-automated ancestry estimation platform for cattle and the objective of this study was to evaluate the utility of extant public software for ancestry estimation and determine the effects of reference population size and composition and number of utilized single nucleotide polymorphism loci on ancestry estimation. We also sought to develop an analysis pipeline that would simplify this process for members of the livestock genomics research community. We developed and tested a tool, “CRUMBLER”, to estimate the global ancestry of cattle using ADMIXTURE and SNPweights based on a defined reference panel. CRUMBLER, was developed and evaluated in cattle, but is a species agnostic pipeline that facilitates the streamlined estimation of breed composition for individuals with potentially complex ancestries using publicly available global ancestry software and a specified reference population SNP dataset. We developed the reference panel from a large cattle genotype data set and breed association pedigree information using iterative analyses to identify purebred individuals that were representative of each breed. We also evaluated the numbers of markers necessary for breed composition estimation and simulated genotypes for advanced generation composite animals to evaluate the precision of the developed tool. The developed CRUMBLER pipeline extracts a specified subset of genotypes that is common to all current commercially available genotyping platforms, processes these into the file formats required for the analysis software, and predicts admixture proportions using the specified reference population allele frequencies. [ABSTRACT FROM AUTHOR]

The synaptonemal complex (SC) is a proteinaceous scaffold required for synapsis and recombination between homologous chromosomes during meiosis. Although the SC has been linked to differences in genome-wide crossover rates, the genetic basis of standing variation in SC structure remains unknown. To investigate the possibility that recombination evolves through changes to the SC, we characterized the genetic architecture of SC divergence on two evolutionary timescales. Applying a novel digital image analysis technique to spermatocyte spreads, we measured total SC length in 9,532 spermatocytes from recombinant offspring of wild-derived mouse strains with differences in this fundamental meiotic trait. Using this large dataset, we identified the first known genomic regions involved in the evolution of SC length. Distinct loci affect total SC length divergence between and within subspecies, with the X chromosome contributing to both. Joint genetic analysis of MLH1 foci—immunofluorescent markers of crossovers—from the same spermatocytes revealed that two of the identified loci also confer differences in the genome-wide recombination rate. Causal mediation analysis suggested that one pleiotropic locus acts early in meiosis to designate crossovers prior to SC assembly, whereas a second locus primarily shapes crossover number through its effect on SC length. One genomic interval shapes the relationship between SC length and recombination rate, likely modulating the strength of crossover interference. Our findings pinpoint SC formation as a key step in the evolution of recombination and demonstrate the power of genetic mapping on standing variation in the context of the recombination pathway. [ABSTRACT FROM AUTHOR]

In recent years it became clear that in eukaryotic genome evolution gene loss is prevalent over gene gain. However, the absence of genes in an annotated genome is not always equivalent to the loss of genes. Due to sequencing issues, or incorrect gene prediction, genes can be falsely inferred as absent. This implies that loss estimates are overestimated and, more generally, that falsely inferred absences impact genomic comparative studies. However, reliable estimates of how prevalent this issue is are lacking. Here we quantified the impact of gene prediction on gene loss estimates in eukaryotes by analysing 209 phylogenetically diverse eukaryotic organisms and comparing their predicted proteomes to that of their respective six-frame translated genomes. We observe that 4.61% of domains per species were falsely inferred to be absent for Pfam domains predicted to have been present in the last eukaryotic common ancestor. Between phylogenetically different categories this estimate varies substantially: for clade-specific loss (ancestral loss) we found 1.30% and for species-specific loss 16.88% to be falsely inferred as absent. For BUSCO 1-to-1 orthologous families, 18.30% were falsely inferred to be absent. Finally, we showed that falsely inferred absences indeed impact loss estimates, with the number of losses decreasing by 11.78%. Our work strengthens the increasing number of studies showing that gene loss is an important factor in eukaryotic genome evolution. However, while we demonstrate that on average inferring gene absences from predicted proteomes is reliable, caution is warranted when inferring species-specific absences. [ABSTRACT FROM AUTHOR]

Splicing expands, reshapes, and regulates the transcriptome of eukaryotic organisms. Despite its importance, key questions remain unanswered, including the following: Can splicing evolve when organisms adapt to new challenges? How does evolution optimize inefficiency of introns’ splicing and of the splicing machinery? To explore these questions, we evolved yeast cells that were engineered to contain an inefficiently spliced intron inside a gene whose protein product was under selection for an increased expression level. We identified a combination of mutations in Cis (within the gene of interest) and in Trans (in mRNA-maturation machinery). Surprisingly, the mutations in Cis resided outside of known intronic functional sites and improved the intron’s splicing efficiency potentially by easing tight mRNA structures. One of these mutations hampered a protein’s domain that was not under selection, demonstrating the evolutionary flexibility of multi-domain proteins as one domain functionality was improved at the expense of the other domain. The Trans adaptations resided in two proteins, Npl3 and Gbp2, that bind pre-mRNAs and are central to their maturation. Interestingly, these mutations either increased or decreased the affinity of these proteins to mRNA, presumably allowing faster spliceosome recruitment or increased time before degradation of the pre-mRNAs, respectively. Altogether, our work reveals various mechanistic pathways toward optimizations of intron splicing to ultimately adapt gene expression patterns to novel demands. [ABSTRACT FROM AUTHOR]

Body mass index (BMI) is influenced by genetic, behavioral and environmental factors, while interactions between genetic and socioeconomic factors have been suggested. Aim of the study was to investigate whether socioeconomic position (SEP) interacts with a BMI-related genetic sum score (GRSBMI) to affect BMI in a population-based cohort. SEP-related health behaviors and a GRS associated with educational attainment (GRSEdu) were included in the analysis to explore potential interactions underlying the GRSBMIxSEP effect. Baseline information on SEP indicators (education, income), BMI, smoking, physical activity, alcohol consumption and genetic risk factors were available for 4,493 participants of the Heinz Nixdorf Recall Study. Interaction analysis was based on linear regression as well as on stratified analyses. In SEP-stratified analyses, the highest genetic effects were observed in the lowest educational group with a 0.24 kg/m2 higher BMI (95%CI: 0.16; 0.31) and in the lowest income quartile with a 0.14 kg/m2 higher BMI (95%CI: 0.09; 0.18) per additional risk allele. Indication for a GRSBMIxSEP interaction was observed for education (ßGRSbmixeducation = -0.02 [95%CI:-0.03; -0.01]) and income (ßGRSbmixincome = -0.05 [95%CI: -0.08; -0.02]). When adjusting for interactions with the GRSEdu and SEP-related health behaviors, effect size estimates of the GRSBMIxSEP interaction remained virtually unchanged. Results gave indication for an interaction of BMI-related genetic risk factors with SEP indicators, showing substantially stronger genetic effects in low SEP groups. This supports the hypothesis that expression of genetic risks is higher in socioeconomically disadvantaged environments. No indication was observed that the GRSBMIxSEP interaction was affected by other SEP-related factors included in the analysis. [ABSTRACT FROM AUTHOR]

Background: In the last four decades, the Asian tiger mosquito, Aedes albopictus, vector of several human arboviruses, has spread from its native range in South-East Asia to all over the world, largely through the transportation of its eggs via the international trade in used tires. Albania was the first country invaded in Europe in 1979, followed by Italy in 1990 and other Mediterranean countries after 2000. Methods/Principal findings: We here inferred the invasion history and migration patterns of Ae. albopictus in Italy (today the most heavily-infested country in Europe), Greece and Albania, by analyzing a panel of >100,000 single nucleotide polymorphisms (SNPs) obtained by sequencing of double-digest Restriction site-Associated DNA (ddRADseq). The obtained dataset was combined with samples previously analyzed from both the native and invasive range worldwide to interpret the results using a broader spatial and historical context. The emerging evolutionary scenario complements the results of other studies in showing that the extraordinary worldwide expansion of Ae. albopictus has occurred thanks to multiple independent invasions by large numbers of colonists from multiple geographic locations in both native and previously invaded areas, consistently with the role of used tires shipments to move large numbers of eggs worldwide. By analyzing mosquitoes from nine sites across ~1,000-km transect in Italy, we were able to detect a complex interplay of drift, isolation by distance mediated divergence, and gene flow in shaping the species very recent invasion and range expansion, suggesting overall high connectivity, likely due to passive transportation of adults via ground transportation, as well as specific adaptations to local conditions. Conclusions/Significance: Results contribute to characterize one of the most successful histories of animal invasion, and could be used as a baseline for future studies to track epidemiologically relevant characters (e.g. insecticide resistance). [ABSTRACT FROM AUTHOR]

Rapid improvements in DNA sequencing technology have resulted in long genome sequences for a large number of similar isolates with a wide range of single nucleotide polymorphism (SNP) rates, where some isolates can have thousands of times lower SNP rates than others. Genome sequences of this kind are a challenge to existing methods for construction of phylogenetic trees. We address the issues by developing a hierarchical approach to phylogeny construction. In this method, the construction is performed at multiple levels, where at each level, groups of isolates with similar levels of similarity are identified and their phylogenetic trees are constructed. Time savings are achieved by using a sufficiently large number of columns from the input alignment, instead of all its columns. Our results show that the new approach is 20-60 times more efficient than existing programs and more accurate in situations where highly similar isolates have a wide range of SNP rates. [ABSTRACT FROM AUTHOR]

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these “hungry spliceosome” conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of “intronization”, whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function. [ABSTRACT FROM AUTHOR]

Phenotypic plasticity is a vital strategy for plants to deal with changing conditions by inducing phenotypes favourable in different environments. Understanding how natural selection acts on variation in phenotypic plasticity in plants is therefore a central question in ecology, but is often ignored in modelling studies. Here we present a new modelling approach that allows for the analysis of selection for variation in phenotypic plasticity as a response strategy. We assess selection for shade avoidance strategies of Arabidopsis thaliana in response to future neighbour shading signalled through a decrease in red:far-red (R:FR) ratio. For this, we used a spatially explicit 3D virtual plant model that simulates individual Arabidopsis plants competing for light in different planting densities. Plant structure and growth were determined by the organ-specific interactions with the light environment created by the vegetation structure itself. Shade avoidance plastic responses were defined by a plastic response curve relating petiole elongation and lamina growth to R:FR perceived locally. Different plasticity strategies were represented by different shapes of the response curve that expressed different levels of R:FR sensitivity. Our analyses show that the shape of the selected shade avoidance strategy varies with planting density. At higher planting densities, more sensitive response curves are selected for than at lower densities. In addition, the balance between lamina and petiole responses influences the sensitivity of the response curves selected for. Combining computational virtual plant modelling with a game theoretical analysis represents a new step towards analysing how natural selection could have acted upon variation in shade avoidance as a response strategy, which can be linked to genetic variation and underlying physiological processes. [ABSTRACT FROM AUTHOR]

Gypsy moth Lymantria dispar L. 1758 (Lepidoptera: Erebidae) is one of the most dangerous forest pests of the Holarctic region. Outbreaks of gypsy moth populations lead to significant defoliation of local forests. Within the vast territory of the West Siberian Plain, we noted the outbreak front movement in the north-east direction with a speed 100–200 km per year. The reason for the outbreak's movement is still unclear because L. dispar females are characterised by limited flight ability, which is not enough to support that movement per se. Herein, we analysed the mtDNA divergence pattern among L. dispar populations collected from the vast territory of the West Siberian Plain to determine the boundaries of populations and reveal the effect of the outbreak's front movement on mtDNA patterns of populations. The 590-bp region of the cytochrome oxidase subunit I gene of the mitochondrial genome was sequenced for 220 specimens that were collected from 18 localities along a transect line (~ 1400 km). Our results clearly show that the gypsy moth populations of the vast Siberian territory are not subdivided. This result can be explained by extensive genetic exchange among local populations. Taking into account that the flight ability of L. dispar females is rather limited, we suggest that spreading occurs through ballooning of early instar larvae. This hypothesis was confirmed by the coincidence of the outbreaks' movement direction with that of the dominant winds, complemented by the observation of ballooned larvae far from a forest edge. [ABSTRACT FROM AUTHOR]

Gene expression analysis of rare or heterogeneous cell populations such as disseminated cancer cells (DCCs) requires a sensitive method allowing reliable analysis of single cells. Therefore, we developed and explored the feasibility of a quantitative PCR (qPCR) assay to analyze single-cell cDNA pre-amplified using a previously established whole transcriptome amplification (WTA) protocol. We carefully selected and optimized multiple steps of the protocol, e.g. re-amplification of WTA products, quantification of amplified cDNA yields and final qPCR quantification, to identify the most reliable and accurate workflow for quantitation of gene expression of the ERBB2 gene in DCCs. We found that absolute quantification outperforms relative quantification. We then validated the performance of our method on single cells of established breast cancer cell lines displaying distinct levels of HER2 protein. The different protein levels were faithfully reflected by transcript expression across the tested cell lines thereby proving the accuracy of our approach. Finally, we applied our method to breast cancer DCCs of a patient undergoing anti-HER2-directed therapy. Here, we were able to measure ERBB2 expression levels in all HER2-protein-positive DCCs. In summary, we developed a reliable single-cell qPCR assay applicable to measure distinct levels of ERBB2 in DCCs. [ABSTRACT FROM AUTHOR]

In this study, we report the plastome of Eriocaulon decemflorum (Eriocaulaceae) and make an effort to understand the genome evolution, structural rearrangements and gene content of the order Poales by comparing it with other available plastomes. The size of complete E. decemflorum plastome is 151,671 bp with an LSC (81,477bp), SSC (17,180bp) and a pair of IRs (26,507 bp). The plastome exhibits GC content of 35.8% and 134 protein-coding genes with 19 genes duplicated in the IR region. The Eriocaulaceae plastome is characterized by the presence of accD, ycf1 and ycf2 genes and presence of introns in clpP and rpoC1 genes which have been lost in the Graminid plastomes. Phylogenomic analysis based on 81 protein-coding genes placed Eriocaulaceae sister to Mayacaceae. The present study enhances our understanding of the evolution of Poales by analyzing the plastome data from the order. [ABSTRACT FROM AUTHOR]

In Brassica oleracea, heterosis is the most efficient tool providing impetus to hybrid vegetable industry. In this context, we presented the first report on identifying superior heterotic crosses for yield and commercial traits in cauliflower involving cytoplasmic male sterile (CMS) and doubled haploid (DH) lines as parents. We studied the suitability of genomic-SSRs and EST-SSRs based genetic distance (GD) and agronomic trait based phenotypic distance (PD) for predicting heterosis in F1 hybrids using CMS and DH based parents. 120 F1 hybrids derived from 20Ogura based CMS lines and 6 DH based testers were evaluated for 16 agronomic traits along with the 26 parental lines and 4 commercial standard checks. The genomic-SSRs and EST-SSRs based genetic structure analysis grouped the 26 parental lines into 4 distinct clusters. The CMS lines Ogu118-6A, Ogu33A, Ogu34-1A were good general combiner for developing early maturity hybrids. The SCA effects were significantly associated with heterosis suggesting non-additive gene effects for the heterotic response of hybrids. Less than unity value of σ2A/D coupled with σ2gca/σ2sca indicated the predominance of non-additive gene action in the expression of studied traits. The correlation analysis of genetic distance with heterosis for commercial traits suggested that microsatellites based genetic distance estimates can be helpful in heterosis prediction to some extent. [ABSTRACT FROM AUTHOR]

Allele age has long been a focus of population genetic research, primarily because it can be an important clue to the fitness effects of an allele. By virtue of their effects on fitness, alleles under directional selection are expected to be younger than neutral alleles of the same frequency. We developed a new coalescent-based estimator of a close proxy for allele age, the time when a copy of an allele first shares common ancestry with other chromosomes in a sample not carrying that allele. The estimator performs well, including for the very rarest of alleles that occur just once in a sample, with a bias that is typically negative. The estimator is mostly insensitive to population demography and to factors that can arise in population genomic pipelines, including the statistical phasing of chromosomes. Applications to 1000 Genomes Data and UK10K genome data confirm predictions that singleton alleles that alter proteins are significantly younger than those that do not, with a greater difference in the larger UK10K dataset, as expected. The 1000 Genomes populations varied markedly in their distributions for singleton allele ages, suggesting that these distributions can be used to inform models of demographic history, including recent events that are only revealed by their impacts on the ages of very rare alleles. [ABSTRACT FROM AUTHOR]

Drug resistance is a rapidly emerging concern, thus prompting the development of novel therapeutics or combinatorial therapy. Currently, combinatorial therapy targets are based on knowledge of drug mode of action and/or resistance mechanisms, constraining the number of target proteins. Unbiased genome-wide screens could reveal novel genetic components within interaction networks as potential targets in combination therapies. Testing this, in the context of antimicrobial resistance, we implemented an unbiased genome-wide screen, performed in Saccharomyces cerevisiae expressing a Candida glabrata PDR1+ gain-of-function allele. Gain-of-function mutations in this gene are the principal mediators of fluconazole resistance in this human fungal pathogen. Eighteen synthetically lethal S. cerevisiae genetic mutants were identified in cells expressing C. glabrata PDR1+. One mutant, lacking the histone acetyltransferase Gcn5, was investigated further. Deletion or drug-mediated inhibition of Gcn5 caused a lethal phenotype in C. glabrata cells expressing PDR1+ alleles. Moreover, deletion or drug-mediated inactivation of Gcn5, inhibited the emergence of fluconazole-resistant C. glabrata isolates in evolution experiments. Thus, taken together, the data generated in this study provides proof of concept that synthetically lethal genetic screens can identify novel candidate proteins that when therapeutically targeted could allow effective treatment of drug-resistant infections. [ABSTRACT FROM AUTHOR]

Arapaima, pirarucu or paiche (Arapaima gigas) is one of the largest freshwater fish in the world, and has a long history of commercial exploitation in the Amazon region. To estimate levels of genetic variability and historical and recent connectivity in Arapaima, we examined variation in eleven microsatellite DNA markers in individuals from 22 localities in Brazil, Colombia, and Peru. The results of analysis of molecular variance, Bayesian clustering and discriminant analysis of principal components showed that Arapaima in our samples represents two major populations, one in the Amazonas and one in the Araguaia-Tocantins River basins. The Amazonas population is further structured by isolation-by-distance with the hydrologically largely unconnected Amapá locality representing the eastern-most extreme of this continuum; gene flow predominates at distances of less than 1500 km with localities separated by over 2000 km dominated by genetic drift and effectively forming different populations. We saw no evidence of multiple species of Arapaima in the Amazonas basin, and analysis of pairwise genetic divergence (FST) with Mantel tests and correlograms indicated that this largest population exhibits a large-scale pattern of isolation-by-distance, with which results from MIGRATE-N agreed. The degree and significance of genetic divergence indicates that most sampled localities represent demographically independent sub-populations, although we did identify several recent migration events between both proximal and more distant localities. The levels of genetic diversity were heterogeneous across sites, including low genetic diversity, effective population sizes, and evidence of genetic bottlenecks in several places. On average the levels of gene diversity and rarefied allelic richness were higher for localities along the Amazonas mainstem than in the tributaries, despite these being the areas of highest fishing pressure, while the lowest values were found in tributary headwaters, where landscape modification is a significant threat. We recommend that managers consider the regional and local threats to these populations and tailor strategies accordingly, strategies which should ensure the ability of young A. gigas to disperse through floodplain corridors to maintain genetic diversity among otherwise sedentary adult sub-populations. [ABSTRACT FROM AUTHOR]

The bonnethead shark, Sphyrna tiburo, is a small elasmobranch distributed in the Eastern Pacific from southern California to Ecuador, and along the Western Atlantic, with preferences for continental margins of North, Central and South America, the Gulf of Mexico, and the Caribbean. Recent studies have suggested that it could be under a process of cryptic speciation, with the possibility to find different species in similar geographic locations. Here we assessed the population structure and genetic diversity of this highly philopatric and non-dispersive species in the Bocas del Toro Archipelago, Panama. Fragments of the mitochondrial genes cytochrome oxidase I and control region, were used to test the genetic structure of adult and juvenile S. tiburo in this area, and were compared with other locations of the Western Atlantic and Belize. We found significant genetic differentiation between Caribbean bonnethead sharks from Bocas del Toro and Belize, when compared to bonnetheads from other locations of the Western Atlantic. These results also suggest that Bocas del Toro could constitute a different genetic population unit for this species, whereby bonnethead sharks in this area could belong to a unique stock. The information obtained in this study could improve our understanding of the population dynamics of the bonnethead shark throughout its distribution range, and may be used as a baseline for future conservation initiatives for coastal sharks in Central America, a poorly studied an often overlooked region for shark conservation and research. [ABSTRACT FROM AUTHOR]

We describe a new genome alignment-based model for understanding the diversity of viruses based on evolutionary genetic relationships. This approach uses information theory and a physical model to determine the information shared by the genes in two genomes. Pairwise comparisons of genes from the viruses are created from alignments using NCBI BLAST, and their match scores are combined to produce a metric between genomes, which is in turn used to determine a global classification using the 5,817 viruses on RefSeq. In cases where there is no measurable alignment between any genes, the method falls back to a coarser measure of genome relationship: the mutual information of 4-mer frequency. This results in a principled model which depends only on the genome sequence, which captures many interesting relationships between viral families, and which creates clusters which correlate well with both the Baltimore and ICTV classifications. The incremental computational cost of classifying a novel virus is low and therefore newly discovered viruses can be quickly identified and classified. The model goes beyond alignment-free classifications by producing a full phylogeny similar to those constructed by virologists using qualitative features, while relying only on objective genes. These results bolster the case for mathematical models in microbiology which can characterize organisms using only their genetic material and provide an independent check for phylogenies constructed by humans, considerably faster and more cheaply than less modern approaches. [ABSTRACT FROM AUTHOR]