Drosophila melanogaster Transcriptome Response to Different Wolbachia Strains.

Wolbachia is a maternally inherited, intercellular bacterial symbiont of insects and some other invertebrates. Here, we investigated the effect of two different Wolbachia strains, differing in a large chromosomal inversion, on the differential expression of genes in D. melanogaster females. We revealed significant changes in the transcriptome of the infected flies compared to the uninfected ones, as well as in the transcriptome of flies infected with the Wolbachia strain, wMelPlus, compared to flies infected with the wMelCS112 strain. We linked differentially expressed genes (DEGs) from two pairwise comparisons, "uninfected-wMelPlus-infected" and "uninfected-wMelCS112-infected", into two gene networks, in which the following functional groups were designated: "Proteolysis", "Carbohydrate transport and metabolism", "Oxidation-reduction process", "Embryogenesis", "Transmembrane transport", "Response to stress" and "Alkaline phosphatases". Our data emphasized similarities and differences between infections by different strains under study: a wMelPlus infection results in more than double the number of upregulated DEGs and half the number of downregulated DEGs compared to a wMelCS112 infection. Thus, we demonstrated that Wolbachia made a significant contribution to differential expression of host genes and that the bacterial genotype plays a vital role in establishing the character of this contribution.

Population Structure and Genetic Diversity of the 175 Soybean Breeding Lines and Varieties Cultivated in West Siberia and Other Regions of Russia.

Soybean is a leguminous plant cultivated in many countries and is considered important in the food industry due to the high levels of oil and protein content in the beans. The high demand for soybeans and its products in the industry requires the expansion of cultivation areas. Despite climatic restrictions, West Siberia is gradually expanding its area of soybean cultivation. In this study, we present the first analysis of the population structure and genetic diversity of the 175 soybean Glycine max breeding lines and varieties cultivated in West Siberia (103 accessions) and other regions of Russia (72 accessions), and we compare them with the cultivated soybean varieties from other geographical locations. Principal component analysis revealed several genetic clusters with different levels of genetic heterogeneity. Studied accessions are genetically similar to varieties from China, Japan, and the USA and are genetically distant to varieties from South Korea. Admixture analysis revealed four ancestry groups based on genetic ancestry and geographical origin, which are consistent with the regions of cultivation and origin of accessions and correspond to the principal component analysis result. Population statistics, including nucleotide diversity, Tajima's D, and linkage disequilibrium, are comparatively similar to those observed for studied accessions of a different origin. This study provides essential population and genetic information about the unique collection of breeding lines and varieties cultivated in West Siberia and other Russian regions to foster further evolutionary, genome-wide associations and functional breeding studies.

WAX INDUCER 1 Regulates ß-Diketone Biosynthesis by Mediating Expression of the Cer-cqu Gene Cluster in Barley.

Plant surface properties are crucial determinants of resilience to abiotic and biotic stresses. The outer layer of the plant cuticle consists of chemically diverse epicuticular waxes. The WAX INDUCER1/SHINE subfamily of APETALA2/ETHYLENE RESPONSIVE FACTORS regulates cuticle properties in plants. In this study, four barley genes homologous to the Arabidopsis thaliana AtWIN1 gene were mutated using RNA-guided Cas9 endonuclease. Mutations in one of them, the HvWIN1 gene, caused a recessive glossy sheath phenotype associated with ß-diketone deficiency. A complementation test for win1 knockout (KO) and cer-x mutants showed that Cer-X and WIN1 are allelic variants of the same genomic locus. A comparison of the transcriptome from leaf sheaths of win1 KO and wild-type plants revealed a specific and strong downregulation of a large gene cluster residing at the previously known Cer-cqu locus. Our findings allowed us to postulate that the WIN1 transcription factor in barley is a master mediator of the ß-diketone biosynthesis pathway acting through developmental stage- and organ-specific transactivation of the Cer-cqu gene cluster.

Analysis of Genome Structure and Its Variations in Potato Cultivars Grown in Russia.

Solanum tuberosum L. (common potato) is one of the most important crops produced almost all over the world. Genomic sequences of potato opens the way for studying the molecular variations related to diversification. We performed a reconstruction of genomic sequences for 15 tetraploid potato cultivars grown in Russia using short reads. Protein-coding genes were identified; conserved and variable parts of pan-genome and the repertoire of the NBS-LRR genes were characterized. For comparison, we used additional genomic sequences for twelve South American potato accessions, performed analysis of genetic diversity, and identified the copy number variations (CNVs) in two these groups of potato. Genomes of Russian potato cultivars were more homogeneous by CNV characteristics and have smaller maximum deletion size in comparison with South American ones. Genes with different CNV occurrences in two these groups of potato accessions were identified. We revealed genes of immune/abiotic stress response, transport and five genes related to tuberization and photoperiod control among them. Four genes related to tuberization and photoperiod were investigated in potatoes previously (phytochrome A among them). A novel gene, homologous to the poly(ADP-ribose) glycohydrolase (PARG) of Arabidopsis, was identified that may be involved in circadian rhythm control and contribute to the acclimatization processes of Russian potato cultivars.

New Wolbachia pipientis Genotype Increasing Heat Stress Resistance of Drosophila melanogaster Host Is Characterized by a Large Chromosomal Inversion.

The maternally transmitted endocellular bacteria Wolbachia is a well-known symbiont of insects, demonstrating both negative and positive effects on host fitness. The previously found Wolbachia strain wMelPlus is characterized by a positive effect on the stress-resistance of its host Drosophila melanogaster, under heat stress conditions. This investigation is dedicated to studying the genomic underpinnings of such an effect. We sequenced two closely related Wolbachia strains, wMelPlus and wMelCS112, assembled their complete genomes, and performed comparative genomic analysis engaging available Wolbachia genomes from the wMel and wMelCS groups. Despite the two strains under study sharing very close gene-composition, we discovered a large (>1/6 of total genome) chromosomal inversion in wMelPlus, spanning through the region that includes the area of the inversion earlier found in the wMel group of Wolbachia genotypes. A number of genes in unique inversion blocks of wMelPlus were identified that might be involved in the induction of a stress-resistant phenotype in the host. We hypothesize that such an inversion could rearrange established genetic regulatory-networks, causing the observed effects of such a complex fly phenotype as a modulation of heat stress resistance. Based on our findings, we propose that wMelPlus be distinguished as a separate genotype of the wMelCS group, named wMelCS3.

Easy and Effective Method for Extracting and Purifying Wolbachia Genomic DNA.

A number of methods for extracting the DNA of maternally inherited obligate intracellular bacteria Wolbachia from an insect host and its subsequent purification have been described in previous scholarship. As Wolbachia is present in the hosts' organisms in rather low quantities, these techniques used to be quite labor-intensive. For this paper, we analyzed them in detail, searched for a possibility to simplify and accelerate the protocol, and proposed an easy and effective method for isolating Wolbachia DNA from Drosophila melanogaster with a purity sufficient for genomic sequencing. Our method involves the centrifugation of homogenized flies or just their ovaries, as the most Wolbachia-enriched tissue, followed by the filtration of homogenate and extraction of DNA using a modified version of the Livak buffer protocol. The proportion of Wolbachia DNA in the total DNA was quantified based on the results of sequencing with the use of the Illumina MiSeq platform and a pipeline of bioinformatic analysis. For the two analyzed D. melanogaster lines infected with two different Wolbachia strains, the proportion was at least 68 and 94%, respectively.

Androgen Receptor Gene CAG Repeat Length Varies and Affects Semen Quality in an Ethnic-Specific Fashion in Young Men from Russia.

Male infertility is a multi-factorial and multi-genetic disorder, and the prevalence of male infertility in the world is estimated at 5-35%. The search for the causes of male infertility allowed for identifying a number of genetic factors including a single X-linked gene of the androgen receptor (AR), and some of its alleles are assumed to negatively affect male fertility. Our aim was (1) to study the variability of the length of CAG repeats of the AR gene and possible associations in the AR CAG genetic variants with semen quality and reproductive hormone levels in a population-based cohort of men and (2) to estimate distributions of AR CAG repeat alleles and associations with semen parameters in different ethnic subgroups. The cohort of 1324 young male volunteers of different ethnicities (median age 23.0 years) was recruited from the general population of five cities of the Russian Federation, regardless of their fertility status. Semen quality (sperm concentration, motility and morphology), reproductive hormone levels (testosterone, estradiol, LH, FSH and inhibin B) and trinucleotide (CAG) n repeat polymorphism in exon 1 of the AR gene were evaluated. The semen samples were analyzed according to the WHO laboratory manual (WHO, 2010), serum hormones were measured by enzyme immunoassay, and the AR CAG repeat length was analyzed by direct sequencing of leukocyte DNA. The median AR CAG repeat length in men of our multi-ethnic population was 23 (range 6-39). In the entire study population, a significant difference (p ≤ 0.05) was found in the frequency distribution and the mean values for the CAG repeat length between the groups with normal (23.2 ± 3.3) and impaired semen quality (23.9 ± 3.2). Additionally, we demonstrated that the total sperm count, sperm concentration, progressive motility and normal morphology were lower in the category of long CAG repeats (CAG ≥ 25) compared with those in the category of short CAG repeats (CAG ≤ 19); however, hormonal parameters did not differ between the long and short CAG categories, with the exception of estradiol. Significant differences were observed in the AR CAG repeat length between the most common ethnic cohorts of Slavs (Caucasians), Buryats (Asians), and Yakuts (Asians). The Buryats and Yakuts had a higher number of CAG repeats than the Slavs (medians: Slavs-23; Buryats-24; Yakuts-25). The range of alleles differed among ethnicities, with the Slavs having the largest range (7-36 repeats, 24 alleles total), the Yakuts having the smallest range (18-32 repeats, 14 alleles total) and the Buryats having the middle range (11-39 repeats, 20 alleles total). The longer CAG repeats were associated with an impaired semen quality within the Slavic (CAG ≥ 25) and Buryat (CAG ≥ 28) groups, but this effect was not found in Yakuts. Hormonal parameters did not differ between the three CAG repeat categories in men of all ethnic groups. This is the largest Russian study of the distribution of AR CAG repeats and the search for association between length of AR CAG repeat tract and impaired spermatogenesis in men from the general population. Our results confirmed the association of longer CAG repeats with a risk of impaired semen quality, but this association can be modified by ethnic origin. Identification of the number of AR CAG repeats can be an effective tool to assess the risk of male subfertility and the control of androgen hormone therapy of reproductive diseases.

Features of Activity of the Phenylpropanoid Biosynthesis Pathway in Melanin-Accumulating Barley Grains.

Barley (Hordeum vulgare L.) grain pigmentation is caused by two types of phenolic compounds: anthocyanins (which are flavonoids) give a blue or purple color, and melanins (which are products of enzymatic oxidation and polymerization of phenolic compounds) give a black or brown color. Genes Ant1 and Ant2 determine the synthesis of purple anthocyanins in the grain pericarp, whereas melanins are formed under the control of the Blp1 gene in hulls and pericarp tissues. Unlike anthocyanin synthesis, melanin synthesis is poorly understood. The objective of the current work was to reveal features of the phenylpropanoid biosynthesis pathway functioning in melanin-accumulating barley grains. For this purpose, comparative transcriptomic and metabolomic analyses of three barley near-isogenic lines accumulating anthocyanins, melanins, or both in the grain, were performed. A comparative analysis of mRNA libraries constructed for three stages of spike development (booting, late milk, and early dough) showed transcriptional activation of genes encoding enzymes of the general phenylpropanoid pathway in all the lines regardless of pigmentation; however, as the spike matured, unique transcriptomic patterns associated with melanin and anthocyanin synthesis stood out. Secondary activation of transcription of the genes encoding enzymes of the general phenylpropanoid pathway together with genes of monolignol synthesis was revealed in the line accumulating only melanin. This pattern differs from the one observed in the anthocyanin-accumulating lines, where - together with the genes of general phenylpropanoid and monolignol synthesis pathways - flavonoid biosynthesis genes were found to be upregulated, with earlier activation of these genes in the line accumulating both types of pigments. These transcriptomic shifts may underlie the observed differences in concentrations of phenylpropanoid metabolites analyzed in the grain at a late developmental stage by high-performance liquid chromatography. Both melanin-accumulating lines showed an increased total level of benzoic acids. By contrast, anthocyanin-accumulating lines showed higher concentrations of flavonoids and p-coumaric and ferulic acids. A possible negative effect of melanogenesis on the total flavonoid content and a positive influence on the anthocyanin content were noted in the line accumulating both types of pigments. As a conclusion, redirection of metabolic fluxes in the phenylpropanoid biosynthesis pathway occurs when melanin is synthesized.

The Expression of miRNAs Involved in Long-Term Memory Formation in the CNS of the Mollusk Helix lucorum.

Mollusks are unique animals with a relatively simple central nervous system (CNS) containing giant neurons with identified functions. With such simple CNS, mollusks yet display sufficiently complex behavior, thus ideal for various studies of behavioral processes, including long-term memory (LTM) formation. For our research, we use the formation of the fear avoidance reflex in the terrestrial mollusk Helix lucorum as a learning model. We have shown previously that LTM formation in Helix requires epigenetic modifications of histones leading to both activation and inactivation of the specific genes. It is known that microRNAs (miRNAs) negatively regulate the expression of genes; however, the role of miRNAs in behavioral regulation has been poorly investigated. Currently, there is no miRNAs sequencing data being published on Helix lucorum, which makes it impossible to investigate the role of miRNAs in the memory formation of this mollusk. In this study, we have performed sequencing and comparative bioinformatics analysis of the miRNAs from the CNS of Helix lucorum. We have identified 95 different microRNAs, including microRNAs belonging to the MIR-9, MIR-10, MIR-22, MIR-124, MIR-137, and MIR-153 families, known to be involved in various CNS processes of vertebrates and other species, particularly, in the fear behavior and LTM. We have shown that in the CNS of Helix lucorum MIR-10 family (26 miRNAs) is the most representative one, including Hlu-Mir-10-S5-5p and Hlu-Mir-10-S9-5p as top hits. Moreover, we have shown the involvement of the MIR-10 family in LTM formation in Helix. The expression of 17 representatives of MIR-10 differentially changes during different periods of LTM consolidation in the CNS of Helix. In addition, using comparative analysis of microRNA expression upon learning in normal snails and snails with deficient learning abilities with dysfunction of the serotonergic system, we identified a number of microRNAs from several families, including MIR-10, which expression changes only in normal animals. The obtained data can be used for further fundamental and applied behavioral research.

Fitness Analysis and Transcriptome Profiling Following Repeated Mild Heat Stress of Varying Frequency in Drosophilamelanogaster Females.

Understanding how repeated stress affects metabolic and physiological functions in the long run is of crucial importance for evaluating anthropogenic pressure on the environment. We investigated fertility, longevity and metabolism in D. melanogaster females exposed to short-term heat stress (38 °C, 1 h) repeated daily or weekly. Daily stress was shown to cause a significant decrease in both fertility and longevity, as well as in body mass and triglyceride (fat) content, but a significant increase in trehalose and glucose content. Weekly stress did not affect longevity and carbohydrate metabolism but resulted in a significant decrease in body mass and fat content. Weekly stress did not affect the total level of fertility, despite sharp fertility drops on the exact days of stressing. However, stressing insects weekly, only in the first two weeks after eclosion, caused a significant increase in the total level of fertility. The analysis of differentially expressed genes in the fat bodies and adjacent tissues of researched groups with the use of RNA-Seq profiling revealed changes in signal pathways related to proteolysis/digestion, heat shock protein 23, and in the tightly linked stress-inducible humoral factor Turandot gene network.

E. coli aggregation and impaired cell division after terahertz irradiation.

In this study we demonstrated that exposure of Escherichia coli (E. coli) to terahertz (THz) radiation resulted in a change in the activities of the tdcABCDEFGR and matA-F genes (signs of cell aggregation), gene yjjQ (signs of suppression of cell motility), dicABCF, FtsZ, and minCDE genes (signs of suppression of cell division), sfmACDHF genes (signs of adhesin synthesis), yjbEFGH and gfcA genes (signs of cell envelope stabilization). Moreover, THz radiation induced E. coli csg operon genes of amyloid biosynthesis. Electron microscopy revealed that the irradiated bacteria underwent increased aggregation; 20% of them formed bundle-like structures consisting of two to four pili clumped together. This could be the result of changes in the adhesive properties of the pili. We also found aberrations in cell wall structure in the middle part of the bacterial cell; these aberrations impaired the cell at the initial stages of division and resulted in accumulation of long rod-like cells. Overall, THz radiation was shown to have adverse effects on bacterial populations resulting in cells with abnormal morphology.

Microorganisms of Two Thermal Pools on Kunashir Island, Russia.

The Kuril Archipelago is a part of the Circum-Pacific Belt (Ring of Fire). These islands have numerous thermal springs. There are very few studies on these microbial communities, and none of them have been conducted by modern molecular biological methods. Here we performed the first metagenomic study on two thermophilic microbial communities of Kunashir Island. Faust Lake is hot (48 °C) and highly acidic (pH 2.0). We constructed 28 metagenome-assembled genomes as well as 17 16S ribosomal RNA sequences. We found that bottom sediments of Faust Lake are dominated by a single species of red algae belonging to the Cyanidiaceae family. Archaeans in Faust Lake are more diverse than bacteria but less abundant. The Tretyakovsky Thermal Spring is also hot (52 °C) but only weakly acidic (pH 6.0). It has much higher microbial diversity (233 metagenome-assembled genomes; 93 16S ribosomal RNAs) and is dominated by bacteria, with only several archaeans and one fungus. Despite their geographic proximity, these two thermal springs were found to not share any species. A comparison of these two lakes with other thermal springs of the Circum-Pacific Belt revealed that only a few members of the communities are shared among different locations.

Whole-Exome Sequencing Analysis of Human Semen Quality in Russian Multiethnic Population.

The global trend toward the reduction of human spermatogenic function observed in many countries, including Russia, raised the problem of extensive screening and monitoring of male fertility and elucidation of its genetic and ethnic mechanisms. Recently, whole-exome sequencing (WES) was developed as a powerful tool for genetic analysis of complex traits. We present here the first Russian WES study for identification of new genes associated with semen quality. The experimental 3 × 2 design of the WES study was based on the analysis of 157 samples including three ethnic groups-Slavs (59), Buryats (n = 49), and Yakuts (n = 49), and two different semen quality groups-pathozoospermia (n = 95) and normospermia (n = 62). Additionally, our WES study group was negative for complete AZF microdeletions of the Y-chromosome. The normospermia group included men with normal sperm parameters in accordance with the WHO-recommended reference limit. The pathozoospermia group included men with impaired semen quality, namely, with any combined parameters of sperm concentration <15 × 106/ml, and/or progressive motility <32%, and/or normal morphology <4%. The WES was performed for all 157 samples. Subsequent calling and filtering of variants were carried out according to the GATK Best Practices recommendations. On the genotyping stage, the samples were combined into four cohorts: three sets corresponded to three ethnic groups, and the fourth set contained all the 157 whole-exome samples. Association of the obtained polymorphisms with semen quality parameters was investigated using the χ2 test. To prioritize the obtained variants associated with pathozoospermia, their effects were determined using Ensembl Variant Effect Predictor. Moreover, polymorphisms located in genes expressed in the testis were revealed based on the genomic annotation. As a result, the nine potential SNP markers rs6971091, rs557806, rs610308, rs556052, rs1289658, rs278981, rs1129172, rs12268007, and rs17228441 were selected for subsequent verification on our previously collected population sample (about 1,500 males). The selected variants located in seven genes FAM71F1, PPP1R15A, TRIM45, PRAME, RBM47, WDFY4, and FSIP2 that are expressed in the testis and play an important role in cell proliferation, meiosis, and apoptosis.

A Bioinformatics Model of Human Diseases on the Basis of Differentially Expressed Genes (of Domestic Versus Wild Animals) That Are Orthologs of Human Genes Associated with Reproductive-Potential Changes.

Earlier, after our bioinformatic analysis of single-nucleotide polymorphisms of TATA-binding protein-binding sites within gene promoters on the human Y chromosome, we suggested that human reproductive potential diminishes during self-domestication. Here, we implemented bioinformatics models of human diseases using animal in vivo genome-wide RNA-Seq data to compare the effect of co-directed changes in the expression of orthologous genes on human reproductive potential and during the divergence of domestic and wild animals from their nearest common ancestor (NCA). For example, serotonin receptor 3A (HTR3A) deficiency contributes to sudden death in pregnancy, consistently with Htr3a underexpression in guinea pigs (Cavia porcellus) during their divergence from their NCA with cavy (C. aperea). Overall, 25 and three differentially expressed genes (hereinafter, DEGs) in domestic animals versus 11 and 17 DEGs in wild animals show the direction consistent with human orthologous gene-markers of reduced and increased reproductive potential. This indicates a reliable association between DEGs in domestic animals and human orthologous genes reducing reproductive potential (Pearson's χ2 test p < 0.001, Fisher's exact test p < 0.05, binomial distribution p < 0.0001), whereas DEGs in wild animals uniformly match human orthologous genes decreasing and increasing human reproductive potential (p > 0.1; binomial distribution), thus enforcing the norm (wild type).

Species Delimitation of the Eisenia nordenskioldi Complex (Oligochaeta, Lumbricidae) Using Transcriptomic Data.

Eisenia nordenskioldi (Eisen, 1879) is the only autochthonous Siberian earthworm with a large distribution that ranges from tundra to steppe and broadleaved forests. This species has a very high morphological, ecological, karyological, and genetic diversity, so it was proposed that E. nordenskioldi should be split into several species. However, the phylogeny of the complex was unclear due to the low resolution of the methods used and the high diversity that should have been taken into account. We investigated this question by (1) studying the diversity of the COI gene of E. nordenskioldi throughout its range and (2) sequencing transcriptomes of different genetic lineages to infer its phylogeny. We found that E. nordenskioldi is monophyletic and is split into two clades. The first one includes the pigmented genetic lineages widespread in the northern and western parts of the distribution, and the second one originating from the southern and southeastern part of the species' range and representing both pigmented and non-pigmented forms. We propose to split the E. nordenskioldi complex into two species, E. nordenskioldi and Eisenia sp. 1 (aff. E. nordenskioldi), corresponding to these two clades. The currently recognized non-pigmented subspecies E. n. pallida will be abolished as a polyphyletic and thus a non-natural taxon, while Eisenia sp. 1 will be expanded to include several lineages earlier recognized as E. n. nordenskioldi and E. n. pallida.

Disruptive natural selection by male reproductive potential prevents underexpression of protein-coding genes on the human Y chromosome as a self-domestication syndrome.

BACKGROUND: In population ecology, the concept of reproductive potential denotes the most vital indicator of chances to produce and sustain a healthy descendant until his/her reproductive maturity under the best conditions. This concept links quality of life and longevity of an individual with disease susceptibilities encoded by his/her genome. Female reproductive potential has been investigated deeply, widely, and comprehensively in the past, but the male one has not received an equal amount of attention. Therefore, here we focused on the human Y chromosome and found candidate single-nucleotide polymorphism (SNP) markers of male reproductive potential. RESULTS: Examining in silico (i.e., using our earlier created Web-service SNP_TATA_Z-tester) all 1206 unannotated SNPs within 70 bp proximal promoters of all 63 Y-linked genes, we found 261 possible male-reproductive-potential SNP markers that can significantly alter the binding affinity of TATA-binding protein (TBP) for these promoters. Among them, there are candidate SNP markers of spermatogenesis disorders (e.g., rs1402972626), pediatric cancer (e.g., rs1483581212) as well as male anxiety damaging family relationships and mother's and children's health (e.g., rs187456378). First of all, we selectively verified in vitro both absolute and relative values of the analyzed TBP-promoter affinity, whose Pearson's coefficients of correlation between predicted and measured values were r = 0.84 (significance p <  0.025) and r = 0.98 (p <  0.025), respectively. Next, we found that there are twofold fewer candidate SNP markers decreasing TBP-promoter affinity relative to those increasing it, whereas in the genome-wide norm, SNP-induced damage to TBP-promoter complexes is fourfold more frequent than SNP-induced improvement (p <  0.05, binomial distribution). This means natural selection against underexpression of these genes. Meanwhile, the numbers of candidate SNP markers of an increase and decrease in male reproductive potential were indistinguishably equal to each other (p <  0.05) as if male self-domestication could have happened, with its experimentally known disruptive natural selection. Because there is still not enough scientific evidence that this could have happened, we discuss the human diseases associated with candidate SNP markers of male reproductive potential that may correspond to domestication-related disorders in pets. CONCLUSIONS: Overall, our findings seem to support a self-domestication syndrome with disruptive natural selection by male reproductive potential preventing Y-linked underexpression of a protein.

The first mitochondrial genome of the relic Acanthobdella peledina (Annelida, Acanthobdellida).

In this study, we assembled a complete mitochondrial genome of the Acanthobdella peledina sample from the Pitea River, Sweden. Thirty-six genes of the mitogenome sequence were identified, including 13 protein-coding genes (PCGs), 2 ribosomal genes (12S and 16S), 21 transport RNA genes, and 1 control region. The complete mitogenome is 14,640 bp long and A + T biased (70.59%).

Transcriptomic analysis confirms differences among nuclear genomes of cryptic earthworm lineages living in sympatry.

BACKGROUND: Many earthworm species demonstrate significant cryptic diversity, with several highly diverged mitochondrial lineages found within most of the taxa studied to date. The status of differences between these lineages on the nuclear level is still unclear. Because of widespread polyploidy in earthworms, most studies were limited to two nuclear loci, the ribosomal and the histone clusters. Here we attempted to elucidate the status of a set of genetic lineages within Eisenia nordenskioldi nordenskioldi, an earthworm species from Northern Asia with high intraspecific diversity. We performed RNA-seq on an IonTorrent platform for five specimens of this species belonging to five genetic lineages, as well as two outgroups from the family Lumbricidae, the congenetic E. andrei, and Lumbricus rubellus. RESULTS: We de novo assembled transcriptomes and constructed datasets of genes present in all seven specimens using broad (ProteinOrtho; 809 genes) and narrow (HaMStR; 203 genes) ortholog assignment. The majority of orthologs had identical amino acid sequences in all studied specimens, which we believe was due to strong bias towards the most conserved genes. However, for the rest of genes the differences among the lineages were lower than those between them and the congeneric E. andrei. Both datasets yielded phylogenetic trees with the same topology. E. n. nordenskioldi was found to be monophyletic. The differences on the genetic level had no concordance with geography, implying complex history of dispersal. CONCLUSIONS: We found that genetic lineages of E. n. nordenskioldi are genetically distinct on nuclear level and probably diverged long ago. Current data implies that they might even represent distinct species within the E. nordenskioldi species complex.

Metabolic pathways and genes identified by RNA-seq analysis of barley near-isogenic lines differing by allelic state of the Black lemma and pericarp (Blp) gene.

BACKGROUND: Some plant species have 'melanin-like' black seed pigmentation. However, the chemical and genetic nature of this 'melanin-like' black pigment have not yet been fully explored due to its complex structure and ability to withstand almost all solvents. Nevertheless, identification of genetic networks participating in trait formation is key to understanding metabolic processes involved in the expression of 'melanin-like' black seed pigmentation. The aim of the current study was to identify differentially expressed genes (DEGs) in barley near-isogenic lines (NILs) differing by allelic state of the Blp (black lemma and pericarp) locus. RESULTS: RNA-seq analysis of six libraries (three replicates for each line) was performed. A total of 957 genome fragments had statistically significant changes in expression levels between lines BLP and BW, with 632 fragments having increased expression levels in line BLP and 325 genome fragments having decreased expression. Among identified DEGs, 191 genes were recognized as participating in known pathways. Among these were metabolic pathways including 'suberin monomer biosynthesis', 'diterpene phytoalexins precursors biosynthesis', 'cutin biosynthesis', 'cuticular wax biosynthesis', and 'phenylpropanoid biosynthesis, initial reactions'. Differential expression was confirmed by real-time PCR analysis of selected genes. CONCLUSIONS: Metabolic pathways and genes presumably associated with black lemma and pericarp colour as well as Blp-associated resistance to oxidative stress and pathogens, were revealed. We suggest that the black pigmentation of lemmas and pericarps is related to increased level of phenolic compounds and their oxidation. The effect of functional Blp on the synthesis of ferulic acid and other phenolic compounds can explain the increased antioxidant capacity and biotic and abiotic stress tolerance of black-grained cereals. Their drought tolerance and resistance to diseases affecting the spike may also be related to cuticular wax biosynthesis. In addition, upregulated synthesis of phytoalexins, suberin and universal stress protein (USP) in lemmas and pericarps of the Blp carriers may contribute to their increased disease resistance. Further description of the DEGs haplotypes and study of their association with physiological characteristics may be useful for future application in barley pre-breeding.

Computer Analysis of Glioma Transcriptome Profiling: Alternative Splicing Events.

Here we present the analysis of alternative splicing events on an example of glioblastoma cell culture samples using a set of computer tools in combination with database integration. The gene expression profiles of glioblastoma were obtained from cell culture samples of primary glioblastoma which were isolated and processed for RNA extraction. Transcriptome profiling of normal brain samples and glioblastoma were done by Illumina sequencing. The significant differentially expressed exon-level probes and their corresponding genes were identified using a combination of the splicing index method. Previous studies indicated that tumor-specific alternative splicing is important in the regulation of gene expression and corresponding protein functions during cancer development. Multiple alternative splicing transcripts have been identified as progression markers, including generalized splicing abnormalities and tumor- and stage-specific events. We used a set of computer tools which were recently applied to analysis of gene expression in laboratory animals to study differential splicing events. We found 69 transcripts that are differentially alternatively spliced. Three cancer-associated genes were considered in detail, in particular: APP (amyloid beta precursor protein), CASC4 (cancer susceptibility candidate 4) and TP53. Such alternative splicing opens new perspectives for cancer research.