Genetic diversity and virulence of phytopathogenic Burkholderia glumae strains isolated from rice cultivars in valleys of the high jungle of Peru.

Burkholderia glumae causes bacterial leaf blight in rice, and its global spread has been exacerbated by climate change. To understand the genetic diversity and virulence of B. glumae strains isolated from rice cultivars in Peru, 47 isolates were obtained from infected rice fields, all belonging to B. glumae, and confirmed by recA and toxB sequences. The BOX-PCR typing group 38 genomic profiles, and these turn into 7 Variable Number Tandem Repeats (VNTR) haplotypes. There was no correlation between clustering and geographical origin. Nineteen strains were selected for phenotypic characterization and virulence, using both the maceration level of the onion bulb proxy and inoculation of seeds of two rice cultivars. Several strains produced pigments other than toxoflavin, which correlated with onion bulb maceration. In terms of virulence at the seed level, all strains produced inhibition at the root and coleoptile level, but the severity of symptoms varied significantly between strains, revealing significant differences in pathogenicity. There is no correlation between maceration and virulence scores, probably reflecting different virulence mechanisms depending on the host infection stage. This is the first study to evaluate the VNTR diversity and virulence of Peruvian strains of B. glumae in two commercial cultivars.

Unlocking the Complete Chloroplast Genome of a Native Tree Species from the Amazon Basin, Capirona (Calycophyllum Spruceanum, Rubiaceae), and Its Comparative Analysis with Other Ixoroideae Species.

Capirona (Calycophyllum spruceanum Benth.) belongs to subfamily Ixoroideae, one of the major lineages in the Rubiaceae family, and is an important timber tree. It originated in the Amazon Basin and has widespread distribution in Bolivia, Peru, Colombia, and Brazil. In this study, we obtained the first complete chloroplast (cp) genome of capirona from the department of Madre de Dios located in the Peruvian Amazon. High-quality genomic DNA was used to construct libraries. Pair-end clean reads were obtained by PE 150 library and the Illumina HiSeq 2500 platform. The complete cp genome of C. spruceanum has a 154,480 bp in length with typical quadripartite structure, containing a large single copy (LSC) region (84,813 bp) and a small single-copy (SSC) region (18,101 bp), separated by two inverted repeat (IR) regions (25,783 bp). The annotation of C. spruceanum cp genome predicted 87 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, 37 transfer RNA (tRNA) genes, and one pseudogene. A total of 41 simple sequence repeats (SSR) of this cp genome were divided into mononucleotides (29), dinucleotides (5), trinucleotides (3), and tetranucleotides (4). Most of these repeats were distributed in the noncoding regions. Whole chloroplast genome comparison with the other six Ixoroideae species revealed that the small single copy and large single copy regions showed more divergence than inverted regions. Finally, phylogenetic analyses resolved that C. spruceanum is a sister species to Emmenopterys henryi and confirms its position within the subfamily Ixoroideae. This study reports for the first time the genome organization, gene content, and structural features of the chloroplast genome of C. spruceanum, providing valuable information for genetic and evolutionary studies in the genus Calycophyllum and beyond.

Characterization of the complete chloroplast genome of a Peruvian landrace of Capsicum chinense Jacq. (Solanaceae), arnaucho chili pepper.

In this study, we sequenced the first complete chloroplast (cp) genome of a Peruvian chili pepper landrace, "arnacucho" (Capsicum chinense). This cp genome has a 156,931 bp in length with typical quadripartite structure, containing a large single copy (LSC) region (87,325 bp) and a 17,912 bp small single-copy (SSC) region, separated by two inverted repeat (IR) regions (25,847 bp); and the percentage of GC content was 37.71%. Arnaucho chili pepper chloroplast genome possesses 133 genes that consists of 86 protein-coding genes, 37 tRNA, eight rRNA, and two pseudogenes. Phylogenetic analysis revealed that this Peruvian chili pepper landrace is closely related to the undomesticated species C. galapagoense; all belong to the Capsiceae tribe.

Chloroplast genome of Tillandsia landbeckii Phil. (Bromeliaceae) a species adapted to the hyper-arid conditions of the Atacama and Peruvian desert.

Tillandsia landbeckii Phil. is a vulnerable species belonging to the terrestrial Bromeliaceae family; it is highly adapted to extremely hyper-arid conditions of the Atacama Desert and Peruvian deserts. In this study, we sequenced, assembled, and annotated its chloroplast genome. T. landbeckii chloroplast genome is 159,131 bp in length, containing a large single-copy region of 87,164 bp, a small single-copy region of 18,521 bp, and a pair of inverted repeat regions of 26,723 bp. The GC content of the chloroplast genome is 37.33%. It encodes a total of 132 genes, including 86 protein-coding genes, 38 tRNA genes and 8 rRNA genes. The phylogenetic tree indicates that T. landbeckii is placed within the Bromeliaceae family and has a close relationship with T. marconae with 100% support.

The complete chloroplast genome of the national tree of Peru, quina (Cinchona officinalis L., Rubiaceae).

Here, we report the first complete chloroplast (cp) genome of Cinchona officinalis. This cp genome has a 156,984 bp in length with typical quadripartite structure, containing a large single copy (LSC) region (83,929 bp) and an 18,051 bp small single-copy (SSC) region, separated by two inverted repeat (IR) regions (27,502 bp). The total GC content was 37.75%. Quina tree chloroplast genome possesses 135 genes that consisted of 89 protein-coding genes, 37 tRNA, eight rRNA, and one pseudogene. Phylogenetic analysis showed that C. officinalis is sister to C. pubescens and sister to them is Isertia laevis; all belong to the Cinchonoideae sub-family.

Chloroplast genome of Tillandsia marconae till & Vitek (Bromeliaceae), a hyperarid desert endangered species.

Tillandsia marconae Till & Vitek (Bromeliaceae) is a rare plant native species that grows over sand, in the coastal desert from Perú and Chile and is considered an endangered species. In this study, we assembled its chloroplast genome. The draft chloroplast genome of T. marconae is ca. 158,873 bp in length, containing a large single-copy region of 86,937 bp, a small single-copy region of 18,506 bp, and a pair of inverted repeat regions of 26,715 bp. The GC content of the draft chloroplast genome is 37.4%. It encodes a total of 135 genes, including 86 protein-coding genes, 38 tRNA genes, 8 rRNA genes, and three pseudogenes. The phylogenetic tree indicated that T. marconae is placed within the Bromeliaceae family and a close relationship with Tillandsia usneoides with 100% support.

Mitochondrial DNA Variation in Peruvian Honey Bee (Apis mellifera L.) Populations Using the tRNAleu-cox2 Intergenic Region.

Mitochondrial DNA variations of Peruvian honey bee populations were surveyed by using the tRNAleu-cox2 intergenic region. Only two studies have characterized these populations, indicating the presence of Africanized honey bee colonies in different regions of Peru and varied levels of Africanization, but the current status of its genetic diversity is unknown. A total of 512 honey bee colonies were sampled from three regions to characterize them. Our results revealed the presence of European and African haplotypes: the African haplotypes identified belong to sub-lineage AI (13) and sub-lineage AIII (03), and the European haplotypes to lineages C (06) and M (02). Of 24 haplotypes identified, 15 new sequences are reported here (11 sub-lineage AI, 2 sub-lineage AIII, and 2 lineage M). Peruvian honey bee populations presented a higher proportion from African than European haplotypes. High proportions of African haplotype were reported for Piura and Junín, unlike Lima, which showed more European haplotypes from lineage C. Few colonies belonging to lineage M would represent accidental purchase or traces of the introduction into Peru in the 19th century.

The honeybee (Apis mellifera) developmental state shapes the genetic composition of the deformed wing virus-A quasispecies during serial transmission.

The main biological threat to the western honeybee (Apis mellifera) is the parasitic mite Varroa destructor, largely because it vectors lethal epidemics of honeybee viruses that, in the absence of this mite, are relatively innocuous. The severe pathology is a direct consequence of excessive virus titres caused by this novel transmission route. However, little is known about how the virus adapts genetically during transmission and whether this influences the pathology. Here, we show that upon injection into honeybee pupae, the deformed wing virus type-A (DWV-A) quasispecies undergoes a rapid, extensive expansion of its sequence space, followed by strong negative selection towards a uniform, common shape by the time the pupae have completed their development, with no difference between symptomatic and asymptomatic adults in either DWV titre or genetic composition. This suggests that the physiological and molecular environment during pupal development has a strong, conservative influence on shaping the DWV-A quasispecies in emerging adults. There was furthermore no evidence of any progressive adaptation of the DWV-A quasispecies to serial intra-abdominal injection, simulating mite transmission, despite the generation of ample variation immediately following each transmission, suggesting that the virus either had already adapted to transmission by injection, or was unaffected by it.

Wing Geometric Morphometrics of Workers and Drones and Single Nucleotide Polymorphisms Provide Similar Genetic Structure in the Iberian Honey Bee (Apis mellifera iberiensis).

Wing geometric morphometrics has been applied to honey bees (Apis mellifera) in identification of evolutionary lineages or subspecies and, to a lesser extent, in assessing genetic structure within subspecies. Due to bias in the production of sterile females (workers) in a colony, most studies have used workers leaving the males (drones) as a neglected group. However, considering their importance as reproductive individuals, the use of drones should be incorporated in these analyses in order to better understand diversity patterns and underlying evolutionary processes. Here, we assessed the usefulness of drone wings, as well as the power of wing geometric morphometrics, in capturing the signature of complex evolutionary processes by examining wing shape data, integrated with geographical information, from 711 colonies sampled across the entire distributional range of Apis mellifera iberiensis in Iberia. We compared the genetic patterns reconstructed from spatially-explicit shape variation extracted from wings of both sexes with that previously reported using 383 genome-wide SNPs (single nucleotide polymorphisms). Our results indicate that the spatial structure retrieved from wings of drones and workers was similar (r = 0.93) and congruent with that inferred from SNPs (r = 0.90 for drones; r = 0.87 for workers), corroborating the clinal pattern that has been described for A. m. iberiensis using other genetic markers. In addition to showing that drone wings carry valuable genetic information, this study highlights the capability of wing geometric morphometrics in capturing complex genetic patterns, offering a reliable and low-cost alternative for preliminary estimation of population structure.

Draft Genome Sequence of Bacillus thuringiensis Strain UNMSM10RA, Isolated from Potato Crop Soil in Peru.

The 5.5-Mb genome sequence of Bacillus thuringiensis strain UNMSM10RA, isolated from potato crop soil, is reported in this study. The strain UNMSM10RA contains 5,347 protein-coding sequences, 105 tRNA genes, 15 rRNA genes, and 5 noncoding RNA (ncRNA) genes, with an average G+C content of 35.1%.

Whole genome SNP-associated signatures of local adaptation in honeybees of the Iberian Peninsula.

The availability of powerful high-throughput genomic tools, combined with genome scans, has helped identifying genes and genetic changes responsible for environmental adaptation in many organisms, including the honeybee. Here, we resequenced 87 whole genomes of the honeybee native to Iberia and used conceptually different selection methods (Samßada, LFMM, PCAdapt, iHs) together with in sillico protein modelling to search for selection footprints along environmental gradients. We found 670 outlier SNPs, most of which associated with precipitation, longitude and latitude. Over 88.7% SNPs laid outside exons and there was a significant enrichment in regions adjacent to exons and UTRs. Enrichment was also detected in exonic regions. Furthermore, in silico protein modelling suggests that several non-synonymous SNPs are likely direct targets of selection, as they lead to amino acid replacements in functionally important sites of proteins. We identified genomic signatures of local adaptation in 140 genes, many of which are putatively implicated in fitness-related functions such as reproduction, immunity, olfaction, lipid biosynthesis and circadian clock. Our genome scan suggests that local adaptation in the Iberian honeybee involves variations in regions that might alter patterns of gene expression and in protein-coding genes, which are promising candidates to underpin adaptive change in the honeybee.

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera).

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M-lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C-lineage). Effective conservation actions require reliable molecular tools to identify pure-bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q-values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.

Revisiting the Iberian honey bee (Apis mellifera iberiensis) contact zone: maternal and genome-wide nuclear variations provide support for secondary contact from historical refugia.

Dissecting diversity patterns of organisms endemic to Iberia has been truly challenging for a variety of taxa, and the Iberian honey bee is no exception. Surveys of genetic variation in the Iberian honey bee are among the most extensive for any honey bee subspecies. From these, differential and complex patterns of diversity have emerged, which have yet to be fully resolved. Here, we used a genome-wide data set of 309 neutrally tested single nucleotide polymorphisms (SNPs), scattered across the 16 honey bee chromosomes, which were genotyped in 711 haploid males. These SNPs were analysed along with an intergenic locus of the mtDNA, to reveal historical patterns of population structure across the entire range of the Iberian honey bee. Overall, patterns of population structure inferred from nuclear loci by multiple clustering approaches and geographic cline analysis were consistent with two major clusters forming a well-defined cline that bisects Iberia along a northeastern-southwestern axis, a pattern that remarkably parallels that of the mtDNA. While a mechanism of primary intergradation or isolation by distance could explain the observed clinal variation, our results are more consistent with an alternative model of secondary contact between divergent populations previously isolated in glacial refugia, as proposed for a growing list of other Iberian taxa. Despite current intense honey bee management, human-mediated processes have seemingly played a minor role in shaping Iberian honey bee genetic structure. This study highlights the complexity of the Iberian honey bee patterns and reinforces the importance of Iberia as a reservoir of Apis mellifera diversity.

Reduced SNP panels for genetic identification and introgression analysis in the dark honey bee (Apis mellifera mellifera).

Beekeeping activities, especially queen trading, have shaped the distribution of honey bee (Apis mellifera) subspecies in Europe, and have resulted in extensive introductions of two eastern European C-lineage subspecies (A. m. ligustica and A. m. carnica) into the native range of the M-lineage A. m. mellifera subspecies in Western Europe. As a consequence, replacement and gene flow between native and commercial populations have occurred at varying levels across western European populations. Genetic identification and introgression analysis using molecular markers is an important tool for management and conservation of honey bee subspecies. Previous studies have monitored introgression by using microsatellite, PCR-RFLP markers and most recently, high density assays using single nucleotide polymorphism (SNP) markers. While the latter are almost prohibitively expensive, the information gained to date can be exploited to create a reduced panel containing the most ancestry-informative markers (AIMs) for those purposes with very little loss of information. The objective of this study was to design reduced panels of AIMs to verify the origin of A. m. mellifera individuals and to provide accurate estimates of the level of C-lineage introgression into their genome. The discriminant power of the SNPs using a variety of metrics and approaches including the Weir & Cockerham's FST, an FST-based outlier test, Delta, informativeness (In), and PCA was evaluated. This study shows that reduced AIMs panels assign individuals to the correct origin and calculates the admixture level with a high degree of accuracy. These panels provide an essential tool in Europe for genetic stock identification and estimation of admixture levels which can assist management strategies and monitor honey bee conservation programs.

Signatures of selection in the Iberian honey bee (Apis mellifera iberiensis) revealed by a genome scan analysis of single nucleotide polymorphisms.

Understanding the genetic mechanisms of adaptive population divergence is one of the most fundamental endeavours in evolutionary biology and is becoming increasingly important as it will allow predictions about how organisms will respond to global environmental crisis. This is particularly important for the honey bee, a species of unquestionable ecological and economical importance that has been exposed to increasing human-mediated selection pressures. Here, we conducted a single nucleotide polymorphism (SNP)-based genome scan in honey bees collected across an environmental gradient in Iberia and used four FST -based outlier tests to identify genomic regions exhibiting signatures of selection. Additionally, we analysed associations between genetic and environmental data for the identification of factors that might be correlated or act as selective pressures. With these approaches, 4.4% (17 of 383) of outlier loci were cross-validated by four FST -based methods, and 8.9% (34 of 383) were cross-validated by at least three methods. Of the 34 outliers, 15 were found to be strongly associated with one or more environmental variables. Further support for selection, provided by functional genomic information, was particularly compelling for SNP outliers mapped to different genes putatively involved in the same function such as vision, xenobiotic detoxification and innate immune response. This study enabled a more rigorous consideration of selection as the underlying cause of diversity patterns in Iberian honey bees, representing an important first step towards the identification of polymorphisms implicated in local adaptation and possibly in response to recent human-mediated environmental changes.

Virulence attenuation of Candida albicans genetic variants isolated from a patient with a recurrent bloodstream infection.

The incidence of Candida albicans infections and the relapse episodes after antifungal treatment have increased in recent decades. Recurrences are mainly due to the persistence of the original infecting strain that may present genetic and genomic rearrangements during interaction with the host, reflecting strain adaptation. In this study, four isolates recovered from a patient during recurrent candidemia episodes were genotyped by microsatellite length polymorphism (MLP) and by multilocus sequence typing (MLST) and found to be genetic variants of the same strain. Using experimental mouse infections, a progressive reduction in the virulence of the four isolates was observed, with the first two isolates more virulent than the third and fourth. Additionally, in the mouse model, the first isolate resisted host control more efficiently, resulting in higher kidney fungal burdens and necrosis as compared to the third isolate. The resolution of inflammation was delayed in mice challenged with the first isolate and the message for IFN-gamma and TNF-alpha in the spleen was lower within the first few hours post-infection. Original and recurrent isolates also displayed different phenotypes regarding activity of secreted enzymes and response to stress agents. Overall, the comparative analysis indicated that the virulence decrease of these isolates was related to a lower ability to resist to the host anticandida effector mechanisms. We showed for the first time that C. albicans genetic variants of the same strain, sequentially isolated from an immunocompromised patient, underwent adaptations in the human host that resulted in virulence attenuation when tested in mice.

Microsatellite typing identifies the major clades of the human pathogen Candida albicans.

Candida albicans population studies showed that this species could be divided into sub-groups of closely related strains, designated by clades. Since the emergence of microsatellite analysis as a PCR based method, this technique has been successfully used as a tool to differentiate C. albicans isolates but has never been tested regarding clustering of the five major clades. In this study we tested microsatellite length polymorphism (MLP) ability to group 29 C. albicans strains previously defined as belonging to clades I, II, III, E and SA, nine atypical strains from Angola and Madagascar, and 78 Portuguese clinical isolates. MLP typing of the total 116 strains analyzed yielded 87 different multilocus allelic combinations which resulted in a high discriminatory power index, of 0.987, with only two markers, CA1 and CEF3. Cluster analysis of the 29 previously defined strains grouped them according to their clade designation with a matrix cophenetic correlation of r=0.963 after a normalized Mantel statistic. Clustering analysis of the 116 strains maintained the same groupings, clearly defining the five major C. albicans clades. The cophenetic value obtained was of r=0.839, and the one-tail probability of the normalized Mantel statistic out of 1000 random permutations was P=0.0020. The proportion of Portuguese isolates in the groups I, II, III and SA was of 2.7%, 15.4%, 3.8% and 0%, respectively. None of the isolates co-clustered with the atypical strains. These results confirm MLP typing as a good method both to type and differentiate C. albicans isolates and to group isolates, identifying the major C. albicans clades, similarly to Ca3 fingerprinting and multilocus sequence typing (MLST).