Differences In Gene Copy Number Research Articles

Identification and biochemical characterization of a carboxylesterase gene associated with β-cypermethrin resistance in Dermanyssus gallinae

Dermanyssus gallinae is a major hematophagous ectoparasite in layer hens. Although the acaricide β-cypermethrin has been used to control mites worldwide, D. gallinae has developed resistance to this compound. Carboxylesterases (CarEs) are important detoxification enzymes that confer resistance to β-cypermethrin in arthropods. However, CarEs associated with β-cypermethrin resistance in D. gallinae have not yet been functionally characterized. Here, we isolated a CarE gene (Deg-CarE) from D. gallinae and assayed its activity. The results revealed significantly higher expression of Deg-CarE in the β-cypermethrin-resistant strain (RS) than in the susceptible strain (SS) toward α-naphthyl acetate (α-NA) and β-naphthyl acetate (β-NA). These findings suggest that enhanced esterase activities might have contributed to β-cypermethrin resistance in D. gallinae. Quantitative real-time PCR analysis revealed that Deg-CarE expression levels were significantly higher in adults than in other life stages. Although Deg-CarE was upregulated in the RS, significant differences in gene copy numbers were not observed. Additionally, Deg-CarE expression was significantly induced by β-cypermethrin in both the SS and RS. Moreover, silencing Deg-CarE via RNA interference decreased the enzyme activity and increased the susceptibility of the RS to β-cypermethrin, confirming that Deg-CarE is crucial for β-cypermethrin detoxification. Finally, recombinant Deg-CarE (rDeg-CarE) expressed in Escherichia coli displayed high enzymatic activity toward α/β-NA. However, metabolic analysis indicated that rDeg-CarE did not directly metabolize β-cypermethrin. The collective findings indicate that D. gallinae resistance to β-cypermethrin is associated with elevated CarEs protein activity and increased Deg-CarE expression levels. These findings provide insights into the metabolic resistance of D. gallinae and offer scientific guidance for the management and control of D. gallinae.

Development of Primers Specific for Detection of Litylenchus crenatae, the Causal Agent of Beech Leaf Disease, in Plant Tissue.

Beech leaf disease (BLD), an emerging threat to American beech (Fagus grandifolia) in the northern United States and Canada, was recently confirmed to be caused by the nematode Litylenchus crenatae subsp. mccannii (hereafter L. crenatae). Consequently, there is a need for a rapid, sensitive, and accurate method for detecting L. crenatae for both diagnostic as well as control purposes. This research developed a new set of DNA primers that specifically amplify L. crenatae and allow for accurate detection of the nematode in plant tissue. These primers have also been used in quantitative PCR (qPCR) to determine relative differences in gene copy number between samples. This primer set provides an improved, effective tool for monitoring and detecting L. crenatae in temperate tree leaf tissue which is necessary to understand the spread of this emerging forest pest and to develop management strategies.

Comparing antimicrobial resistant genes and phenotypes across multiple sequencing platforms and assays for Enterobacterales clinical isolates

IntroductionWhole genome sequencing (WGS) of bacterial isolates can be used to identify antimicrobial resistance (AMR) genes. Previous studies have shown that genotype-based AMR has variable accuracy for predicting carbapenem resistance in carbapenem-resistant Enterobacterales (CRE); however, the majority of these studies used short-read platforms (e.g. Illumina) to generate sequence data. In this study, our objective was to determine whether Oxford Nanopore Technologies (ONT) long-read WGS would improve detection of carbapenem AMR genes with respect to short-read only WGS for nine clinical CRE samples. We measured the minimum inhibitory breakpoint (MIC) using two phenotype assays (MicroScan and ETEST) for six antibiotics, including two carbapenems (meropenem and ertapenem) and four non-carbapenems (gentamicin, ciprofloxacin, cefepime, and trimethoprim/sulfamethoxazole). We generated short-read data using the Illumina NextSeq and long-read data using the ONT MinION. Four assembly methods were compared: ONT-only assembly; ONT-only assembly plus short-read polish; ONT + short-read hybrid assembly plus short-read polish; short-read only assembly.ResultsConsistent with previous studies, our results suggest that the hybrid assembly produced the highest quality results as measured by gene completeness and contig circularization. However, ONT-only methods had minimal impact on the detection of AMR genes and plasmids compared to short-read methods, although, notably, differences in gene copy number differed between methods. All four assembly methods showed identical presence/absence of the blaKPC-2 carbapenemase gene for all samples. The two phenotype assays showed 100% concordant results for the non-carbapenems, but only 65% concordance for the two carbapenems. The presence/absence of AMR genes was 100% concordant with AMR phenotypes for all four non-carbapenem drugs, although only 22%—50% sensitivity for the carbapenems.ConclusionsOverall, these findings suggest that the lack of complete correspondence between CRE AMR genotype and phenotype for carbapenems, while concerning, is independent of sequencing platform/assembly method.

Population differentiated copy number variation between Eurasian wild boar and domesticated pig populations

Sus scrofa is a globally distributed livestock species that still maintains two different ways of life: wild and domesticated. Herein, we detected copy number variation (CNV) of 328 animals using short read alignment on Sscrofa11.1. We compared CNV among five groups of porcine populations: Asian domesticated (AD), European domesticated (ED), Asian wild (AW), European wild (EW), and Near Eastern wild (NEW). In total, 21,673 genes were identified on 154,872 copy number variation region (CNVR). Differences in gene copy numbers between populations were measured by considering the variance-based value {V}_{ST} and the one-way ANOVA test followed by Scheffe test. As a result, 111 genes were suggested as copy number variable genes. Abnormally gained copy number on EEA1 in all populations was suggested the presence of minor CNV in the reference genome assembly, Sscrofa11.1. Copy number variable genes were related to meat quality, immune response, and reproduction traits. Hierarchical clustering of all individuals and mean pairwise {V}_{ST} in breed level were visualized genetic relationship of 328 individuals and 56 populations separately. Our findings have shown how the complex history of pig evolution appears in genome-wide CNV of various populations with different regions and lifestyles.

Disparate Inventories of Hypoxia Gene Sets Across Corals Align With Inferred Environmental Resilience

Aquatic deoxygenation has been flagged as an overlooked but key factor contributing to mass bleaching-induced coral mortality. During deoxygenation events triggered by coastal nutrient pollution and ocean warming, oxygen supplies lower to concentrations that can elicit an aerobic metabolic crisis i.e., hypoxia. Surprisingly little is known of the fundamental hypoxia gene set inventory that corals possess to respond to lowered oxygen (i.e., deoxygenation). For instance, it is unclear whether gene copy number differences exist across species that may affect the efficacy of a measured transcriptomic stress response. Therefore, we conducted an ortholog-based meta-analysis to investigate how hypoxia gene inventories differ amongst coral species to assess putative copy number variations (CNVs). We specifically elucidated CNVs for a compiled list of 32 hypoxia genes across 24 protein sets from species with a sequenced genome spanning corals from the robust and complex clade. We found approximately a third of the investigated genes exhibited copy number differences, and these differences were species-specific rather than attributable to the robust-complex split. Interestingly, we consistently found the highest gene expansion present in Porites lutea, which is considered to exhibit inherently greater stress tolerance than other species. Consequently, our analysis suggests that hypoxia stress gene expansion may coincide with increased stress tolerance. As such, the unevenly expanded (or reduced) hypoxia genes presented here provide key genes of interest to target in examining (or diagnosing) coral stress responses. Important next steps will involve determining to what extent such gene copy differences align with certain coral traits.

Contribution of multiple overexpressed carboxylesterase genes to indoxacarb resistance in Spodoptera litura.

As an important family of detoxification enzymes, carboxylesterases (CarEs) have important roles in the development of insecticide resistance in almost all agricultural pests. Previous studies have suggested that enhancement of CarE activity is an important mechanism mediating indoxacarb resistance in Spodoptera litura, and several CarE genes have been found to be overexpressed in indoxacarb-resistant strains. However, the functions of these CarE genes in indoxacarb resistance needs to be further investigated. The synergist triphenyl phosphate effectively reduced the resistance of S. litura to indoxacarb, suggesting an involvement of CarEs in indoxacarb resistance. Among seven identified S. litura CarE genes (hereafter SlituCOE), six were overexpressed in two indoxacarb-resistant strains, but there were no significant differences in gene copy number. Knockdown of SlituCOE009 and SlituCOE050 enhanced indoxacarb sensitivity in both susceptible and resistant strains, whereas knockdown of SlituCOE090, SlituCOE093 and SlituCOE074 enhanced indoxacarb sensitivity in only the resistant strain. Knockdown of the sixth gene, SlituCOE073, did not have any effect. Furthermore, simultaneous knockdown of the five SlituCOE genes had a greater effect on increasing indoxacarb sensitivity than silencing them individually. By contrast, overexpression of the five SlituCOE genes individually in Drosophila melanogaster significantly decreased the toxicity of indoxacarb to transgenic fruit flies. Furthermore, modeling and docking analysis indicated that the catalytic pockets of SlituCOE009 and SlituCOE074 were ideally shaped for indoxacarb and N-decarbomethoxylated metabolite (DCJW), but the binding affinity for DCJW was stronger than for indoxacarb. This study reveals that multiple overexpressed CarE genes are involved in indoxacarb resistance in S. litura.

Evolution of dosage compensation does not depend on genomic background.

Organisms have evolved various mechanisms to cope with the differences in the gene copy numbers between sexes caused by degeneration of Y and W sex chromosomes. Complete dosage compensation or at least expression balance between sexes has been reported predominantly in XX/XY systems, but rarely in ZZ/ZW systems. However, this often-reported pattern is based on comparisons of lineages where sex chromosomes evolved from nonhomologous genomic regions, potentially differing in sensitivity to differences in gene copy numbers. Here we document that two reptilian lineages (XX/XY iguanas and ZZ/ZW softshell turtles), which independently co-opted the same ancestral genomic region for the function of sex chromosomes, evolved different gene dose regulatory mechanisms. The independent co-option of the same genomic region for the role of sex chromosomes as in the iguanas and the softshell turtles offers great opportunity for testing evolutionary scenarios on sex chromosome evolution under the explicit control of the genomic background and gene identity. We show that the parallel loss of functional genes from the Y chromosome of the green anole and the W chromosome of the Florida softshell turtle led to different dosage compensation mechanisms. Our approach controlling for genetic background thus does not support that the variability in the regulation of gene dose differences is a consequence of ancestral autosomal gene content.

Measurement of gene amplifications related to drug resistance in Plasmodium falciparum using droplet digital PCR

BackgroundCopy number variations (CNVs) of the Plasmodium falciparum multidrug resistance 1 (pfmdr1), P. falciparum plasmepsin2 (pfplasmepsin2) and P. falciparum GTP cyclohydrolase 1 (pfgch1) genes are associated with anti-malarial drug resistance in P. falciparum malaria. Droplet digital PCR (ddPCR) assays have been developed for accurate assessment of CNVs in several human genes. The aim of the present study was to develop and validate ddPCR assays for detection of the CNVs of P. falciparum genes associated with resistance to anti-malarial drugs.MethodsA multiplex ddPCR assay was developed to detect the CNVs in the pfmdr1 and pfplasmepsin2 genes, while a duplex ddPCR assay was developed to detect CNV in the pfgch1 gene. The gene copy number (GCN) quantification limit, as well as the accuracy and precision of the ddPCR assays were determined and compared to conventional quantitative PCR (qPCR). In order to reduce the cost of testing, a multiplex ddPCR assay of two target genes, pfmdr1 and pfplasmepsin2, was validated. In addition, the CNVs of genes of field samples collected from Thailand from 2015 to 2019 (n = 84) were assessed by ddPCR and results were compared to qPCR as the reference assay.ResultsThere were no significant differences between the GCN results obtained from uniplex and multiplex ddPCR assays for detection of CNVs in the pfmdr1 and pfplasmepsin2 genes (p = 0.363 and 0.330, respectively). Based on the obtained gene copy number quantification limit, the accuracy and percent relative standard deviation (%RSD) value of the multiplex ddPCR assay were 95% and 5%, respectively, for detection of the CNV of the pfmdr1 gene, and 91% and 5% for detection of the CNV of the pfplasmepsin2 gene. There was no significant difference in gene copy numbers assessed by uniplex or duplex ddPCR assays regarding CNV in the pfgch1 gene (p = 0.276). The accuracy and %RSD value of the duplex ddPCR assay were 95% and 4%, respectively, regarding pfgch1 GCN. In the P. falciparum field samples, pfmdr1 and pfplasmepsin2 GCNs were amplified in 15% and 27% of samples from Ubon Ratchathani, Thailand, while pfgch1 GCN was amplified in 50% of samples from Yala, Thailand. There was 100% agreement between the GCN results obtained from the ddPCR and qPCR assays (κ = 1.00). The results suggested that multiplex ddPCR assay is the optional assay for the accurate detection of gene copy number without requiring calibration standards, while the cost and required time are reduced. Based on the results of this study, criteria for GCN detection by ddPCR analysis were generated.ConclusionsThe developed ddPCR assays are simple, accurate, precise and cost-effective tools for detection of the CNVs in the pfmdr1, pfplasmepsin2 and pfgch1 genes of P. falciparum. The ddPCR assay is a useful additional tool for the surveillance of anti-malarial drug resistance.

Expansion and Accelerated Evolution of 9-Exon Odorant Receptors in Polistes Paper Wasps.

Independent origins of sociality in bees and ants are associated with independent expansions of particular odorant receptor (OR) gene subfamilies. In ants, one clade within the OR gene family, the 9-exon subfamily, has dramatically expanded. These receptors detect cuticular hydrocarbons (CHCs), key social signaling molecules in insects. It is unclear to what extent 9-exon OR subfamily expansion is associated with the independent evolution of sociality across Hymenoptera, warranting studies of taxa with independently derived social behavior. Here, we describe OR gene family evolution in the northern paper wasp, Polistes fuscatus, and compare it to four additional paper wasp species spanning ∼40 million years of evolutionary divergence. We find 200 putatively functional OR genes in P. fuscatus, matching predictions from neuroanatomy, and more than half of these are in the 9-exon subfamily. Most OR gene expansions are tandemly arrayed at orthologous loci in Polistes genomes, and microsynteny analysis shows species-specific gain and loss of 9-exon ORs within tandem arrays. There is evidence of episodic positive diversifying selection shaping ORs in expanded subfamilies. Values of omega (dN/dS) are higher among 9-exon ORs compared to other OR subfamilies. Within the Polistes OR gene tree, branches in the 9-exon OR clade experience relaxed negative (relaxed purifying) selection relative to other branches in the tree. Patterns of OR evolution within Polistes are consistent with 9-exon OR function in CHC perception by combinatorial coding, with both natural selection and neutral drift contributing to interspecies differences in gene copy number and sequence.

Complete genome sequence of lovastatin producer Aspergillus terreus ATCC 20542 and evaluation of genomic diversity among A. terreus strains

In the present study, the complete genome of a filamentous fungus Aspergillus terreus ATCC 20542 was sequenced, assembled, and annotated. This strain is mainly recognized for being a model wild-type lovastatin producer and a parental strain of high-yielding industrial mutants. It is also a microorganism with a rich repertoire of secondary metabolites that has been a subject of numerous bioprocess-related studies. In terms of continuity, the genomic sequence provided in this work is of the highest quality among all the publicly available genomes of A. terreus strains. The comparative analysis revealed considerable diversity with regard to the catalog of biosynthetic gene clusters found in A. terreus. Even though the cluster of lovastatin biosynthesis was found to be well-conserved at the species level, several unique genes putatively associated with metabolic functions were detected in A. terreus ATCC 20542 that were not detected in other investigated genomes. The analysis was conducted also in the context of the primary metabolic pathways (sugar catabolism, biomass degradation potential, organic acid production), where the visible differences in gene copy numbers were detected. However, the species-level genomic diversity of A. terreus was more evident for secondary metabolism than for the well-conserved primary metabolic pathways. The newly sequenced genome of A. terreus ATCC 20542 was found to harbor several unique sequences, which can be regarded as interesting subjects for future experimental efforts on A. terreus metabolism and fungal biosynthetic capabilities.Key points• The high-quality genome of Aspergillus terreus ATCC 20542 has been assembled and annotated.• Comparative analysis with other sequenced Aspergillus terreus strains has revealed considerable diversity in biosynthetic gene repertoire, especially related to secondary metabolism.• The unique genomic features of A. terreus ATCC 20542 are discussed.

Independent Evolution of Sex Chromosomes in Eublepharid Geckos, A Lineage with Environmental and Genotypic Sex Determination.

Geckos demonstrate a remarkable variability in sex determination systems, but our limited knowledge prohibits accurate conclusions on the evolution of sex determination in this group. Eyelid geckos (Eublepharidae) are of particular interest, as they encompass species with both environmental and genotypic sex determination. We identified for the first time the X-specific gene content in the Yucatán banded gecko, Coleonyx elegans, possessing X1X1X2X2/X1X2Y multiple sex chromosomes by comparative genome coverage analysis between sexes. The X-specific gene content of Coleonyx elegans was revealed to be partially homologous to genomic regions linked to the chicken autosomes 1, 6 and 11. A qPCR-based test was applied to validate a subset of X-specific genes by comparing the difference in gene copy numbers between sexes, and to explore the homology of sex chromosomes across eleven eublepharid, two phyllodactylid and one sphaerodactylid species. Homologous sex chromosomes are shared between Coleonyx elegans and Coleonyx mitratus, two species diverged approximately 34 million years ago, but not with other tested species. As far as we know, the X-specific gene content of Coleonyx elegans / Coleonyx mitratus was never involved in the sex chromosomes of other gecko lineages, indicating that the sex chromosomes in this clade of eublepharid geckos evolved independently.

Differential Activation of Immune Effector Processes in Mature Compared to Immature Sacrococcygeal Teratomas

Objective: This study aims to characterize the molecular signatures of sacrococcygeal teratomas (SCTs). Methods: Three SCTs were analyzed with whole genome genotyping. RNA sequencing of 10 SCTs dominated by mature, immature and neuroglial elements was analyzed. Expression in SCT-samples with different elements were compared to each other and to a reference group of malignant pediatric tumors. Macrophages, T- and B-lymphocytes were detected by immunohistochemistry. Results: No chromosomal imbalances were detected. SCTs showed overexpression of genes involved in neurosignaling, DNA-binding molecules and pathways of early germ cells. Genes associated with immune effector processes were overexpressed in mature compared to immature SCTs, and immune cell infiltration was found predominantly around mature epithelial elements. Conclusion: The broad repertoire of histological elements in SCTs reflects differences in transcriptional regulation rather than differences in gene copy numbers. A paucity of immune response in immature SCTs may be a factor contributing to their uninhibited growth.

Gene expression and target-site mutations are associated with resistance to ALS inhibitors in annual sedge (Cyperus compressus) biotypes from Georgia

AbstractAnnual sedge (Cyperus compressus L.) populations with resistance to halosulfuron were identified in turfgrass at two new locations in Georgia. Research was conducted to evaluate (1) resistance levels to two acetolactate synthase (ALS) inhibitors, (2) ALS enzyme susceptibility, (3) genetic differences associated with resistance, and (4) differential levels of ALS gene expression in these biotypes. In dose–response experiments, the biotypes were >160 times resistant to halosulfuron but only 12 times more resistant to imazaquin compared with a susceptible biotype. In vitro enzyme assays indicated that resistant (R) biotypes required 6.1-fold greater concentrations of imazaquin to reduce ALS activity 50% compared with the susceptible (S) biotype. Both R biotypes had a similar Pro-197-Ser amino acid substitution in the ALS gene that confers resistance to sulfonylureas. Compared with the S biotype, R biotypes had 4.4 times higher ALS gene expression than the S biotype. No differences in gene copy number were found between any biotypes for the ALS gene. Overall, ALS-resistant C. compressus selected by halosulfuron use in turfgrass may be exhibiting partial susceptibility to imazaquin but complete resistance to sulfonylureas. Differential levels of susceptibility to ALS inhibitors for these biotypes are associated with the Pro-197-Ser substitution and enhanced expression of the ALS gene.

Characterisation of the gut microbial community of rhesus macaques in high-altitude environments

BackgroundThe mammal intestinal microbiota is involved in various physiological processes and plays a key role in host environment adaption. However, for non-human primates (NHPs), little is known about their gut microbial community in high-altitude environments and even less about their adaption to such habitats. We characterised the gut microbial community of rhesus macaques from multiple high-altitude environments and compared it to those of low-altitude populations.ResultsWe collected faecal samples of rhesus macaques from four high-altitude populations (above 3000 m) and three low-altitude populations (below 500 m). By calculating the alpha diversity index, we found that high-altitude populations exhibited a higher diversity. Statistical analysis of beta diversity indicated significant differences between high- and low-altitude populations. Significant differences were also detected at the phylum and family levels. At the phylum level, the high-altitude gut microbial community was dominated by Firmicutes (63.42%), while at low altitudes, it was dominated by Bacteroidetes (47.4%). At the family level, the high-altitude population was dominated by Ruminococcaceae (36.2%), while the low-altitude one was dominated by Prevotellaceae (39.6%). Some families, such as Christensenellaceae and Rikenellaceae, were consistently higher abundant in all high-altitude populations. We analysed the overlap of operational taxonomic units (OTUs) in high-altitude populations and determined their core OTUs (shared by all four high-altitude populations). However, when compared with the low-altitude core OTUs, only 65% were shared, suggesting a divergence in core OTUs. Function prediction indicated a significant difference in gene copy number of 35 level-2 pathways between high- and low-altitude populations; 29 of them were higher in high altitudes, especially in membrane transport and carbohydrate metabolism.ConclusionsThe gut microbial community of high-altitude rhesus macaques was significantly distinct from that of low-altitude populations in terms of diversity, composition and function. High-altitude populations were dominated by Firmicutes and Ruminococcace, while in low-altitude populations, Bacteroidetes and Prevotellaceae were dominant. The difference in gut microbiota between these two populations may be caused by differences in host diet, environmental temperature and oxygen pressure. These differentiated gut microbial microorganisms may play a critical role in the adaptive evolution of rhesus macaques to high-altitude environments.

Polar bear evolution is marked by rapid changes in gene copy number in response to dietary shift

Polar bear (Ursus maritimus) and brown bear (Ursus arctos) are recently diverged species that inhabit vastly differing habitats. Thus, analysis of the polar bear and brown bear genomes represents a unique opportunity to investigate the evolutionary mechanisms and genetic underpinnings of rapid ecological adaptation in mammals. Copy number (CN) differences in genomic regions between closely related species can underlie adaptive phenotypes and this form of genetic variation has not been explored in the context of polar bear evolution. Here, we analyzed the CN profiles of 17 polar bears, 9 brown bears, and 2 black bears (Ursus americanus). We identified an average of 318 genes per individual that showed evidence of CN variation (CNV). Nearly 200 genes displayed species-specific CN differences between polar bear and brown bear species. Principal component analysis of gene CN provides strong evidence that CNV evolved rapidly in the polar bear lineage and mainly resulted in CN loss. Olfactory receptors composed 47% of CN differentiated genes, with the majority of these genes being at lower CN in the polar bear. Additionally, we found significantly fewer copies of several genes involved in fatty acid metabolism as well as AMY1B, the salivary amylase-encoding gene in the polar bear. These results suggest that natural selection shaped patterns of CNV in response to the transition from an omnivorous to primarily carnivorous diet during polar bear evolution. Our analyses of CNV shed light on the genomic underpinnings of ecological adaptation during polar bear evolution.

Comparative Evolution of Sand Fly Salivary Protein Families and Implications for Biomarkers of Vector Exposure and Salivary Vaccine Candidates

Sand fly salivary proteins that produce a specific antibody response in humans and animal reservoirs have been shown to be promising biomarkers of sand fly exposure. Furthermore, immunity to sand fly salivary proteins were shown to protect rodents and non-human primates against Leishmania infection. We are missing critical information regarding the divergence amongst sand fly salivary proteins from different sand fly vectors, a knowledge that will support the search of broad or specific salivary biomarkers of vector exposure and those for vaccines components against leishmaniasis. Here, we compare the molecular evolution of the salivary protein families in New World and Old World sand flies from 14 different sand fly vectors. We found that the protein families unique to OW sand flies are more conserved than those unique to NW sand flies regarding both sequence polymorphisms and copy number variation. In addition, the protein families unique to OW sand flies do not display as many conserved cysteine residues as the one unique to the NW group (28.5% in OW vs. 62.5% in NW). Moreover, the expression of specific protein families is restricted to the salivary glands of unique sand fly taxon. For instance, the ParSP15 family is unique to the Larroussius subgenus whereas phospholipase A2 is only expressed in member of Larroussius and Adlerius subgenera. The SP2.5-like family is only expressed in members of the Phlebotomus and Paraphlebotomus subgenera. The sequences shared between OW and NW sand flies have diverged at similar rates (38.7 and 45.3% amino acid divergence, respectively), yet differences in gene copy number were evident across protein families and sand fly species. Overall, this comparative analysis sheds light on the different modes of sand fly salivary protein family divergence. Also, it informs which protein families are unique and conserved within taxon for the choice of taxon-specific biomarkers of vector exposure, as well as those families more conserved across taxa to be used as pan-specific vaccines for leishmaniasis.

LAPTM4B gene copy number gain is associated with inferior response to anthracycline-based chemotherapy in hormone receptor negative breast carcinomas.

To determine the associations between lysosomal-associated transmembrane protein 4b (LAPTM4B) gene copy number and response to different chemotherapy regimens in hormone receptor negative (HR-) primary breast carcinomas. Two cohorts were analyzed: (1) 69 core biopsies from HR-breast carcinomas treated with neoadjuvant chemotherapy (anthracycline based in 72.5% of patients and non-anthracycline based in 27.5% of patients). (2) Tissue microarray (TMA) of 74 HR-breast carcinomas treated with adjuvant therapy (77.0% of the patients received anthracycline, 17.6% of the patients non-anthracycline-based therapy, and in 5.4% of the cases, no treatment data are available). Interphase FISH technique was applied on pretreatment core biopsies (cohort I) and on TMAs (cohort II) using custom-made dual-labelled FISH probes (LAPTM4B/CEN8q FISH probe Abnova Corp.). In the neoadjuvant cohort in the anthracycline-treated group, we observed a significant difference (p = 0.029) of average LAPTM4B copy number between the non-responder and pathological complete responder groups (4.1 ± 1.1 vs. 2.6 ± 0.1). In the adjuvant setting, the anthracycline-treated group of metastatic breast carcinomas was characterized by higher LAPTM4B copy number comparing to the non-metastatic ones (p = 0.046). In contrast, in the non-anthracycline-treated group of patients, we did not find any LAPTM4B gene copy number differences between responder vs. non-responder groups or between metastatic vs. non-metastatic groups. Our results confirm the possible role of the LAPTM4B gene in anthracycline resistance in HR- breast cancer. Analyzing LAPTM4B copy number pattern may support future treatment decision.

VITCOMIC2: visualization tool for the phylogenetic composition of microbial communities based on 16S rRNA gene amplicons and metagenomic shotgun sequencing

BackgroundThe 16S rRNA gene-based amplicon sequencing analysis is widely used to determine the taxonomic composition of microbial communities. Once the taxonomic composition of each community is obtained, evolutionary relationships among taxa are inferred by a phylogenetic tree. Thus, the combined representation of taxonomic composition and phylogenetic relationships among taxa is a powerful method for understanding microbial community structure; however, applying phylogenetic tree-based representation with information on the abundance of thousands or more taxa in each community is a difficult task. For this purpose, we previously developed the tool VITCOMIC (VIsualization tool for Taxonomic COmpositions of MIcrobial Community), which is based on the genome-sequenced microbes’ phylogenetic information. Here, we introduce VITCOMIC2, which incorporates substantive improvements over VITCOMIC that were necessary to address several issues associated with 16S rRNA gene-based analysis of microbial communities.ResultsWe developed VITCOMIC2 to provide (i) sequence identity searches against broad reference taxa including uncultured taxa; (ii) normalization of 16S rRNA gene copy number differences among taxa; (iii) rapid sequence identity searches by applying the graphics processing unit-based sequence identity search tool CLAST; (iv) accurate taxonomic composition inference and nearly full-length 16S rRNA gene sequence reconstructions for metagenomic shotgun sequencing; and (v) an interactive user interface for simultaneous representation of the taxonomic composition of microbial communities and phylogenetic relationships among taxa. We validated the accuracy of processes (ii) and (iv) by using metagenomic shotgun sequencing data from a mock microbial community.ConclusionsThe improvements incorporated into VITCOMIC2 enable users to acquire an intuitive understanding of microbial community composition based on the 16S rRNA gene sequence data obtained from both metagenomic shotgun and amplicon sequencing.

Between-species differences in gene copy number are enriched among functions critical for adaptive evolution in Arabidopsis halleri.

BackgroundGene copy number divergence between species is a form of genetic polymorphism that contributes significantly to both genome size and phenotypic variation. In plants, copy number expansions of single genes were implicated in cultivar- or species-specific tolerance of high levels of soil boron, aluminium or calamine-type heavy metals, respectively. Arabidopsis halleri is a zinc- and cadmium-hyperaccumulating extremophile species capable of growing on heavy-metal contaminated, toxic soils. In contrast, its non-accumulating sister species A. lyrata and the closely related reference model species A. thaliana exhibit merely basal metal tolerance.ResultsFor a genome-wide assessment of the role of copy number divergence (CND) in lineage-specific environmental adaptation, we conducted cross-species array comparative genome hybridizations of three plant species and developed a global signal scaling procedure to adjust for sequence divergence. In A. halleri, transition metal homeostasis functions are enriched twofold among the genes detected as copy number expanded. Moreover, biotic stress functions including mostly disease Resistance (R) gene-related genes are enriched twofold among genes detected as copy number reduced, when compared to the abundance of these functions among all genes.ConclusionsOur results provide genome-wide support for a link between evolutionary adaptation and CND in A. halleri as shown previously for Heavy metal ATPase4. Moreover our results support the hypothesis that elemental defences, which result from the hyperaccumulation of toxic metals, allow the reduction of classical defences against biotic stress as a trade-off.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3319-5) contains supplementary material, which is available to authorized users.

Evaluation of peat and sawdust as permeable reactive barrier materials for stimulating in situ biodegradation of trichloroethene

Two low cost solid organic materials, sawdust and peat, were tested in laboratory batch microcosm and flow-through column experiments to determine their suitability for application in permeable reactive barriers (PRBs) supporting biodegradation of trichloroethene (TCE). In microcosms with peat, TCE (∼30μM) was sequentially and completely degraded to cis-dichloroethene (cDCE), vinyl chloride, and ethene through reductive dechlorination. In microcosms with sawdust, reductive dechlorination of TCE stopped at cDCE and high methane production (up to 3000μM) was observed. 16S rRNA gene copy numbers of Dehalobacter and Archaea were higher (1000 and 10 times, respectively) in sawdust microcosms than those in peat microcosms. Dehalococcoides and vcrA gene copy numbers were 10 times higher in peat microcosms than in sawdust microcosms. These gene copy number differences are consistent with the extent of TCE degradation and production of methane in the microcosms. Flow-through column experiments showed that hydraulic conductivity reduction with time was consistently greater in the sawdust column compared to the peat column. The greater conductivity reduction was likely due to biofouling and methane gas bubble formation. The experimental observations indicate that peat has potential to be a better solid organic material than sawdust to support reductive dechlorination of TCE in PRB applications.