The draft genome of the pest tephritid fruit fly Bactrocera tryoni: resources for the genomic analysis of hybridising species

Background

The tephritid fruit flies include a number of economically important pests of horticulture, with a large accumulated body of research on their biology and control. Amongst the Tephritidae, the genus Bactrocera, containing over 400 species, presents various species groups of potential utility for genetic studies of speciation, behaviour or pest control. In Australia, there exists a triad of closely-related, sympatric Bactrocera species which do not mate in the wild but which, despite distinct morphologies and behaviours, can be force-mated in the laboratory to produce fertile hybrid offspring. To exploit the opportunities offered by genomics, such as the efficient identification of genetic loci central to pest behaviour and to the earliest stages of speciation, investigators require genomic resources for future investigations.

Results

We produced a draft de novo genome assembly of Australia’s major tephritid pest species, Bactrocera tryoni. The male genome (650 -700 Mbp) includes approximately 150Mb of interspersed repetitive DNA sequences and 60Mb of satellite DNA. Assessment using conserved core eukaryotic sequences indicated 98% completeness. Over 16,000 MAKER-derived gene models showed a large degree of overlap with other Dipteran reference genomes. The sequence of the ribosomal RNA transcribed unit was also determined. Unscaffolded assemblies of B. neohumeralis and B. jarvisi were then produced; comparison with B. tryoni showed that the species are more closely related than any Drosophila species pair. The similarity of the genomes was exploited to identify 4924 potentially diagnostic indels between the species, all of which occur in non-coding regions.

Conclusions

This first draft B. tryoni genome resembles other dipteran genomes in terms of size and putative coding sequences. For all three species included in this study, we have identified a comprehensive set of non-redundant repetitive sequences, including the ribosomal RNA unit, and have quantified the major satellite DNA families. These genetic resources will facilitate the further investigations of genetic mechanisms responsible for the behavioural and morphological differences between these three species and other tephritids. We have also shown how whole genome sequence data can be used to generate simple diagnostic tests between very closely-related species where only one of the species is scaffolded.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1153) contains supplementary material, which is available to authorized users.     

Sulfate reduction and S-oxidation in a moorland pool sediment

In an oligotrophic moorland pool in The Netherlands, S cycling near the sediment/water boundary was investigated by measuring (1) SO₄ ^2– reduction rates in the sediment, (2) depletion of SO₄ ^2– in the overlying water column and (3) release of³⁵S from the sediment into the water column. Two locations differing in sediment type (highly organic and sandy) were compared, with respect to reduction rates and depletion of SO₄ ^2– in the overlying water.Sulfate reduction rates in sediments of an oligotrophic moorland pool were estimated by diagenetic modelling and whole core³⁵SO₄ ^2– injection. Rates of SO₄ ^2– consumption in the overlying water were estimated by changes in SO₄ ^2– concentration over time in in situ enclosures. Reduction rates ranged from 0.27–11.2 mmol m^–2 d^–1. Rates of SO₄ ^2– uptake from the enclosed water column varied from –0.5, –0.3 mmol m^–2 d^–1 (November) to 0.43–1.81 mmol m^–2 d^–1 (July, August and April). Maximum rates of oxidation to SO₄ ^2– in July 1990 estimated by combination of SO₄ ^2– reduction rates and rates of in situ SO₄ ^2– uptake in the enclosed water column were 10.3 and 10.5 mmol m^–2 d^–1 at an organic rich and at a sandy site respectively.Experiments with³⁵S^2– and³⁵SO₄ ^2– tracer suggested (1) a rapid formation of organically bound S from dissimilatory reduced SO₄ ^2– and (2) the presence of mainly non SO₄ ^2–-S derived from reduced S transported from the sediment into the overlying water. A³⁵S^2– tracer experiment showed that about 7% of³⁵S^2– injected at 1 cm depth in a sediment core was recovered in the overlying water column.Sulfate reduction rates in sediments with higher volumetric mass fraction of organic matter did not significantly differ from those in sediments with a lower mass fraction of organic matter.Corresponding author     

Studies on the incorporation of precursors into purine and pyrimidine nucleotides via 'de novo' and 'salvage' pathways in normal lymphocytes and lymphoblastic cell-line cells

The incorporation of radioactive precursors into purine and pyrimidine nucleotides via 'de novo' and 'salvage' pathways was measured in normal lymphocytes, resting as well as proliferating, and lymphoblastic cell-line cells (MOLT-3). Lymphocytes stimulated with anti-CD3 were taken as actively proliferating lymphocytes (35% in the S-phase, 40 h after stimulation). The incorporation of the precursors in the purine and pyrimidine ribonucleotides was measured by a combination of anion-exchange high-performance liquid chromatography (HPLC) and on-line radioactivity measurement. The actively proliferating normal lymphocytes and MOLT-3 cells incorporated 30-500 times more of the various precursors in the ribonucleotides compared to normal resting lymphocytes. The imbalance in the nucleotide pool found in proliferating normal and lymphoblastic cells was reflected in the incorporation pattern of the various precursors. The activities of the branch-point enzymes IMP dehydrogenase and CTP synthetase most likely determine the differences in the composition of the nucleotide pools between resting and proliferating cells.     

Molecular evolution of olfactomedin

Olfactomedin is a secreted polymeric glycoprotein of unknown function,originally discovered at the mucociliary surface of the amphibian olfactoryneuroepithelium and subsequently found throughout the mammalian brain. As afirst step toward elucidating the function of olfactomedin, itsphylogenetic history was examined to identify conserved structural motifs.Such conserved motifs may have functional significance and provide targetsfor future mutagenesis studies aimed at establishing the function of thisprotein. Previous studies revealed 33% amino acid sequence identity betweenrat and frog olfactomedins in their carboxyl terminal segments. Furtheranalysis, however, reveals more extensive homologies throughout themolecule. Despite significant sequence divergence, cysteines essential forhomopolymer formation such as the CXC motif near the amino terminus areconserved, as is the characteristic glycosylation pattern, suggesting thatthese posttranslational modifications are essential for function.Furthermore, evolutionary analysis of a region of 53 amino acids of fish,frog, rat, mouse, and human olfactomedins indicates that an ancestralolfactomedin gene arose before the evolution of terrestrial vertebrates andevolved independently in teleost, amphibian, and mammalian lineages.Indeed, a distant olfactomedin homolog was identified in Caenorhabditiselegans. Although the amino acid sequence of this invertebrate protein islonger and highly divergent compared with its vertebrate homologs, theprotein from C. elegans shows remarkable similarities in terms of conservedmotifs and posttranslational modification sites. Six universally conservedmotifs were identified, and five of these are clustered in the carboxylterminal half of the protein. Sequence comparisons indicate that evolutionof the N-terminal half of the molecule involved extensive insertions anddeletions; the C-terminal segment evolved mostly through point mutations,at least during vertebrate evolution. The widespread occurrence ofolfactomedin among vertebrates and invertebrates underscores the notionthat this protein has a function of universal importance. Furthermore,extensive modification of its N-terminal half and the acquisition of aC-terminal SDEL endoplasmic-reticulum- targeting sequence may have enabledolfactomedin to adopt new functions in the mammalian central nervoussystem.     

Molecular phylogeny and divergence times of drosophilid species

The phylogenetic relationships and divergence times of 39 drosophilid species were studied by using the coding region of the Adh gene. Four genera--Scaptodrosophila, Zaprionus, Drosophila, and Scaptomyza (from Hawaii)--and three Drosophila subgenera--Drosophila, Engiscaptomyza, and Sophophora--were included. After conducting statistical analyses of the nucleotide sequences of the Adh, Adhr (Adh-related gene), and nuclear rRNA genes and a 905-bp segment of mitochondrial DNA, we used Scaptodrosophila as the outgroup. The phylogenetic tree obtained showed that the first major division of drosophilid species occurs between subgenus Sophophora (genus Drosophila) and the group including subgenera Drosophila and Engiscaptomyza plus the genera Zaprionus and Scaptomyza. Subgenus Sophophora is then divided into D. willistoni and the clade of D. obscura and D. melanogaster species groups. In the other major drosophilid group, Zaprionus first separates from the other species, and then D. immigrans leaves the remaining group of species. This remaining group then splits into the D. repleta group and the Hawaiian drosophilid cluster (Hawaiian Drosophila, Engiscaptomyza, and Scaptomyza). Engiscaptomyza and Scaptomyza are tightly clustered. Each of the D. repleta, D. obscura, and D. melanogaster groups is monophyletic. The splitting of subgenera Drosophila and Sophophora apparently occurred about 40 Mya, whereas the D. repleta group and the Hawaiian drosophilid cluster separated about 32 Mya. By contrast, the splitting of Engiscaptomyza and Scaptomyza occurred only about 11 Mya, suggesting that Scaptomyza experienced a rapid morphological evolution. The D. obscura and D. melanogaster groups apparently diverged about 25 Mya. Many of the D. repleta group species studied here have two functional Adh genes (Adh-1 and Adh-2), and these duplicated genes can be explained by two duplication events.      

Validated predictive modelling of the environmental resistome

Multi-drug-resistant bacteria pose a significant threat to public health. The role of the environment in the overall rise in antibiotic-resistant infections and risk to humans is largely unknown. This study aimed to evaluate drivers of antibiotic-resistance levels across the River Thames catchment, model key biotic, spatial and chemical variables and produce predictive models for future risk assessment. Sediment samples from 13 sites across the River Thames basin were taken at four time points across 2011 and 2012. Samples were analysed for class 1 integron prevalence and enumeration of third-generation cephalosporin-resistant bacteria. Class 1 integron prevalence was validated as a molecular marker of antibiotic resistance; levels of resistance showed significant geospatial and temporal variation. The main explanatory variables of resistance levels at each sample site were the number, proximity, size and type of surrounding wastewater-treatment plants. Model 1 revealed treatment plants accounted for 49.5% of the variance in resistance levels. Other contributing factors were extent of different surrounding land cover types (for example, Neutral Grassland), temporal patterns and prior rainfall; when modelling all variables the resulting model (Model 2) could explain 82.9% of variations in resistance levels in the whole catchment. Chemical analyses correlated with key indicators of treatment plant effluent and a model (Model 3) was generated based on water quality parameters (contaminant and macro- and micro-nutrient levels). Model 2 was beta tested on independent sites and explained over 78% of the variation in integron prevalence showing a significant predictive ability. We believe all models in this study are highly useful tools for informing and prioritising mitigation strategies to reduce the environmental resistome.     

心脏特异表达LMNA^E82K转基因小鼠的建立

目的建立心脏特异表达LMNAE82K转基因小鼠,为研究LMNAE82K与心肌病发病机制的关系提供工具动物。方法把LMNAE82K基因插入α-MHC启动子下游,构建转基因表达载体,显微注射法建立C57BL/6JLMNAE82K转基因小鼠,PCR鉴定转基因小鼠的基因型,采用Western Blot鉴定LMNAE82K在心脏组织中的表达,H＆E染色和超声检测转基因小鼠心脏的病理改变。结果建立了2个心脏组织特异表达LMNAE82K的转基因小鼠品系。超声检查显示转基因小鼠心室壁变薄,收缩期容积和舒张期容积增加,射血分数及短轴缩短率降低。结论LMNAE82K转基因小鼠具有LMNAE82K引起的家族性扩心病有类似的病理变化,为研究LMNAE82K与心肌病发病机制的关系的研究提供了有价值的疾病动物模型。