Genomic innovations, transcriptional plasticity and gene loss underlying the evolution and divergence of two highly polyphagous and invasive Helicoverpa pest species.

BACKGROUND: Helicoverpa armigera and Helicoverpa zea are major caterpillar pests of Old and New World agriculture, respectively. Both, particularly H. armigera, are extremely polyphagous, and H. armigera has developed resistance to many insecticides. Here we use comparative genomics, transcriptomics and resequencing to elucidate the genetic basis for their properties as pests. RESULTS: We find that, prior to their divergence about 1.5 Mya, the H. armigera/H. zea lineage had accumulated up to more than 100 more members of specific detoxification and digestion gene families and more than 100 extra gustatory receptor genes, compared to other lepidopterans with narrower host ranges. The two genomes remain very similar in gene content and order, but H. armigera is more polymorphic overall, and H. zea has lost several detoxification genes, as well as about 50 gustatory receptor genes. It also lacks certain genes and alleles conferring insecticide resistance found in H. armigera. Non-synonymous sites in the expanded gene families above are rapidly diverging, both between paralogues and between orthologues in the two species. Whole genome transcriptomic analyses of H. armigera larvae show widely divergent responses to different host plants, including responses among many of the duplicated detoxification and digestion genes. CONCLUSIONS: The extreme polyphagy of the two heliothines is associated with extensive amplification and neofunctionalisation of genes involved in host finding and use, coupled with versatile transcriptional responses on different hosts. H. armigera's invasion of the Americas in recent years means that hybridisation could generate populations that are both locally adapted and insecticide resistant.

Markov invariants, plethysms, and phylogenetics.

We explore model-based techniques of phylogenetic tree inference exercising Markov invariants. Markov invariants are group invariant polynomials and are distinct from what is known in the literature as phylogenetic invariants, although we establish a commonality in some special cases. We show that the simplest Markov invariant forms the foundation of the Log-Det distance measure. We take as our primary tool group representation theory, and show that it provides a general framework for analyzing Markov processes on trees. From this algebraic perspective, the inherent symmetries of these processes become apparent, and focusing on plethysms, we are able to define Markov invariants and give existence proofs. We give an explicit technique for constructing the invariants, valid for any number of character states and taxa. For phylogenetic trees with three and four leaves, we demonstrate that the corresponding Markov invariants can be fruitfully exploited in applied phylogenetic studies.

Population Genetic Structure and Potential Incursion Pathways of the Bluetongue Virus Vector Culicoides brevitarsis (Diptera: Ceratopogonidae) in Australia.

Culicoides brevitarsis is a vector of the bluetongue virus (BTV), which infects sheep and cattle. It is an invasive species in Australia with an assumed Asian/South East Asian origin. Using one mitochondrial marker (i.e., part of the cytochrome oxidase subunit I gene) and six nuclear markers, we inferred population genetic structure and possible incursion pathways for Australian C. brevitarsis. Nine mitochondrial haplotypes, with low nucleotide sequence diversity (0.0-0.7%) among these, were identified in a sample of 70 individuals from seven sites. Both sets of markers revealed a homogeneous population structure, albeit with evidence of isolation by distance and two genetically distinct clusters distributed along a north-to-south cline. No evidence of a cryptic species complex was found. The geographical distribution of the mitochondrial haplotypes is consistent with at least two incursion pathways into Australia since the arrival of suitable livestock hosts. By contrast, 15 mitochondrial haplotypes, with up to four times greater nucleotide sequence diversity (0.0-2.9%) among these, were identified in a sample of 16 individuals of the endemic C. marksi (sampled from a site in South Australia and another in New South Wales). A phylogenetic tree inferred using the mitochondrial marker revealed that the Australian and Japanese samples of C. brevitarsis are as evolutionarily different from one another as some of the other Australian species (e.g., C. marksi, C. henryi, C. pallidothorax) are. The phylogenetic tree placed four of the species endemic to Australia (C. pallidothorax, C. bundyensis, C. marksi, C. henryi) in a clade, with a fifth such species (C. bunrooensis) sharing a common ancestor with that clade and a clade comprising two Japanese species (C. verbosus, C. kibunensis).

A novel transcription factor, ELF5, belongs to the ELF subfamily of ETS genes and maps to human chromosome 11p13-15, a region subject to LOH and rearrangement in human carcinoma cell lines.

The ETS transcription factors are a large family implicated in the control of cellular proliferation and tumorigenesis. In addition, chromosomal translocations involving ETS family members are associated with a range of different human cancers. Given the extensive involvement of ETS factors in tumorigenesis, it becomes important to identify any additional ETS genes that may also play oncogenic roles. We identify a novel gene, ELF5, that appears to belong to the ELF (E74-like-factor) subfamily of the ETS transcription factor family, based upon similarity within the 'ETS domain'. ELF5 displays a similar, but more restricted, expression pattern to that of the newly isolated epithelium-specific ETS gene, ELF3. Unlike most other ETS family members, ELF5 is not expressed in hematopoietic compartments, but is restricted to organs such as lung, stomach, kidney, prostate, bladder and mammary gland. ELF5 is localized to human chromosome 11p13-15, a region that frequently undergoes loss of heterozygosity (LOH) in several types of carcinoma, including those of breast, kidney and prostate. We find that ELF5 expression is not detectable in a number of carcinoma cell lines, some of which display loss or rearrangement of an ELF5 allele. Similar to other ETS family members, ELF5 displays specific binding to DNA sequences containing a GGAA-core. In addition, ELF5 is able to transactivate through these ETS sequences, present upstream from a minimal promoter. Our data suggest that ELF5 may play roles in mammary, lung, prostate and/or kidney function, and possibly also in tumorigenesis.

INSECT PHYLOGENOMICS. Response to Comment on "Phylogenomics resolves the timing and pattern of insect evolution".

Tong et al. comment on the accuracy of the dating analysis presented in our work on the phylogeny of insects and provide a reanalysis of our data. They replace log-normal priors with uniform priors and add a "roachoid" fossil as a calibration point. Although the reanalysis provides an interesting alternative viewpoint, we maintain that our choices were appropriate.

Polymorphism of phase II enzymes: identification of new enzymes and polymorphic variants by database analysis.

The Phase II enzymes of xenobiotic metabolism are characterized by a high level of substrate diversity and genetic polymorphism. Genetic polymorphism of the Phase II enzymes can be of substantial clinical significance as some variants have differences in substrate specificity, stability and levels of expression. Variation in these factors can give rise to abnormal drug metabolism and susceptibility to carcinogens and toxins. A new approach to the discovery of additional members of Phase II enzyme families and the identification of polymorphic variants using searches of the EST databases has been investigated. The examples provided demonstrate that relatively simple search strategies can be highly productive.

The cytochrome b region in the mitochondrial DNA of the ant Tetraponera rufoniger: sequence divergence in Hymenoptera may be associated with nucleotide content.

Polymerase chain reaction (PCR) followed by sequencing of single-stranded DNA yielded sequence information from the cytochrome b (cyt b) region in mitochondrial DNA from the ant Tetraponera rufoniger. Compared with the cyt b genes from Apis mellifera, Drosophila melanogaster, and D. yakuba, the overall A+T content (A+T%) of that of T. rufoniger is lower (69.9% vs 80.7%, 74.2%, and 73.9%, respectively) than those of the other three. The codon usage in the cyt b gene of T. rufoniger is biased although not as much as in A. mellifera, D. melanogaster, and D. yakuba; T. rufoniger has eight unused codons whereas D. melanogaster, D. yakuba, and A. mellifera have 21, 20, and 23, respectively. The inferred cyt b polypeptide chain (PPC) of T. rufoniger has diverged at least as much from a common ancestor with D. yakuba as has that of A. mellifera (approximately 3.5 vs approximately 2.9). Despite the lower A+T%, the relative frequencies of amino acids in the cyt b PPC of T. rufoniger are significantly (P < 0.05) associated with the content of adenine and thymine (A+T%) and size of codon families. The mitochondrially located cytochrome oxidase subunit II genes (CO-II) of endopterygote insects have significantly higher average A+T% (approximately 75%) than those of exopterygous (approximately 69%) and paleopterous (approximately 69%) insects. The increase in A+T% of endopterygote insects occurred in Upper Carboniferous and coincided with a significant acceleration of PPC divergence. However, acceleration of PPC divergence is not significantly correlated with the increase of the A+T% (P > 0.1). The high A+T%, the biased codon usage, and the increased PPC divergence of Hymenoptera can in that respect most easily be explained by directional mutation pressure which began in the Upper Carboniferous and still occurs in most members of the order. Given the roughly identical A+T% of the cyt b and CO-II genes from the other insects whose DNA sequences are known (A. mellifera, D. melanogaster, and D. yakuba), it seems most likely that the A+T% of T. rufoniger declined secondarily within the last 100 Myr as a result of a reduced directional mutation pressure.

A molecular framework for the phylogeny of the ant subfamily dolichoderinae.

Partial sequences are reported for the mitochondrial genes for cytochrome oxidase subunits 2 and 3 and for cytochrome b, and the entire sequence of the gene for tRNA(Leu)(UUR) for species from 14 genera of dolichoderine ants and from three outgroup genera. Considerable variation was observed between tRNA genes in the size of the TPsiC arm and the DHU and anticodon loops and whether or not the TPsiC stem possesses a GC pair. The outgroup taxa showed complete TAA CO1 stop codons, but dolichoderines have either TA or T. The outgroup taxa showed a noncoding gap between the CO1 and the tRNA(Leu)(UUR) genes. A phylogeny-independent compatibility test using the amino acid sequences showed differences between the genes consistent with variation in evolutionary rates, according with other studies. Base compositions proved heterogeneous between species, hence phylogenetic analysis was restricted to the protein sequences using maximum likelihood and the mtREV24 replacement matrix. A maximum-likelihood consensus tree has similarities to those from morphological studies with some exceptions such Leptomyrmex falling within the dolichoderine genera rather than basally, and the accretion of genera formerly included under Iridomyrmex. Features of the tRNA genes and the CO1 termination codons agree quite well with the molecular phylogeny.

Accelerated evolution of cytochrome b in simian primates: adaptive evolution in concert with other mitochondrial proteins?

We have sequenced the cytochrome b gene of Horsfield's tarsier, Tarsius bancanus, to complete a data set of sequences for this gene from representatives of each primate infraorder. These primate cytochrome b sequences were combined with those from representatives of three other mammalian orders (cat, whale, and rat) in an analysis of relative evolutionary rates. The nonsynonymous nucleotide substitution rate of the cytochrome b gene has increased approximately twofold along lineages leading to simian primates compared to that of the tarsier and other primate and nonprimate mammalian species. However, the rate of transversional substitutions at fourfold degenerate sites has remained uniform among all lineages. This increase in the evolutionary rate of cytochrome b is similar in character and magnitude to that described previously for the cytochrome c oxidase subunit II gene. We propose that the evolutionary rate increase observed for cytochrome b and cytochrome c oxidase subunit II may underlie an episode of coadaptive evolution of these two proteins in the mitochondria of simian primates.

Nucleotide composition bias affects amino acid content in proteins coded by animal mitochondria.

We show that in animal mitochondria homologous genes that differ in guanine plus cytosine (G + C) content code for proteins differing in amino acid content in a manner that relates to the G + C content of the codons. DNA sequences were analyzed using square plots, a new method that combines graphical visualization and statistical analysis of compositional differences in both DNA and protein. Square plots divide codons into four groups based on first and second position A + T (adenine plus thymine) and G + C content and indicate differences in amino acid content when comparing sequences that differ in G + C content. When sequences are compared using these plots, the amino acid content is shown to correlate with the nucleotide bias of the genes. This amino acid effect is shown in all protein-coding genes in the mitochondrial genome, including cox I, cox II, and cyt b, mitochondrial genes which are commonly used for phylogenetic studies. Furthermore, nucleotide content differences are shown to affect the content of all amino acids with A + T- and G + C-rich codons. We speculate that phylogenetic analysis of genes so affected may tend erroneously to indicate relatedness (or lack thereof) based only on amino acid content.

TrExML: a maximum-likelihood approach for extensive tree-space exploration.

MOTIVATION: Maximum-likelihood analysis of nucleotide and amino acid sequences is a powerful approach for inferring phylogenetic relationships and for comparing evolutionary hypotheses. Because it is a computationally demanding and time-consuming process, most algorithms explore only a minute portion of tree-space, with the emphasis on finding the most likely tree while ignoring the less likely, but not significantly worse, trees. However, when such trees exist, it is equally important to identify them to give due consideration to the phylogenetic uncertainty. Consequently, it is necessary to change the focus of these algorithms such that near optimal trees are also identified. RESULTS: This paper presents the Advanced Stepwise Addition Algorithm for exploring tree-space and two algorithms for generating all binary trees on a set of sequences. The Advanced Stepwise Addition Algorithm has been implemented in TrExML, a phylogenetic program for maximum-likelihood analysis of nucleotide sequences. TrExML is shown to be more effective at finding near optimal trees than a similar program, fastDNAml, implying that TrExML offers a better approach to account for phylogenetic uncertainty than has previously been possible. A program, TreeGen, is also described; it generates binary trees on a set of sequences allowing for extensive exploration of tree-space using other programs. AVAILABILITY: TreeGen, TrExML, and the sequence data used to test the programs are available from the following two WWW sites: http://whitetail.bemidji.msus. edu/trexml/and http://jcsmr.anu.edu.au/dmm/humgen.+ ++html.

Analysis of directional mutation pressure and nucleotide content in mitochondrial cytochrome b genes.

We present a new approach for analyzing directional mutation pressure and nucleotide content in protein-coding genes. Directional mutation pressure, the heterogenicity in the likelihood of different nucleotide substitutions, is used to explain the increasing or decreasing guanine-cytosine content (GC%) in DNA and is represented by microD, in agreement with Sueoka (1962, Proc Natl Acad Sci USA 48:582-592). The new method uses simulation to facilitate identification of significant A+T or G+C pressure as well as the comparison of directional mutation pressure among genes, even when they are translated by different genetic codes. We use the method to analyze the evolution of directional mutation pressure and nucleotide content of mitochondrial cytochrome b genes. Results from a survey of 110 taxa indicate that the cytochrome b genes of most taxa are subjected to significant directional mutation pressure and that the gene is subject to A+T pressure in most cases. Only in the anseriform bird Cairina moschata is the cytochrome b gene subject to significant G+C pressure. The GC% at nonsynonymous codon sites decreases proportionately with increasing A+T pressure, and with a slope less than one, indicating a presence of selective constraints. The cytochrome b genes of insects, nematodes, and eumycotes are subject to extreme A+T pressures (microD = 0.123, 0.224, and 0.130) and, in parallel, the GC% of the nonsynonymous codon sites has decreased from about 0.44 in organisms that are not subjected to A+T or G+C pressure to about 0.332, 0.323, and 0.367, respectively. The distribution of taxa according to the GC% at nonsynonymous codon sites and directional mutation pressure supports the notion that variation in these parameters is a phylogenetic component.

Zeta, a novel class of glutathione transferases in a range of species from plants to humans.

Sequence alignment and phylogenetic analysis has identified a new subgroup of glutathione S-transferase (GST)-like proteins from a range of species extending from plants to humans. This group has been termed the Zeta class. An atomic model of the N-terminal domain suggests that the members of the Zeta class have a similar structure to that of other GSTs, binding glutathione in a similar orientation in the G site. Recombinant human GSTZ1-1 has been expressed in Escherichia coli and characterized. The protein is a dimer composed of 24.2 kDa subunits and has minimal glutathione-conjugating activity with ethacrynic acid and 7-chloro-4-nitrobenz-2-oxa-1, 3-diazole. Although low in comparison with other GSTs, GSTZ1-1 has glutathione peroxidase activity with t-butyl and cumene hydroperoxides. The members of the Zeta class have been conserved over a long evolutionary period, suggesting that they might have a role in the metabolism of a compound that is common in many living cells.

Amylase and 16S rRNA genes from a hyperthermophilic archaebacterium.

A hyperthermophilic and amylolytic prokaryote, designated Rt3, was isolated from a thermal spring near Rotorua, New Zealand. The 16S rRNA gene of Rt3 was cloned and sequenced with the aim of determining its phylogenetic affiliations. The phylogenetic analysis of this sequence, which included a selection of archaebacterial and eubacterial 16S rRNA sequences, indicates that Rt3 most likely belongs to the archaebacterial order Thermococcales. An amylase gene (amyA) from Rt3, encoding a highly thermostable amylase activity, was cloned and its DNA sequence determined. Transcriptional signals typical of archaebacteria were evident in this sequence. The sequence is homologous to a broad range of enzymes from the AMY superfamily and contains a typical N-terminal signal peptide. Phylogenetic analysis and comparison of structural features with other AMY superfamily enzymes reveals that, firstly, the closest homologues of the Rt3 amylase are members of the Bacillus and Plant alpha-amylase groups; and secondly, that the Rt3 amylase is closely related to only one other currently known archaebacterial enzyme, i.e. an (AMY superfamily) alpha-amylase from Natronococcus.

Unbiased estimation of symmetrical directional mutation pressure from protein-coding DNA.

The most generally applicable procedure for obtaining estimates of the symmetrical, or strandnonspecific, directional mutation pressure (microD) on protein-coding DNA sequences is to determine the G+C content at synonymous codon sites (Psyn), and to divide Psyn by twice the arithmetic mean of the G+C content at synonymous codon sites of a large number of randomly generated, synonymously coding DNA sequences (Psyn). Unfortunately, the original procedure yields biased estimates of Psyn and microD and is computationally expensive. We here present a fast procedure for estimating unbiased microD values. The procedure employs direct calculation of Psyn (approximately Psyn) and two normalization procedures, one for Psyn < or = Psyn and another for Psyn > or = Psyn. The normalization removes a bias sometimes caused by codons specifying arginine, asparagine, isoleucine, and leucine. Consequently, comparison of protein-coding genes that are translated using different genetic codes is facilitated.

Mitochondrial DNA sequences in ancient Australians: Implications for modern human origins.

DNA from ancient human remains provides perspectives on the origin of our species and the relationship between molecular and morphological variation. We report analysis of mtDNA from the remains of 10 ancient Australians. These include the morphologically gracile Lake Mungo 3 [ approximately 60 thousand years (ka) before present] and three other gracile individuals from Holocene deposits at Willandra Lakes (<10 ka), all within the skeletal range of living Australians, and six Pleistocene/early Holocene individuals (15 to <8 ka) from Kow Swamp with robust morphologies outside the skeletal range of contemporary indigenous Australians. Lake Mungo 3 is the oldest (Pleistocene) "anatomically modern" human from whom DNA has been recovered. His mtDNA belonged to a lineage that only survives as a segment inserted into chromosome 11 of the nuclear genome, which is now widespread among human populations. This lineage probably diverged before the most recent common ancestor of contemporary human mitochondrial genomes. This timing of divergence implies that the deepest known mtDNA lineage from an anatomically modern human occurred in Australia; analysis restricted to living humans places the deepest branches in East Africa. The other ancient Australian individuals we examined have mtDNA sequences descended from the most recent common ancestor of living humans. Our results indicate that anatomically modern humans were present in Australia before the complete fixation of the mtDNA lineage now found in all living people. Sequences from additional ancient humans may further challenge current concepts of modern human origins.

Identification of novel glutathione transferases and polymorphic variants by expressed sequence tag database analysis.

The human expressed sequence tag (EST) database can be searched by different sequence alignment strategies to identify new members of gene families and allelic variants. To illustrate the value of database analysis for gene discovery, we have focused on the glutathione S-transferase (GST) super family, an approach that has led to the identification of the Zeta class. The Zeta class GSTs catalyze the glutathione-dependent biotransformation of alpha-haloacids and the isomerization of maleylacetoacetic acid to fumarylacetoacetic acid, an essential step in the catabolism of tyrosine. Allelic variants of the GST Z1 and GST A2 genes have also been identified by EST database analysis. One GST Z1 variant (GST Z1A) has significantly higher activity with dichloroacetic acid as a substrate than other GST Z1 isoforms. This variant may be important in the clinical treatment of lactic acidosis where dichloroacetic acid is prescribed. Our experience with the application of EST database searching methods suggests that it may be productively applied to other gene families of pharmacogenetic interest.

Gene targeting of Desrt, a novel ARID class DNA-binding protein, causes growth retardation and abnormal development of reproductive organs.

We have cloned and characterized a novel murine DNA-binding protein Desrt, with a motif characteristic of the ARID (A-T rich interaction domain) family of transcription factors. The Desrt gene encodes an 83-kD protein that is shown to bind DNA and is widely expressed in adult tissues. To examine the in vivo function of Desrt, we have generated mice with a targeted mutation in the ARID domain of Desrt. Homozygous mutants have reduced viability, pronounced growth retardation, and a high incidence of abnormalities of the female and male reproductive organs including cryptorchidism. This may thus serve as a model to dissect the mechanisms involved in the development of the reproductive tract including testicular descent. Gene-targeted mice also display a reduction in the thickness of the zona reticularis of the adrenal gland and transient aberrations of the T and B cell compartments of primary lymphoid organs. These data show that this novel DNA-binding protein, Desrt, has a nonredundant function during growth and in the development of the reproductive system.

Many parallel losses of infA from chloroplast DNA during angiosperm evolution with multiple independent transfers to the nucleus.

We used DNA sequencing and gel blot surveys to assess the integrity of the chloroplast gene infA, which codes for translation initiation factor 1, in >300 diverse angiosperms. Whereas most angiosperms appear to contain an intact chloroplast infA gene, the gene has repeatedly become defunct in approximately 24 separate lineages of angiosperms, including almost all rosid species. In four species in which chloroplast infA is defunct, transferred and expressed copies of the gene were found in the nucleus, complete with putative chloroplast transit peptide sequences. The transit peptide sequences of the nuclear infA genes from soybean and Arabidopsis were shown to be functional by their ability to target green fluorescent protein to chloroplasts in vivo. Phylogenetic analysis of infA sequences and assessment of transit peptide homology indicate that the four nuclear infA genes are probably derived from four independent gene transfers from chloroplast to nuclear DNA during angiosperm evolution. Considering this and the many separate losses of infA from chloroplast DNA, the gene has probably been transferred many more times, making infA by far the most mobile chloroplast gene known in plants.