Species-specific Clades Research Articles

Comparative genomic analysis of chemosensory-related gene families in gastropods

Chemoreception is critical for the survival and reproduction of animals. Except for a reduced group of insects and chelicerates, the molecular identity of chemosensory proteins is poorly understood in invertebrates. Gastropoda is the extant mollusk class with the greatest species richness, including marine, freshwater, and terrestrial lineages, and likely, highly diverse chemoreception systems. Here, we performed a comprehensive comparative genome analysis taking advantage of the chromosome-level information of two Gastropoda species, one of which belongs to a lineage that underwent a whole genome duplication event. We identified thousands of previously uncharacterized chemosensory-related genes, the majority of them encoding G protein-coupled receptors (GPCR), mostly organized into clusters distributed across all chromosomes. We also detected gene families encoding degenerin epithelial sodium channels (DEG-ENaC), ionotropic receptors (IR), sensory neuron membrane proteins (SNMP), Niemann–Pick type C2 (NPC2) proteins, and lipocalins, although with a lower number of members. Our phylogenetic analysis of the GPCR gene family across protostomes revealed: (i) remarkable gene family expansions in Gastropoda; (ii) clades including members from all protostomes; and (iii) species-specific clades with a substantial number of receptors. For the first time, we provide new and valuable knowledge into the evolution of the chemosensory gene families in invertebrates other than arthropods.

Phylogenetic Relationships and Next-Generation Barcodes in the Genus Torreya Reveal a High Proportion of Misidentified Cultivated Plants.

Accurate species identification is key to conservation and phylogenetic inference. Living plant collections from botanical gardens/arboretum are important resources for the purpose of scientific research, but the proportion of cultivated plant misidentification are un-tested using DNA barcodes. Here, we assembled the next-generation barcode (complete plastid genome and complete nrDNA cistron) and mitochondrial genes from genome skimming data of Torreya species with multiple accessions for each species to test the species discrimination and the misidentification proportion of cultivated plants used in Torreya studies. A total of 38 accessions were included for analyses, representing all nine recognized species of genus Torreya. The plastid phylogeny showed that all 21 wild samples formed species-specific clades, except T. jiulongshanensis. Disregarding this putative hybrid, seven recognized species sampled here were successfully discriminated by the plastid genome. Only the T. nucifera accessions grouped into two grades. The species identification rate of the nrDNA cistron was 62.5%. The Skmer analysis based on nuclear reads from genome skims showed promise for species identification with seven species discriminated. The proportion of misidentified cultivated plants from arboreta/botanical gardens was relatively high with four accessions (23.5%) representing three species. Interspecific relationships within Torreya were fully resolved with maximum support by plastomes, where Torreya jackii was on the earliest diverging branch, though sister to T. grandis in the nrDNA cistron tree, suggesting that this is likely a hybrid species between T. grandis and an extinct Torreya ancestor lineage. The findings here provide quantitative insights into the usage of cultivated samples for phylogenetic study.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae).

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes and its only congeneric species, Platycarya strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole-genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae).

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes, and its only congeneric species, P. strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics, and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Evolution and functional analysis of the GRAS family genes in six Rosaceae species

BackgroundGRAS genes formed one of the important transcription factor gene families in plants, had been identified in several plant species. The family genes were involved in plant growth, development, and stress resistance. However, the comparative analysis of GRAS genes in Rosaceae species was insufficient.ResultsIn this study, a total of 333 GRAS genes were identified in six Rosaceae species, including 51 in strawberry (Fragaria vesca), 78 in apple (Malus domestica), 41 in black raspberry (Rubus occidentalis), 59 in European pear (Pyrus communis), 56 in Chinese rose (Rosa chinensis), and 48 in peach (Prunus persica). Motif analysis showed the VHIID domain, SAW motif, LR I region, and PFYRE motif were considerably conserved in the six Rosaceae species. All GRAS genes were divided into 10 subgroups according to phylogenetic analysis. A total of 15 species-specific duplicated clades and 3 lineage-specific duplicated clades were identified in six Rosaceae species. Chromosomal localization presented the uneven distribution of GRAS genes in six Rosaceae species. Duplication events contributed to the expression of the GRAS genes, and Ka/Ks analysis suggested the purification selection as a major force during the evolution process in six Rosaceae species. Cis-acting elements and GO analysis revealed that most of the GRAS genes were associated with various environmental stress in six Rosaceae species. Coexpression network analysis showed the mutual regulatory relationship between GRAS and bZIP genes, suggesting the ability of the GRAS gene to regulate abiotic stress in woodland strawberry. The expression pattern elucidated the transcriptional levels of FvGRAS genes in various tissues and the drought and salt stress in woodland strawberry, which were verified by RT-qPCR analysis.ConclusionsThe evolution and functional analysis of GRAS genes provided insights into the further understanding of GRAS genes on the abiotic stress of Rosaceae species.

Adaptive Evolution of the Fox Coronavirus Based on Genome-Wide Sequence Analysis.

Purpose To report the first complete fox coronavirus (CoV) genome sequence obtained through genome-wide amplifications and to understand the adaptive evolution of fox CoV. Methods Anal swab samples were collected from 35 foxes to detect the presence of CoV and obtain the virus sequence. Phylogenetic analysis was conducted using MrBayes. The possibility of recombination within these sequences was assessed using GARD. Analysis of the levels of selection pressure experienced by these sequences was assessed using methods on both the PAML and Data Monkey platforms. Results Of the 35 samples, two were positive, and complete genome sequences for the viruses were obtained. Phylogenetic analysis, using Bayesian methods, of these sequences, together with other CoV sequences, revealed that the fox CoV sequences clustered with canine coronavirus (CCoV) sequences, with sequences from other carnivores more distantly related. In contrast to the feline, ferret and mink CoV sequences that clustered into species-specific clades, the fox CoV fell within the CCoV clade. Minimal evidence for recombination was found among the sequences. A total of 7, 3, 14, and 2 positively selected sites were identified in the M, N, S, and 7B genes, respectively, with 99, 111, and 581 negatively selected sites identified in M, N, and S genes, respectively. Conclusion The complete genome sequence of fox CoV has been obtained for the first time. The results suggest that the genome sequence of fox CoV may have experienced adaptive evolution in the genes replication, entry, and virulence. The number of sites in each gene that experienced negative selection is far greater than the number that underwent positive selection, suggesting that most of the sequence is highly conserved and important for viral survive. However, positive selection at a few sites likely aided these viruses to adapt to new environments.

Molecular signatures of Janthinobacterium lividum from Trinidad support high potential for crude oil metabolism

BackgroundJanthinobacterium lividum is considered to be a psychrotrophic bacterial species. For the first time in the literature, J. lividum strains were isolated from Trinidad presenting with atypical features - hydrocarbonoclastic and able to survive in a tropical environment.MethodsIdentification of the Trinidad strains was carried out through 16S rRNA phylogenetic analysis. Gene-specific primers were designed to target the VioA which encodes violacein pigment and the EstA/B gene which encodes secreted extracellular lipase. Bioinformatics analyses were carried out on the nucleotide and amino acid sequences of VioA and EstA/B genes of the Trinidad Janthinobacterium strains to assess functionality and phylogenetic relatedness to other Janthinobacterium sequences specifically and more broadly, to other members of the Oxalobacteraceae family of betaproteobacteria.Results16S rRNA confirmed the identity of the Trinidad strains as J. lividum and resolved three of the Trinidad strains at the intra-specific level. Typical motility patterns of this species were recorded. VioAp sequences were highly conserved, however, synonymous substitutions located outside of the critical sites for enzyme function were detected for the Trinidad strains. Comparisons with PDB 6g2p model from aa231 to aa406 further indicated no functional disruption of the VioA gene of the Trinidad strains. Phylogeny of the VioA protein sequences inferred placement of all J. lividum taxa into a highly supported species-specific clade (bs = 98%). EstA/Bp sequences were highly conserved, however, synonymous substitutions were detected that were unique to the Trinidad strains. Phylogenetic inference positioned the Trinidad consensus VioA and EstA protein sequences in a clearly distinct branch.ConclusionsThe findings showed that the primary sequence of VioAp and EstA/Bp were unique to the Trinidad strains and these molecular signatures were reflected in phylogenetic inference. Our results supported chemotaxis, possible elective inactivation of VioA gene expression and secreted lipase activity as survival mechanisms of the Trinidad strains in petrogenic conditions.

Identification and expression analysis of cuticular protein genes in the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae)

Structural cuticular proteins (CPs) are major components of the insect cuticle, and they play critical roles in insect development and insecticide resistance. Here, a total of 196 CP genes were successfully annotated in the Plutella xylostella genome. On the basis of motif analysis, these CPs were classified into 10 different families, including 122 CPR, 12 CPAP1, 8 CPAP3, 9 CPLCP, 2 Tweedle, 1 CPF, 1 CPFL, 1 CPCFC, 17 CPG and 2 18 aa proteins, and the remaining 21 unclassified CPs were classed as cuticular proteins hypothetical (CPH). A phylogenetic analysis of CPs from different insects revealed species-specific clades of RR-1 and RR-2 genes, suggesting that CP gene duplication might occur independently among insect taxa, while we also found that some other CPs (such as CPAP1 and CPAP3) had a closer relationship based on their conserved domain architecture. Using available RNAseq libraries, the expression profiles of the CPs were analyzed over the four developmental stages of the insect (i.e., egg, larva, pupa, and adult), revealing stage-specific expression patterns for the CPs. In a chlorpyrifos resistant strain, 18 CP genes were found to be more than two-fold upregulated compared to the susceptible control strain, and qRT-PCR analysis showed that these CP genes were overexpressed after exposure to chlorpyrifos, suggesting a potential role in the molecular mechanism of insecticide resistance in P. xylostella. This study provides the tools and molecular basis to study the role of CPs in the post-embryonal development and the mechanisms of insecticide resistance of P. xylostella.

Viromics of extant insect orders unveil the evolution of the flavi-like superfamily.

Insects are the most diversified and species-rich group of animals and harbor an immense diversity of viruses. Several taxa in the flavi-like superfamily, such as the genus Flavivirus, are associated with insects; however, systematic studies on insect virus genetic diversity are lacking, limiting our understanding of the evolution of the flavi-like superfamily. Here, we examined the diversity of flavi-like viruses within the most complete and up-to-date insect transcriptome collection comprising 1,243 insect species by employing a Flaviviridae RdRp profile hidden Markov model search. We identified seventy-six viral sequences in sixty-one species belonging to seventeen insect, one entognathan, and one arachnidan orders. Phylogenetic analyses revealed that twenty-seven sequences fell within the Flaviviridae phylogeny but did not group with established genera. Despite the large diversity of insect hosts studied, we only detected one virus in a blood-feeding insect, which branched within the genus Flavivirus, indicating that this genus likely diversified only in hematophagous arthropods. Nine new jingmenviruses with novel host associations were identified. One of the jingmenviruses established a deep rooting lineage additional to the insect- and tick-associated clades. Segment co-segregation phylogenies support the separation of tick- and insect-associated groups within jingmenviruses, with evidence for segment reassortment. In addition, fourteen viruses grouped with unclassified flaviviruses encompassing genome length of up to 20 kb. Species-specific clades for Hymenopteran- and Orthopteran-associated viruses were identified. Forty-nine viruses populated three highly diversified clades in distant relationship to Tombusviridae, a plant-infecting virus family, suggesting the detection of three previously unknown insect-associated families that contributed to tombusvirus evolution.

Phylogenetic Analysis and Molecular Diversity of Capsicum Based on rDNA-ITS Region

The genus Capsicum is comprised of 5 domesticated and more than 30 wild species. The region of nuclear ribosomal DNA internal transcribed spacers (rDNA-ITS) has widely been used for species identification, but has rarely been used in Capsicum. In this study, the evaluation of genetic diversity and a phylogenetic analysis were conducted using rDNA-ITS of 28 Capsicum accessions, including five domesticated and two wild species. We surveyed six conventional keys of domesticated species and another five traits in Capsicum accessions. Specific morphological characteristics were found in C. annuum, C. baccatum, and C.pubescens. Three subclones of each accession were sequenced, and rDNA-ITS polymorphisms were detected in all accessions excluding C. annuum, suggesting that incomplete concerted evolution occurred in rDNA-ITS of Capsicum. The genetic diversity was evaluated using nucleotide polymorphism and diversity. C. annuum had the lowest genetic diversity of all species in this study. The phylogenetic tree formed a species-specific clade for C. annuum, C. baccatum, and C. pubescens. The C. chinense clade existed in the C. frutescens clade, implying that it was a cultivated variant of C. frutescens. C. chacoense likely belonged to the C. baccatum complex according to its morphologic and genetic features. This study indicated that the rDNA-ITS region can be used for simple identification of domesticated Capsicum species.

Genome Size Evolution Mediated by Gypsy Retrotransposons in Brassicaceae

The dynamic activity of transposable elements (TEs) contributes to the vast diversity of genome size and architecture among plants. Here, we examined the genomic distribution and transposition activity of long terminal repeat retrotransposons (LTR-RTs) in Arabidopsis thaliana (Ath) and three of its relatives, Arabidopsis lyrata (Aly), Eutrema salsugineum (Esa), and Schrenkiella parvula (Spa), in Brassicaceae. Our analyses revealed the distinct evolutionary dynamics of Gypsyretrotransposons, which reflects the different patterns of genome size changes of the four species over the past million years. The rate of Gypsy transposition in Aly is approximately five times more rapid than that of Ath and Esa, suggesting an expanding Aly genome. Gypsy insertions in Esa are strictly confined to pericentromeric heterochromatin and associated with dramatic centromere expansion. In contrast, Gypsy insertions in Spa have been largely suppressed over the last million years, likely as a result of a combination of an inherent molecular mechanism of preferential DNA removal and purifying selection at Gypsy elements. Additionally, species-specific clades of Gypsy elements shaped the distinct genome architectures of Aly and Esa.

Extensive cytonuclear discordance in a crested newt from the Balkan Peninsula glacial refugium

Abstract Integration of multilocus data and species distribution modelling into phylogeography allows mitochondrial DNA (mtDNA)-based scenarios to be fine-tuned. We address the question of whether extensive mtDNA substructuring in the crested newt Triturus macedonicus from the Balkan Peninsula is matched in the nuclear genome. We determine the intraspecific population structure based on 52 nuclear DNA markers and project a species distribution model on climate layers for the Last Glacial Maximum. We show that T. macedonicus accumulated nuclear DNA population structure in an area predicted to have been climatically stable during the Pleistocene, with four nuclear DNA groups in the western part of the species range. The distribution of these nuclear DNA groups shows little agreement with that of mtDNA structuring, which shows three highly distinct species-specific clades and a fourth one introgressed from another crested newt species. This cytonuclear discordance conveys that historical biogeographical scenarios based on mtDNA exclusively should be interpreted with caution. Our findings further highlight the important role the Balkan Peninsula has played in the evolution and preservation of European biodiversity.

Molecular diversity of chromosomal metallo-β-lactamase genes in Elizabethkingia genus

Elizabethkingia genus is an opportunistic life-threatening pathogen with an intrinsic multidrug-resistant phenotype. It is the only known microorganism with multi-chromosome-borne metallo-β-lactamase (MBL) genes. To determine the diversity and distribution of MBLs BlaBlaB and BlaGOB in this genus, comprehensive bioinformatic screening was applied in 109 available Elizabethkingia genomes. A total of 23 and 32 novel BlaBlaB and BlaGOB variants were found in Elizabethkingia spp., respectively; 12 and 15 clusters were assigned in these BlaBlaB and BlaGOB based on the amino acid identities and phylogenetic studies. Clustering of some variants did not conform to species-specific clades, which indicated potential inter-species dissemination of MBL genes among Elizabethkingia species. Cloning of representative blaBlaB and blaGOB into E. coli DH5α resulted in increased and diverse minimum inhibitory concentrations (MICs) to most β-lactams, including cephalosporins, carbapenems, and β-lactams-inhibitors. This study extends the database of class B carbapenemases, emphasizing the diversity of different MBL genes in the genus Elizabethkingia, which may represent potential reservoirs of acquired MBLs.

Evolution of the R2R3-MYB gene family in six Rosaceae species and expression in woodland strawberry

R2R3-MYB gene family play important roles in plants development, metabolism, and responses to various biotic and abiotic stresses. In this study, 838 R2R3-MYB genes were identified from six Rosaceae species, including 105 in woodland strawberry (Fragaria vesca), 173 in European pear (Pyrus communis), 219 in apple (Malus domestica), 121 in peach (Prunus persica), 121 in Chinese rose (Rosa chinensis), and 99 in black raspberry (Rubus occidentalis). All R2R3-MYB genes in the six Rosaceae species were clustered into 51 species-specific duplicated clades with 109 genes and 50 lineage-specific duplicated clades with 242 genes according to phylogenetic analysis. R2R3-MYB genes were distributed on all chromosomes in each of the six species, with a small amount of tandem duplication events. The proportion of tandem repeat genes ranged from 0 to 25.1%. The R2R3-MYB protein was conserved in a clade and likely to share similar functions. The distribution of Ks showed the duplication times of R2R3-MYB genes in six Rosaceae species. Furthermore, most of the R2R3-MYB genes had Ka/Ks values less than 1, which indicated they were driven by purifying selection during the evolutionary processes. The GO term enrichment analysis revealed that R2R3-MYB genes in strawberry and black raspberry were more divergent than in other Rosaceae species. Analysis of transcriptomes of 42 different tissues and development stages of woodland strawberry showed that high expression levels of R2R3-MYB suggested that the R2R3-MYB genes in strawberry played a key role in growth and development of both vegetative tissues and fruits. The strawberry R2R3-MYB genes in sub-group of S1, S2, S11, S20, and S22 had high expression levels both in young leaves (YL) and old leaves (OL) strawberry tissues under drought treatments.

NGS-based phylogeny of diphtheria-related pathogenicity factors in different Corynebacterium spp. implies species-specific virulence transmission

BackgroundDiphtheria toxin (DT) is produced by toxigenic strains of the human pathogen Corynebacterium diphtheriae as well as zoonotic C. ulcerans and C. pseudotuberculosis. Toxigenic strains may cause severe respiratory diphtheria, myocarditis, neurological damage or cutaneous diphtheria. The DT encoding tox gene is located in a mobile genomic region and tox variability between C. diphtheriae and C. ulcerans has been postulated based on sequences of a few isolates. In contrast, species-specific sequence analysis of the diphtheria toxin repressor gene (dtxR), occurring both in toxigenic and non-toxigenic Corynebacterium species, has not been done yet. We used whole genome sequencing data from 91 toxigenic and 46 non-toxigenic isolates of different pathogenic Corynebacterium species of animal or human origin to elucidate differences in extracted DT, DtxR and tox-surrounding genetic elements by a phylogenetic analysis in a large sample set.ResultsSequences of both DT and DtxR, extracted from whole genome sequencing data, could be classified in four distinct, nearly species-specific clades, corresponding to C. diphtheriae, C. pseudotuberculosis, C. ulcerans and atypical C. ulcerans from a non-toxigenic toxin gene-bearing wildlife cluster. Average amino acid similarities were above 99% for DT and DtxR within the four groups, but lower between them. For DT, subgroups below species level could be identified, correlating with different tox-comprising mobile genetic elements. In most C. diphtheriae, tox genes were located within known prophages. In contrast, in C. ulcerans diverse tox-including mobile elements could be identified: either prophages differing from C. diphtheriae prophages or an alternative pathogenicity island (PAI) described previously. One isolate showed a different, shorter tox-comprising putative PAI. Beyond the tox-overlapping elements, most isolates harbored a variety of additional prophages.ConclusionOur NGS data from 137 isolates indicate the existence of different genetic backgrounds of DT-mediated pathogenicity in different Corynebacterium species and evolution of once acquired pathogenicity features with the strains. Different groups of pathogenicity-related elements within C. ulcerans imply that tox transmission pathways between isolates may differ in the zoonotic species and contribute to their emerging pathogenic potential.

New Data on the Systematics of Comb-fin Squids Chtenopteryx spp. (Cephalopoda: Chtenopterygidae) from the Canary Islands.

Alejandro Escánez, Álvaro Roura, Rodrigo Riera, Ángel Francisco González, and Ángel Guerra (2018) The systematics of the comb-fin squid species is problematic and poorly resolved. In total, 53 specimens of comb-fin squids (Chtenopteryx spp.) were caught at depths ranging from 30 to 800 m off the Canary Islands (NE Atlantic Ocean). Mantle lengths of the individuals ranged from 18 to 43 mm and the sample included immature, mature male and mature female specimens. Two species of comb-fin squids, Chtenopteryx canariensis and C. sicula, were identified by combining traditional morphological characters with a molecular analysis of a fragment of the cytochrome c oxidase subunit I (COI) gene. Intra- and interspecific genetic distances and maximum likelihood tree analyses based on COI sequences available from GenBank suggest the existence of at least four species, two from the Pacific and two from the Atlantic Ocean. Our data expand the current geographic range of C. canariensis from the NE to NW Atlantic. In the GenBank database, several sequences of comb-fin squid in different species-specific clades have been attributed only to C. sicula, indicating the possible existence of cryptic species and the need to re-analyse these data.

Phylogenetic Relatedness Among Plasmids Harbored by Campylobacter jejuni and Campylobacter coli Isolated From Retail Meats.

Campylobacter jejuni and Campylobacter coli are two of the major causes of foodborne illness. In this study, 29 plasmids isolated from 20 retail meat isolates of Campylobacter jejuni and Campylobacter coli were fully-sequenced individually or as a part of a whole genome sequencing approach. The fully-sequenced plasmids ranged in size from 3 to 119 kb. Molecular characterization of the sequenced plasmids was based on pangenomic analysis and types of genes present on these plasmids and similar ones from GenBank. The plasmids were categorized into four different groups. These groups include type-1 that consisted mainly of pTet plasmids with the tetO gene, type-2 plasmids commonly found in C. coli strains, type-3 which has pVir plasmids, and type-4 that consisted mainly of smaller plasmids. The type-2 plasmids were unique, common among C. coli strains, and carried several conjugative transfer genes. The type-2 plasmids were most similar to a plasmid from Helicobacter pullorum. Maximum parsimony analysis and NeighborNet analysis were used to assess the phylogenetic relatedness among the 29 plasmid sequences presented in this study in addition to the other 104 plasmid sequences of Campylobacter species available in GenBank to date. Results from MP analysis revealed multiple lineages among Campylobacter plasmids which was supported by NeighborNet analysis. Clustering of plasmids did not conform to species-specific clades which suggested an intra-species dissemination of plasmids among Campylobacter species. To our knowledge, this is the first extensive phylogenetic analysis of Campylobacter plasmids sequenced to date.

In silico genome-wide identification, phylogeny and expression analysis of the R2R3-MYB gene family in Medicago truncatula

The R2R3-MYB genes make up one of the largest transcription factor families in plants, and play regulatory roles in various biological processes such as development, metabolism and defense response. Although genome-wide analyses of this gene family have been conducted in several species, R2R3-MYB genes have not been systematically analyzed in Medicago truncatula, a sequenced model legume plant. Here, we performed a comprehensive, genome-wide computational analysis of the structural characteristics, phylogeny, functions and expression patterns of M. truncatula R2R3-MYB genes. DNA binding domains are highly conserved among the 155 putative MtR2R3-MYB proteins that we identified. Chromosomal location analysis revealed that these genes were distributed across all eight chromosomes. Results showed that the expansion of the MtR2R3-MYB family was mainly attributable to segmental duplication and tandem duplication. A comprehensive classification was performed based on phylogenetic analysis of the R2R3-MYB gene families in M. truncatula, Arabidopsis thaliana and other plant species. Evolutionary relationships within clades were supported by clade-specific conserved motifs outside the MYB domain. Species-specific clades have been gained or lost during evolution, resulting in functional divergence. Also, tissue-specific expression patterns were investigated. The functions of stress response-related clades were further verified by the changes in transcript levels of representative R2R3-MYB genes upon treatment with abiotic and biotic stresses. This study is the first report on identification and characterization of R2R3-MYB gene family based on the genome of M. truncatula, and will facilitate functional analysis of this gene family in the future.

Recurrent hybridization and recent origin obscure phylogenetic relationships within the ‘white-headed’ gull (Larus sp.) complex

Species complexes that have undergone recent radiations are often characterized by extensive allele sharing due to recent ancestry and (or) introgressive hybridization. This can result in discordant evolutionary histories of genes and heterogeneous genomes, making delineating species limits difficult. Here we examine the phylogenetic relationships among a complex group of birds, the white-headed gulls (Aves: Laridae), which offer a unique window into the speciation process due to their recent evolutionary history and propensity to hybridize. Relationships were examined among 17 species (61 populations) using a multilocus approach, including mitochondrial and nuclear intron DNA sequences and microsatellite genotype information. Analyses of microsatellite and intron data resulted in some species-based groupings, although most species were not represented by a single cluster. Considerable allele and haplotype sharing among white-headed gull species was observed; no locus contained a species-specific clade. Despite this, our multilocus approach provided better resolution among some species than previous studies. Interestingly, most clades appear to correspond to geographic locality: our BEAST analysis recovered strong support for a northern European/Icelandic clade, a southern European/Russian clade, and a western North American/canus clade, with weak evidence for a high latitude clade spanning North America and northwestern Europe. This geographical structuring is concordant with behavioral observations of pervasive hybridization in areas of secondary contact. The extent of allele and haplotype sharing indicates that ecological and sexual selection are likely not strong enough to complete reproductive isolation within several species in the white-headed gull complex. This suggests that just a few genes are driving the speciation process.

Structural and functional annotation of the MADS-box transcription factor family in grapevine.

BackgroundMADS-box genes encode transcription factors that are involved in developmental control and signal transduction in eukaryotes. In plants, they are associated to numerous development processes most notably those related to reproductive development: flowering induction, specification of inflorescence and flower meristems, establishment of flower organ identity, as well as regulation of fruit, seed and embryo development. Genomic analyses of MADS-box genes in different plant species are providing new relevant information on the function and evolution of this transcriptional factor family. We have performed a true genome-wide analysis of the complete set of MADS-box genes in grapevine (Vitis vinifera), analyzed their expression pattern and establish their phylogenetic relationships (including MIKC* and type I MADS-box) with genes from 16 other plant species. This study was integrated to previous works on the family in grapevine.ResultsA total of 90 MADS-box genes were detected in the grapevine reference genome by completing current gene annotations with a genome-wide analysis based on sequence similarity. We performed a thorough in-depth curation of all gene models and combined the results with gene expression information including RNAseq data to clarifying the expression of newly identified genes and improve their functional characterization. Curated data were uploaded to the ORCAE database for grapevine in the frame of the grapevine genome curation effort. This approach resulted in the identification of 30 additional MADS box genes. Among them, ten new MIKCC genes were identified, including a potential new group of short proteins similar to the SVP protein subfamily. The MIKC* subgroup contains six genes in grapevine that can be grouped in the S (4 genes) and P (2 genes) clades, showing less redundancy than that observed in Arabidopsis thaliana. Expression pattern of these genes in grapevine is compatible with a role in male gametophyte development. Most of the identified new genes belong to the type I MADS-box genes and were classified as members of the Mα and Mγ subclasses. Ours analyses indicate that only few members of type I genes in grapevine have homology in other species and that species-specific clades appeared both in the Mα and Mγ subclasses. On the other hand, as deduced from the phylogenetic analysis with other plant species, genes that can be crucial for development of central cell, endosperm and embryos seems to be conserved in plants.ConclusionsThe genome analysis of MADS-box genes in grapevine, the characterization of their pattern of expression and the phylogenetic analysis with other plant species allowed the identification of new MADS-box genes not yet described in other plant species as well as basic characterization of their possible role, particularly in the case of type I and MIKC* genes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2398-7) contains supplementary material, which is available to authorized users.