Transposition burst of mariner-like elements in the sequenced genome of Rhodnius prolixus

Transposable elements (TEs) are widespread in insect's genomes. However, there are wide differences in the proportion of the total DNA content occupied by these repetitive sequences in different species. We have analyzed the TEs present in R. prolixus (vector of the Chagas disease) and showed that 3.0% of this genome is occupied by Class II TEs, belonging mainly to the Tc1-mariner superfamily (1.65%) and MITEs (1.84%). Interestingly, most of this genomic content is due to the expansion of two subfamilies belonging to: irritans himar, a well characterized subfamily of mariners, and prolixus1, one of the two novel subfamilies here described. The high amount of sequences in these subfamilies suggests that bursts of transposition occurred during the life cycle of this family. In an attempt to characterize these elements, we performed an in silico analysis of the sequences corresponding to the DDD/E domain of the transposase gene. We performed an evolutionary analysis including network and Bayesian coalescent-based methods in order to infer the dynamics of the amplification, as well as to estimate the time of the bursts identified in these subfamilies. Given our data, we hypothesized that the TE expansions occurred around the time of speciation of R. prolixus around 1.4 mya. This suggestion lays on the “Transposon Model” of TE evolution, in which the members of a TE population that are replicative active are present at multiple loci in the genome, but their replicative potential varies, and of the “Life Cycle Model” that states that when present-day TEs have been involved in amplification bursts, they share an ancestral copy that dates back to this initial amplification.     

Genetic diversity of wild population of Pyropia haitanensis based on SSR analysis

Genetic diversity analysis was conducted on 80 individuals of 4 populations of non-bred Pyropia haitanensis by simple sequence repeat (SSR) method. Using 15 pairs of microsatellite primers, 37 polymorphic loci were amplified, representing 94.9% of all loci. At the population level, the percentage of polymorphic bands (P) and polymorphism information content (PIC) were 66.67–84.62% and 0.481–0.488, with average value at 73.72% and 0.483, respectively. Expected heterozygosity (H_e) and Shannon's information index (I) were 0.279 and 0.434, respectively, at the species level, and 0.233 and 0.356 at population level. According to the coefficient of gene differentiation (G_ST), a large proportion of genetic variance (83.6%) of P. haitanensis was among individuals within populations, only 16.4% genetic variance was among populations, which was identified with the moderate gene flow value (N_m = 2.542). UPGMA clustered the 4 populations into 3 groups, and no significant correlation was found between the genetic distance and the corresponding geographic distance among the populations.     

The genetic diversity and population structure of wild soybean evaluated by chloroplast and nuclear gene sequences

Glycine soja, also called wild soybean, is the wild ancestor of domesticated soybean (Glycine max), and one of the world's major cultivated crops. Wild soybean is a valuable resource for the breeding of cultivated soybean and harbors useful genes or agronomic traits. To use and conserve this valuable resource, we conducted a study to evaluate the genetic diversity and population structure of wild soybean using the sequencing data of two nuclear loci (AF105221 and PhyB) and one chloroplast locus (trnQ-rps16) of more than 600 individuals representing 53 populations throughout the natural distribution range. The results showed that most of the variation was found within the populations and groups, but significant genetic differentiation was also detected among different eco-geographical groups. Correlations between genetic and geographical distance at all the loci were consistent with the isolation by distance gene flow model. G. soja exhibited the highest genetic diversity in middle and downstream of Yangzi River (MDYR) region, followed by North East China (NEC), and was the lowest in North West China (NWC). We concluded that both in situ and ex situ conservation strategies required for wild soybean populations, especially which are native to MDYR and NEC regions.     

An independent occurrence of the chimeric P450 enzyme CYP337B3 of Helicoverpa armigera confers cypermethrin resistance in Pakistan

The increasing resistance level of insect pest species is a major concern to agriculture worldwide. The cotton bollworm, Helicoverpa armigera, is one of the most important pest species due to being highly polyphagous, geographically widespread, and resistant towards many chemical classes of insecticides. We previously described the mechanism of fenvalerate resistance in Australian populations conferred by the chimeric cytochrome P450 monooxygenase CYP337B3, which arose by unequal crossing-over between CYP337B1 and CYP337B2. Here, we show that this mechanism is also present in the cypermethrin-resistant FSD strain from Pakistan. The Pakistani and the Australian CYP337B3 alleles differ by 18 synonymous and three nonsynonymous SNPs and additionally in the length and sequence of the intron. Nevertheless, the activity of both CYP337B3 proteins is comparable. We demonstrate that CYP337B3 is capable of metabolizing cypermethrin (trans- and especially cis-isomers) to the main metabolite 4'-hydroxycypermethrin, which exhibits no intrinsic toxicity towards susceptible larvae. In a bioassay, CYP337B3 confers a 7-fold resistance towards cypermethrin in FSD larvae compared to susceptible larvae from the Australian TWB strain lacking CYP337B3. Linkage analysis shows that presence of CYP337B3 accounts for most of the cypermethrin resistance in the FSD strain; up-regulation of other P450s in FSD plays no detectable role in resistance. The presence or absence of CYP337B3 can be easily detected by a simple PCR screen, providing a powerful tool to rapidly distinguish resistant from susceptible individuals in the field and to determine the geographical distribution of this resistance gene. Our results suggest that CYP337B3 evolved twice independently by unequal crossing-over between CYP337B2 and two different CYP337B1 alleles.     

Marked differences of haplotype tagging SNP distribution,linkage, and haplotype profile of IL23 receptor gene in Roma and Hungarian population samples

Polymorphisms of the interleukin-23 receptor (IL23R) gene have been found to play an important role in the development of several autoimmune diseases. We examined five susceptible (rs10889677, rs1004819, rs2201841, rs11805303, rs11209032), one protective (rs7517847) and two neutral variants (rs7530511, rs1884444) of the IL23R gene in pooled DNA of healthy Roma (Gipsy) and Hungarian population samples. Our aim was to determine the genetic variability of the major haplotype tagging polymorphisms, and the haplotype profile of IL23R between the two groups. We analyzed 273 healthy Roma and 253 Hungarian DNA samples using PCR/RFLP assay. Comparing the five susceptible conferring alleles, there were significant increase (p < 0.05), while in the protective alleles, there were decrease in the allele frequencies in Roma population (p < 0.05). One of the neutral alleles showed increase, the another one did not differ between the two groups. The haplotype analysis of the SNPs revealed fundamentally different association types of SNPs in the two groups; moreover, the frequencies of the various haplotypes also exhibited strong differences, as of ht4 and ht5 haplotypes were significantly higher, whereas the frequencies of ht2 and ht3 haplotypes were significantly lower in the Roma population than in Hungarians (p < 0.05). The data presented here show profound differences in the IL23R genetic profiles in the Roma population, that likely has also clinical implications in respect their possible role in the development of certain immunological diseases.     

Localization of the ultraviolet-sensor Opn5m and its effect on myopia-related gene expression in the late-embryonic chick eye

Recent studies show that exposure to ultraviolet (UV) light suppresses ocular elongation, which causes myopia development. However, the specific mechanisms of this process have not been elucidated. A UV-sensor, Opsin 5 (Opn5) mRNA was shown to be present in extraretinal tissues. To test the possibility that UV-signals mediated by Opn5 would have a direct effect on the outer connective tissues of the eye, we first examined the expression patterns of a mammalian type Opn5 (Opn5m) in the late-embryonic chicken eye. Quantitative PCR showed Opn5m mRNA expression in the cornea and sclera. The anti-Opn5m antibody stained a small subset of cells in the corneal stroma and fibrous sclera. We next assessed the effect of UV-A (375 nm) irradiation on the chicken fibroblast cell line DF-1 overexpressing chicken Opn5m. UV-A irradiation for 30 min significantly increased the expression of Early growth response 1 (Egr1), known as an immediate early responsive gene, and of Matrix metalloproteinase 2 (Mmp2) in the presence of retinal chromophore 11-cis-retinal. In contrast, expression of Transforming growth factor beta 2 and Tissue inhibitor of metalloproteinase 2 was not significantly altered. These results indicate that UV-A absorption by Opn5m can upregulate the expression levels of Egr1 and Mmp2 in non-neuronal, fibroblasts. Taken together with the presence of Opn5m in the cornea and sclera, it is suggested that UV-A signaling mediated by Opn5 in the extraretinal ocular tissues could influence directly the outer connective tissues of the chicken late-embryonic eye.     

Seasonal patterns of microcystin-producing and non-producing Planktothrix rubescens genotypes in a deep pre-alpine lake

The filamentous cyanobacterium Planktothrix rubescens produces secondary metabolites called microcystins (MC) that are potent toxins for most eukaryotes, including zooplankton grazers, cattle and humans. P. rubescens occurs in many deep and thermally stratified lakes throughout Europe. In Lake Zurich (Switzerland), it re-appeared in the 1970s concomitant with decreasing eutrophication. Since then, P. rubescens has become the dominant species in this major drinking water reservoir, where it forms massive metalimnetic blooms during late summer. These cyanobacteria harbor subpopulations of non-MC producers, but little is known about the environmental factors affecting the success of such genotypes. The non-MC-producing subpopulation of P. rubescens was studied using a quantitative real-time PCR (qPCR) assay on the MC synthetase (mcy) gene cluster that targets a deletion on the mcyH and mcyA genes, which inactivates MC biosynthesis. Two complementary qPCR assays were used to assess the total population abundance (based on the 16S rDNA gene) and the mcy gene copy number (based on a conserved region in the adenylation domain of the mcyB gene). The objective was to evaluate the seasonal patterns of the share of non-MC-producing filaments in the total P. rubescens population. The mcyHA mutants were present in low proportions (up to 14%) throughout the year. Their highest relative abundances occurred during the winter mixis, when total concentrations of P. rubescens were minimal. The MC deficient mutants seemed to better survive in sparse populations, possibly because of lower grazing pressure and a consequently reduced need for MC-mediated protection. Alternatively, the mutants might cope better with the sub-optimal, stressful pressure and light conditions during the winter mixis. Altogether, our results suggest that subtle trade-offs might seasonally determine the proportions of non-MC producers within P. rubescens populations.     

LIM-homeodomain transcription factor Awh is a key component activating all three fibroin genes,fibH, fibL and fhx,in the silk gland of the silkworm,Bombyx mori

In the silkworm Bombyx mori, three fibroin genes, fibroin-heavy-chain (fibH), fibroin-light-chain (fibL) and fibrohexamerin (fhx), are coexpressed only in the posterior silk gland (PSG) cells, while the sericin genes encoding silk glue proteins are expressed in the middle silk gland (MSG) cells. Silk gland factor-2 (SGF-2) is a PSG-specific activator complex of fibH, composed of a LIM-homeodomain protein, Awh, and its cofactors, Ldb and Lcaf. We investigated whether SGF-2 can activate other fibroin genes using transgenic silkworms. The genes for Ldb and Lcaf were expressed ubiquitously in various tissues, while the gene for Awh was expressed strictly specific in PSG of the wild type silkworms. Misexpression of Awh in transgenic silkworms induced ectopic expression of fibL and fhx as well as fibH in MSG. Coincidently with the induction of fibL and fhx by Awh, binding of SGF-2 to the promoter of fibL and fhx was detected in vitro, and SGF-2 binds directly to the fhx core promoter. Ectopic expression of the fibroin genes was observed at high levels in the middle part of MSG. Moreover, fibL and fhx were induced in the anterior silk gland (ASG) of the transgenic silkworms, but fibH was not. These results indicate that Awh is a key activator of all three fibroin genes, and the activity is probably regulated in conjunction with additional factors.     

Integrated modeling of protein-coding genes in the Manduca sexta genome using RNA-Seq data from the biochemical model insect

The genome sequence of Manduca sexta was recently determined using 454 technology. Cufflinks and MAKER2 were used to establish gene models in the genome assembly based on the RNA-Seq data and other species' sequences. Aided by the extensive RNA-Seq data from 50 tissue samples at various life stages, annotators over the world (including the present authors) have manually confirmed and improved a small percentage of the models after spending months of effort. While such collaborative efforts are highly commendable, many of the predicted genes still have problems which may hamper future research on this insect species. As a biochemical model representing lepidopteran pests, M. sexta has been used extensively to study insect physiological processes for over five decades. In this work, we assembled Manduca datasets Cufflinks 3.0, Trinity 4.0, and Oases 4.0 to assist the manual annotation efforts and development of Official Gene Set (OGS) 2.0. To further improve annotation quality, we developed methods to evaluate gene models in the MAKER2, Cufflinks, Oases and Trinity assemblies and selected the best ones to constitute MCOT 1.0 after thorough crosschecking. MCOT 1.0 has 18,089 genes encoding 31,666 proteins: 32.8% match OGS 2.0 models perfectly or near perfectly, 11,747 differ considerably, and 29.5% are absent in OGS 2.0. Future automation of this process is anticipated to greatly reduce human efforts in generating comprehensive, reliable models of structural genes in other genome projects where extensive RNA-Seq data are available.     

LAMMER kinase contributes to genome stability in Ustilago maydis

Here we report identification of the lkh1 gene encoding a LAMMER kinase homolog (Lkh1) from a screen for DNA repair-deficient mutants in Ustilago maydis. The mutant allele isolated results from a mutation at glutamine codon 488 to a stop codon that would be predicted to lead to truncation of the carboxy-terminal kinase domain of the protein. This mutant (lkh1^Q488*) is highly sensitive to ultraviolet light, methyl methanesulfonate, and hydroxyurea. In contrast, a null mutant (lkh1Δ) deleted of the entire lkh1 gene has a less severe phenotype. No epistasis was observed when an lkh1^Q488* rad51Δ double mutant was tested for genotoxin sensitivity. However, overexpressing the gene for Rad51, its regulator Brh2, or the Brh2 regulator Dss1 partially restored genotoxin resistance of the lkh1Δ and lkh1^Q488* mutants. Deletion of lkh1 in a chk1Δ mutant enabled these double mutant cells to continue to cycle when challenged with hydroxyurea. lkh1Δ and lkh1^Q488* mutants were able to complete the meiotic process but exhibited reduced heteroallelic recombination and aberrant chromosome segregation. The observations suggest that Lkh1 serves in some aspect of cell cycle regulation after DNA damage or replication stress and that it also contributes to proper chromosome segregation in meiosis.     

Complete mitochondrial genome of Pycnonotus xanthorrhous (Passeriformes,Pycnonotidae) and phylogenetic consideration

The complete mitochondrial genome (mitogenome) of Pycnonotus xanthorrhous was sequenced via next generation sequencing. The full length of the circular genome is 16,952 bp. It consists of 37 typical animal mitochondrial genes including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) and 2 ribosomal RNA (rRNA) genes. P. xanthorrhous also contains one control region (CR) and one pseudo-control region, and shares the identical gene arrangements with sequenced Pycnonotus spp. which differs from the typical vertebrates gene order. Phylogenetic analyses indicates that Passerida sensu stricto contains three major clades and the core Sylvioidea form a monophyletic group. Furthermore, we investigated the evolution of control region within this lineage and revealed the multiple independent origins of duplicate control region.     

Structural features,evolutionary relationships,and transcriptional regulation of C-type lectin-domain proteins in Manduca sexta

C-type lectins (CTLs) are a large family of Ca²⁺-dependent carbohydrate-binding proteins recognizing various glycoconjugates and functioning primarily in immunity and cell adhesion. We have identified 34 CTLDP (for CTL-domain protein) genes in the Manduca sexta genome, which encode proteins with one to three CTL domains. CTL-S1 through S9 (S for simple) have one or three CTL domains; immulectin-1 through 19 have two CTL domains; CTL-X1 through X6 (X for complex) have one or two CTL domains along with other structural modules. Nine simple CTLs and seventeen immulectins have a signal peptide and are likely extracellular. Five complex CTLs have both an N-terminal signal peptide and a C-terminal transmembrane region, indicating that they are membrane anchored. Immulectins exist broadly in Lepidoptera and lineage-specific gene duplications have generated three clusters of fourteen genes in the M. sexta genome, thirteen of which have similar expression patterns. In contrast to the family expansion, CTL-S1∼S6, S8, and X1∼X6 have 1:1 orthologs in at least four lepidopteran/dipteran/coleopteran species, suggestive of conserved functions in a wide range of holometabolous insects. Structural modeling suggests the key residues for Ca²⁺-dependent or independent binding of certain carbohydrates by CTL domains. Promoter analysis identified putative κB motifs in eighteen of the CTL genes, which did not have a strong correlation with immune inducibility in the mRNA or protein levels. Together, the gene identification, sequence comparisons, structure modeling, phylogenetic analysis, and expression profiling establish a solid foundation for future studies of M. sexta CTL-domain proteins.     

The steroid hormone 20-hydroxyecdysone promotes switching from autophagy to apoptosis by increasing intracellular calcium levels

Autophagy regulates cell survival (or cell death in several cases), whereas apoptosis regulates cell death. However, the relationship between autophagy and apoptosis and the regulative mechanism is unclear. We report that steroid hormone 20-hydroxyecdysone (20E) promotes switching from autophagy to apoptosis by increasing intracellular calcium levels in the midgut of the lepidopteran insect Helicoverpa armigera. Autophagy and apoptosis sequentially occurred during midgut programmed cell death under 20E regulation, in which lower concentrations of 20E induced microtubule-associated protein 1 light chain 3–phosphatidylethanolamine (LC3–II, also known as autophagy-related gene 8, ATG8) expression and autophagy. High concentrations of 20E induced cleavage of ATG5 to NtATG5 and pro-caspase-3 to active caspase-3, which led to a switch from autophagy to apoptosis. Blocking autophagy by knockdown of ATG5, ATG7, or ATG12, or with the autophagy inhibitor 3-methyladenine, inhibited 20E-induced autophagy and apoptosis. Blocking apoptosis by using the apoptosis inhibitor Ac-DEVD-CHO did not prevent 20E-induced autophagy, suggesting that apoptosis relies on autophagy. ATG5 knockdown resulted in abnormal pupation and delayed pupation time. High concentrations of 20E induced high levels of intracellular Ca²⁺, NtATG5, and active caspase-3, which mediated the switch from autophagy to apoptosis. Blocking 20E-mediated increase of cellular Ca²⁺ caused a decrease of NtATG5 and active caspase-3 and repressed the transformation from autophagy to apoptosis, thereby promoting cell survival. 20E induces an increase in the concentration of intracellular Ca²⁺, thereby switching autophagic cell survival to apoptotic cell death.