The complete mitochondrial genome of <Emphasis Type="Italic">Vanessa indica</Emphasis> and phylogenetic analyses of the family Nymphalidae

Vanessa indica is a small butterfly lacking historical molecular and biological research. Vanessa indica belongs to the family Nymphalidae (Lepidoptera: Papilionoidea), which is the largest group of butterflies and are nearly ubiquitous. However, after more than a century of taxonomic and molecular studies, there is no consensus for family classification, and the phylogenetic relationships within Nymphalidae are controversial. The first objective was to sequence and characterize the complete mitochondrial genome of V. indica. The most important objective was to completely reconstruct the phylogenetic relationships for family members within Nymphalidae. The mitochondrial genomic DNA (mtDNA) of V. indica was extracted and amplified by polymerase chain reaction. The complete mitochondrial sequence was annotated and characterized by analyzing sequences with SeqMan program. The phylogenetic analyses were conducted on thirteen protein coding genes (PCGs) in 95 mtDNA of Nymphalidae downloaded from GenBank for reference using the maximum likelihood method and Bayesian inference to ensure the validity of the results. The complete mitogenome was a circular molecule with 15,191 bp consisting of 13 protein coding genes, two ribosomal RNA genes (16S rRNA and 12S rRNA), 22 transfer RNA (tRNA) genes, and an A?+?T-rich region (D-loop). The nucleotide composition of the genome was highly biased for A?+?T content, which accounts for 80.0% of the nucleotides. All the tRNAs have putative secondary structures that are characteristic of mitochondrial tRNAs, except tRNA^Ser(AGN). All the PCGs started with ATN codons, except cytochrome c oxidase subunit 1 (COX1), which was found to start with an unusual CGA codon. Four genes were observed to have unusual codons: COX1 terminated with atypical TT and the other three genes terminated with a single T. The A?+?T rich region of 327 bp consisted of repetitive sequences, including a ATAGA motif, a 19-bp poly-T stretch, and two microsatellite-like regions (TA)₈. The phylogenetic analyses consistently placed Biblidinae as a sister cluster to Heliconiinae and Calinaginae as a sister clade to Satyrinae. Moreover, the phylogenetic tree identified Libytheinae as a monophyletic group within Nymphalidae. The complete mitogenome of V. indica was 15,191 bp with mitochondrial characterizations common for lepidopteran species, which enriched the mitochondria data of Nymphalid species. And the phylogenetic analysis revealed different classifications and relationships than those previously described. Our results are significant because they would be useful in further understanding of the evolutionary biology of Nymphalidae.     

Identification and polymorphism of pectinase genes <Emphasis Type="Italic">PGU</Emphasis> in the <Emphasis Type="Italic">Saccharomyces bayanus</Emphasis> complex

Pectinase (endo-polygalacturonase) is the key enzyme splitting plant pectin. The corresponding single gene PGU1 is documented for the yeast S. cerevisiae. On the basis of phylogenetic analysis of the PGU nucleotide sequence available in the GenBank, a family of divergent PGU genes is found in the species complex S. bayanus: S. bayanus var. uvarum, S. eubayanus, and hybrid taxon S. pastorianus. The PGU genes have different chromosome localization.     

The complete genome sequence and phylogenetic analysis of the mitochondrial DNA of the wood-decaying fungus <Emphasis Type="Italic">Fomitopsis palustris</Emphasis>

The complete mitochondrial genome sequence of the wood-decaying fungus Fomitopsis palustris (Basidiomycete, Agaricomycotina) was determined by next-generation sequencing technology. The complete sequence of the circular mitochondrial DNA of F. palustris was 63,479 bp in length with a 75.98% AT content. The mitochondrial genome encoded 14 conserved proteins, 2 ribosomal RNAs, 26 transfer RNAs, and 19 additional open reading frames. The coxI and cob genes contained six and one group I introns, respectively, and encoded eight open reading frames, including seven intron-encoded endonucleases. The complete mitochondrial genome of F. palustris presented herein represents the first such report for brown rot basidiomycetes. In addition, the BLAST score ratio and phylogenetic analysis may open new avenues to understanding the evolutionary status of this fungus .     

Horizontal transfer of the C-termini of <Emphasis Type="Italic">tccC</Emphasis> genes in <Emphasis Type="Italic">Photorhabdus</Emphasis> and <Emphasis Type="Italic">Xenorhabdus</Emphasis>

The high molecular weight insecticidal toxin complexes (Tcs), including four toxin-complex loci (tca, tcb, tcc and tcd), were first identified in Photorhabdus luminescens W14. Each member of tca, tcb or tcc is required for oral toxicity of Tcs. However, the sequence sources of the C-termini of tccC3, tccC4, tccC6 and tccC7 are unknown. Here, we performed a whole genome survey to identify the orthologs of Tc genes, and found 165 such genes in 14 bacterial genomes, including 40 genes homologous to tccC1-7 in P. luminescens TT01. The sequence sources of the C-termini of tccC2-6 were determined by sequence analysis. Further phylogenetic investigations suggested that the C-termini of 6 tccC genes experienced horizontal gene transfer events.     

Complete mitochondrial genome sequence of <Emphasis Type="Italic">Cucullaea labiata</Emphasis> (Arcoida: Cucullaeidae) and phylogenetic implications

The complete mitochondrial genome of Cucullaea labiata (Arcoida: Cucullaeidae) was firstly determined in this study in order to better understand the phylogenetic relationship between Cucullaeidae and Arcidae. The C. labiata mitochondrial genome was 25,845 bp in size and contained 12 protein-coding genes, 2 rRNA and 22 tRNA genes. The number and the location of the tRNA genes were different from three Arcidae species (Scapharca broughtonii, Scapharca kagoshimensis and Tegillarca granosa). Gene arrangement also differed dramatically. The length of the non-coding regions was 10,559 bp, in which the largest one (6057 bp) included eight point nine copies of a 659 bp repeat motif. The number of repeated sequences was different in different individuals, similar to the findings from the mitochondrial genome of S. broughtonii and Placopecten magellanicus. One intron was found in cox1 gene both in CL_98 and in CL_99 individuals of C. labiata. The reason why mitochondrial introns are retained so scarcely in bivalve taxa needs further research. Phylogenetic analyses based on 12 concatenated amino acid sequences of protein-coding genes supported Cucullaeidae was the sister group of Arcidae.     

Genome-wide survey and expression analysis of the stress-associated protein gene family in desert poplar, <Emphasis Type="Italic">Populus euphratica</Emphasis>

Stress-associated proteins (SAPs) are a novel class of zinc finger proteins that extensively participate in abiotic stress responses. To date, no overall analysis and expression profiling of SAP genes in woody plants have been reported. Populus euphratica is distributed in desert regions and is extraordinarily adaptable to abiotic stresses. Thus, it is regarded as a promising candidate for studying abiotic stress resistance mechanisms of woody plants. In this study, 18 non-redundant SAP genes were identified from the genome of P. euphratica using basic local alignment search tool algorithms and functional domain verification. Among these 18 PeuSAP genes, 15 were intronless. To investigate the evolutionary relationships of SAP genes in P. euphratica and other Salicaceae plants, phylogenetic analyses were performed. Subsequently, the expression profiles of the 18 PeuSAP genes were analyzed in different tissues and under various stresses (drought, salt, heat, cold, and abscisic acid (ABA) treatment) using quantitative real-time PCR. Tissue expression analysis indicated that PeuSAPs showed no tissue specificity. PeuSAPs were induced by multiple abiotic stresses, especially drought, salt, and heat stresses, perhaps because of abundant cis-acting heat shock elements and drought-inducible elements in the promoter regions of the PeuSAPs. Moreover, single nucleotide polymorphisms (SNPs) variant analysis revealed many synonymous and non-synonymous SNPs in PeuSAP genes, but the zinc finger structure was conserved during evolution. These results provide an overview of the SAP gene family in P. euphratica and a reference for further functional research on PeuSAP genes.     

A new primer set for amplification of the cytochrome <Emphasis Type="Italic">b</Emphasis> gene in lantern fishes (Myctophidae)

New universal primers are offered for amplification of complete sequence of mitochondrial cyt b gene in lantern fishes of the family Myctophidae. Analysis of cyt b variability in the species of seven genera (Bolinichthys, Ceratoscopelus, Diaphus, Diogenichthys, Lampanyctus, Lepidophanes, Nannobrachium) demonstrates considerable divergence between species: an average of 18.2% (p-distances). Diversity (nucleotide diversity, number of segregating and phylogenetically informative sites, average number of nucleotide differences) and divergence significantly exceed those of another widely used mitochondrial marker, a fragment of cytochrome c oxidase subunit 1 sequence (cox 1). The mitochondrial cyt b gene amplified with the developed primers can be recommended as an informative tool for phylogenetic and population studies of lantern fishes.     

Characterization and phylogenetic analysis of the complete mitogenome of a rare cavefish, <Emphasis Type="Italic">Sinocyclocheilus multipunctatus</Emphasis> (Cypriniformes: Cyprinidae)

The genus Sinocyclocheilus is a representative group of cave creatures. However, genetic studies on Sinocyclocheilus are rare. The primary objective of this study was to explore the structure and feature of the complete mitochondrial genome of S. multipunctatus, and reconstruct the mitogenomic phylogeny of Sinocyclocheilus. The mitochondrial DNA of S. multipunctatus was amplified by overlapping PCR fragments. The mitogenome was assembled by the SeqMan and annotated using MitoAnnotator. The phylogenetic tree was established using the Bayesian inference and Maximum likelihood methods. The mitogenome of S. multipunctatus is a typical circular molecule of 16,586 bp with base composition A (31.25%), T (25.90%), G (16.35%), and C (26.50%), and consists of 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) genes, two ribosomal RNAs, and a 931 bp control region. Phylogenetic analysis reveals two clades in the Sinocyclocheilus with robust support. S. multipunctatus is close to a newly discovered cavefish, S. ronganensis. We obtained and described the complete mitogenome of S. multipunctatus, and investigated its phylogenetic status, which may provide a valuable resource for future phylogenetic analyses and population genetic studies in Sinocyclocheilus.     

Association between expression levels and growth trait-related SNPs located in promoters of the <Emphasis Type="Italic">MC4R</Emphasis> and <Emphasis Type="Italic">MSTN</Emphasis> genes in <Emphasis Type="Italic">Spinibarbus hollandi</Emphasis>

Melanocortin 4 receptor: (MC4R) and Myostatin (MSTN) are two important growth trait-related genes in animals. In this study, we showed that two SNPs, MC4R-719A>G and MSTN-519C>T, found in the promoters of the MC4R and MSTN genes, respectively, are both associated with growth traits in Spinibarbus hollandi. Furthermore, we observed that there were significant associations between the expression levels of the MC4R and MSTN genes and these two growth trait-related SNPs. The expression level of MC4R gene in brain was lower in GG genotype fish with extremely high growth performance than that in AA genotype fish with extremely low growth performance. Expression level of the MSTN gene in muscle was lower in TT genotype fish with extremely high growth performance than that in CC and CT genotype fish with lower growth performance. The results indicated that these SNPs located in the promoters of MC4R and MSTN are associated with growth-related traits through modification of gene expression levels. The MSTN and MC4R SNPs may have useful application in effective marker-assisted selection aimed to increase output in S. hollandi.     

Identification of <Emphasis Type="Italic">Bacillus</Emphasis> Probiotics Isolated from Soil Rhizosphere Using 16S rRNA, <Emphasis Type="Italic">recA</Emphasis>, <Emphasis Type="Italic">rpoB</Emphasis> Gene Sequencing and RAPD-PCR

Some Bacillus species, especially Bacillus subtilis and Bacillus pumilus groups, have highly similar 16S rRNA gene sequences, which are hard to identify based on 16S rDNA sequence analysis. To conquer this drawback, rpoB, recA sequence analysis along with randomly amplified polymorphic (RAPD) fingerprinting was examined as an alternative method for differentiating Bacillus species. The 16S rRNA, rpoB and recA genes were amplified via a polymerase chain reaction using their specific primers. The resulted PCR amplicons were sequenced, and phylogenetic analysis was employed by MEGA 6 software. Identification based on 16S rRNA gene sequencing was underpinned by rpoB and recA gene sequencing as well as RAPD-PCR technique. Subsequently, concatenation and phylogenetic analysis showed that extent of diversity and similarity were better obtained by rpoB and recA primers, which are also reinforced by RAPD-PCR methods. However, in one case, these approaches failed to identify one isolate, which in combination with the phenotypical method offsets this issue. Overall, RAPD fingerprinting, rpoB and recA along with concatenated genes sequence analysis discriminated closely related Bacillus species, which highlights the significance of the multigenic method in more precisely distinguishing Bacillus strains. This research emphasizes the benefit of RAPD fingerprinting, rpoB and recA sequence analysis superior to 16S rRNA gene sequence analysis for suitable and effective identification of Bacillus species as recommended for probiotic products.     

Genome-wide association study of <Emphasis Type="Italic">Mycobacterium avium</Emphasis> subspecies <Emphasis Type="Italic">Paratuberculosis</Emphasis> infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50?K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50?K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P?<?5?×?10^??5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1?Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5?Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

     

Identification of Japanese <Emphasis Type="Italic">Lymantria</Emphasis> species (Lepidoptera: Lymantriidae) based on DNA sequences

We examined the utility of the cytochrome c oxidase subunit I (COI) and II (COII) genes of mitochondrial DNA (mtDNA) as a tool to identify nine Japanese Lymantria species, including four Asian gypsy moth species [Lymantria dispar japonica (Motschulsky), Lymantria umbrosa (Butler), Lymantria albescens Hori and Umeno, and Lymantria postalba Inoue]. In phylogenetic trees constructed for the COI and COII genes using the maximum likelihood methods, we could identify seven out of the nine Lymantria species (L. albescens and L. postalba could not be identified). These results suggest that the DNA sequences of the COI and COII genes are useful for identifying Japanese Lymantria species.     

Molecular Phylogeny of Serotines (Mammalia,Chiroptera, <Emphasis Type="Italic">Eptesicus</Emphasis>): Evolutionary and Taxonomical Aspects of the <Emphasis Type="Italic">E. serotinus</Emphasis> Species Group

This article provides a phylogenetic analysis of five nuclear and mitochondrial cytochrome b genes of palaearctic serotines. Nuclear data yield five monophyletic clades: Botta’s serotine and the South African long-tailed house bat E. hottentottus; the common serotine bat (including all studied E. serotinus subspecies and andersoni form of undefined status) and the meridional serotine E. isabellinus; the Gobi serotine, including Bobrinski’s serotine; the northern bat; and New World serotines. We found latest taxonomic decisions regarding mirza and pachyomus questionable and needing further revision. The significant inconsistency between mitochondrial and nuclear phylogenies obtained for genes of different inheritance systems suggests repetitive introgression events in the evolution of the genus.     

Nuclear mtDNA pseudogenes as a source of new variants of the mtDNA cytochrome <Emphasis Type="Italic">b</Emphasis> haplotypes: A case study of Siberian rubythroat <Emphasis Type="Italic">Luscinia calliope</Emphasis> (Muscicapidae,Aves)

Sequence polymorphism of the mitochondrial DNA cytochrome b gene fragment was analyzed in 21 specimens of subspecies Luscinia calliope calliope (Pallas, 1776) and two specimens of L. c. anadyrensis (Portenko, 1939). On sequence chromatograms, in 19 specimens of L. c. calliope, double peaks of heteroplasmy type in the taxon-specific positions were revealed. Moreover, two clone variants were identified. The first variant was the calliope mitochondrial cyt b gene and the second was the nuclear cyt b pseudogene, similar to the mitochondrial haplotype anadyrensis-camtschatkensis. In L. c. anadyrensis, four clone variants, represented by the mitochondrial calliope and anadyrensis-camtschatkensis cyt b genes and nuclear calliope and sachalinensis cyt b pseudogenes, were identified. Some nuclear cyt b pseudogenes were highly similar (98–99%) to the mitochondrial genes of the subspecies L. c. anadyrensis, L. c. camtschatkensis, and L. c. sachalinensis. In the same time, the majority of nuclear pseudogene sequences were characterized by a high level of polymorphism, caused by nonsynonymous substitutions (up to five substitutions per sequence), the presence of indels in some of the clones, and TAA and TGA stop codons. In our opinion, the mitochondrial haplotypes anadyrensis-camtschatkensis and sachalinensis occurred as a result of intergenomic homologous recombination. This finding provides a new insight into the colonization history of the northeastern part of the range by L. calliope, according to which populating the territory of Chukotka, Kamchatka, and Sakhalin took place at different times and along the independent pathways.     

Molecular diversity and phylogeny of indigenous <Emphasis Type="Italic">Rhizobium leguminosarum</Emphasis> strains associated with <Emphasis Type="Italic">Trifolium repens</Emphasis> plants in Romania

The symbiotic nitrogen fixing legumes play an essential role in sustainable agriculture. White clover (Trifolium repens L.) is one of the most valuable perennial legumes in pastures and meadows of temperate regions. Despite its great agriculture and economic importance, there is no detailed available information on phylogenetic assignation and characterization of rhizobia associated with native white clover plants in South-Eastern Europe. In the present work, the diversity of indigenous white clover rhizobia originating in 11 different natural ecosystems in North-Eastern Romania were assessed by a polyphasic approach. Initial grouping showed that, 73 rhizobial isolates, representing seven distinct phenons were distributed into 12 genotypes, indicating a wide phenotypic and genotypic diversity among the isolates. To clarify their phylogeny, 44 representative strains were used in sequence analysis of 16S rRNA gene and IGS fragments, three housekeeping genes (atpD, glnII and recA) and two symbiosis-related genes (nodA and nifH). Multilocus sequence analysis (MLSA) phylogeny based on concatenated housekeeping genes delineated the clover isolates into five putative genospecies. Despite their diverse chromosomal backgrounds, test strains shared highly similar symbiotic genes closely related to Rhizobium leguminosarum biovar trifolii. Phylogenies inferred from housekeeping genes were incongruent with those of symbiotic genes, probably due to occurrence of lateral transfer events among native strains. This is the first polyphasic taxonomic study to report on the MLSA-based phylogenetic diversity of indigenous rhizobia nodulating white clover plants grown in various soil types in South-Eastern Europe. Our results provide valuable taxonomic data on native clover rhizobia and may increase the pool of genetic material to be used as biofertilizers.     

Genetic relationships of Chukchi charr <Emphasis Type="Italic">Salvelinus andriashevi</Emphasis> and Taranetz charr <Emphasis Type="Italic">Salvelinus taranetzi</Emphasis>

On the basis of sequence analysis of the mitochondrial DNA (mtDNA) control region (CR), cytochrome b (Cytb), and cytochrome oxidase-1 (CoI) genes, the relationships of endemic species Salvelinus andriashevi Berg, 1948, represented by the only population from Lake Estikhed (Chukotka), were estimated. The data on the genealogical analysis of mtDNA haplotypes supported phylogenetic closeness of S. andriashevi and S. taranetzi. It was also demonstrated that the specimens of Chukchi charr, along with Salvelinus sp. 4 (Lake Nachikinskoe), S. krogiusae (Lake Dal’nee), S. boganidae and S. elgyticus (Lake El’gygytgyn), and S. a. erythrinus from Canada’s Northwest Territories (NWT) belonged to the Arctic group of Taranetz charr. The problem of coincidence of taxonomic differentiation of charrs of the genus Salvelinus based on morphological and genetic analyses is discussed.     

Multilocus sequence typing for phylogenetic view and <Emphasis Type="Italic">vip</Emphasis> gene diversity of <Emphasis Type="Italic">Bacillus thuringiensis</Emphasis> strains of the Assam soil of North East India

An agriculturally important insecticidal bacterium, Bacillus thuringiensis have been isolated from the soil samples of various part of Assam including the Kaziranga National Park. Previously, the isolates were characterized based on morphology, 16S rDNA sequencing, and the presence of the various classes’ crystal protein gene(s). In the present study, the phylogenetic analysis of a few selected isolates was performed by an unambiguous and quick method called the multiple locus sequence typing (MLST). A known B. thuringiensis strain kurstaki 4D4 have been used as a reference strain for MLST. A total of four the MLST locus of housekeeping genes, recF, sucC, gdpD and yhfL were selected. A total of 14 unique sequence types (STs) was identified. A total number of alleles identified for the locus gdpD and sucC was 12, followed by locus yhfL was 11, however, only 6 alleles were detected for the locus recF. The phylogenetic analysis using MEGA 7.0.26 showed three major lineages. Approximately, 87% of the isolates belonged to the STs corresponding to B. thuringiensis, whereas two isolates, BA07 and BA39, were clustered to B. cereus. The isolates were also screened for the diversity of vegetative insecticidal protein (vip) genes. In all, 8 isolates showed the presence of vip1, followed by 7 isolates having vip2 and 6 isolates for vip3 genes. The expression of Vip3A proteins was analyzed by western blot analyses and expression of the Vip3A protein was observed in the isolate BA20. Thus, the phylogenetic relationship and diversity of Bt isolates from Assam soil was established based on MLST, in addition, found isolates having vip genes, which could be used for crop improvement.