Juvenile-mature wood transition in pine: correlation between wood properties and candidate gene expression profiles

There is a strong desire to exploit transcriptomics data from model species for the genetic improvement of non-model crops. Here, we use gene expression profiles from the commercial model Pinus taeda to identify candidate genes implicated in juvenile-mature wood transition in the non-model relative, P. sylvestris. Re-analysis of 'public domain' SAGE data from xylem tissues of P. taeda revealed 283 mature-abundant and 396 juvenile-abundant tags (P < 0.01), of which 70 and 137, respectively matched to genes with known function. Based on sequence similarity, we then isolated 16 putative homologues of genes that in P. taeda exhibited widest divergence in expression between juvenile and mature samples. Candidate expression levels in P. sylvestris were almost invariably differential between juvenile and mature woody tissue samples among two cohorts of five trees collected from the same seed source and selected for genetic uniformity by genetic distance analysis. However, the direction of differential expression was not always consistent with that described in the original P. taeda SAGE data. Correlation was observed between gene expression and juvenile-mature wood anatomical characteristics by OPLS analysis. Four candidates (alpha-tubulin, porin MIP1, lipid transfer protein and aquaporin like protein) apparently had greatest influence on the wood traits measured. Speculative function of these genes in relation to juvenile-mature wood transition is briefly explored. Thus, we demonstrate the feasibility of exploiting SAGE data from a model species to identify consistently differentially expressed candidates in a related non-model species.

Using phylogeny to investigate the origins of the Cape flora: the importance of taxonomic, gene and genome sampling strategies

Phylogenetic methods hold great promise for the reconstruction of the transition from precursor to modern flora and the identification of underlying factors which drive the process. The phylogenetic methods presently used to address the question of the origin of the Cape flora of South Africa are considered here. The sampling requirements of each of these methods, which include dating of diversifications using calibrated molecular trees, sister pair comparisons, lineage through time plots and biogeographical optimizations are reviewed. Sampling of genes, genomes and species are considered. Although increased higher-level studies and increased sampling are required for robust interpretation, it is clear that much progress is already made. It is argued that despite the remarkable richness of the flora, the Cape flora is a valuable model system to demonstrate the utility of phylogenetic methods in determining the history of a modern flora.

Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.

A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data

We describe a general likelihood-based 'mixture model' for inferring phylogenetic trees from gene-sequence or other character-state data. The model accommodates cases in which different sites in the alignment evolve in qualitatively distinct ways, but does not require prior knowledge of these patterns or partitioning of the data. We call this qualitative variability in the pattern of evolution across sites "pattern-heterogeneity" to distinguish it from both a homogenous process of evolution and from one characterized principally by differences in rates of evolution. We present studies to show that the model correctly retrieves the signals of pattern-heterogeneity from simulated gene-sequence data, and we apply the method to protein-coding genes and to a ribosomal 12S data set. The mixture model outperforms conventional partitioning in both these data sets. We implement the mixture model such that it can simultaneously detect rate- and pattern-heterogeneity. The model simplifies to a homogeneous model or a rate- variability model as special cases, and therefore always performs at least as well as these two approaches, and often considerably improves upon them. We make the model available within a Bayesian Markov-chain Monte Carlo framework for phylogenetic inference, as an easy-to-use computer program.

DNA variation in a 13-Mb region including the F9 gene: inferring the genealogical history and causal role of a hemophilia B mutation (IVS 5+13 A -> G)

About 5.5% of all UK hemophilia B patients have the base substitution IVS 5+13 A-->G as the only change in their factor (F)IX gene (F9). This generates a novel donor splice site which fits the consensus better than the normal intron 5 donor splice. Use of the novel splice site should result in a missense mutation followed by the abnormal addition of four amino acids to the patients' FIX. In order to explain the prevalence of this mutation, its genealogical history is examined. Analysis of restriction fragment length polymorphism in the 21 reference UK individuals (from different families) with the above mutation showed identical haplotypes in 19 while two differed from the rest and from each other. In order to investigate the history of the mutation and to verify that it had occurred independently more than once, the sequence variation in 1.5-kb segments scattered over a 13-Mb region including F9 was examined in 18 patients and 15 controls. This variation was then analyzed with a recently developed Bayesian approach that reconstructs the genealogy of the gene investigated while providing evidence of independent mutations that contribute disconnected branches to the genealogical tree. The method also provides minimum estimates of the age of the mutation inherited by the members of coherent trees. This revealed that 17 or 18 mutant genes descend from a founder who probably lived 450 years ago, while one patient carries an independent mutation. The independent recurrence of the IVS5+13 A-->G mutation strongly supports the conclusion that it is the cause of these patients' mild hemophilia.

Isolation and characterization of an AINTEGUMENTA-like gene in different coconut (Cocos nucifera L.) varieties from Sri Lanka

Development of an efficient tissue culture protocol in coconut is hampered by numerous technical constraints. Thus a greater understanding of the fundamental aspects of embryogenesis is essential. The role of AINTEGUMENTA-like genes in embryogenesis has been elucidated not only in model plants but also in economically important crops. A coconut gene, CnANT, that encodes two APETALA2 (AP2) domains and a conserved linker region similar to those of the BABY BOOM transcription factor was cloned, characterized, and its tissue specific expression was examined. The full-length cDNA of 1,780 bp contains a 1,425-bp open reading frame that encodes a putative peptide of 474 amino acids. The genomic DNA sequence includes 2,317 bp and consists of nine exons interrupted by eight introns. The exon/intron organization of CnANT is similar to that of homologous genes in other plant species. Analysis of differential tissue expression by real-time polymerase chain reaction indicated that CnANT is expressed more highly in in vitro grown tissues than in other vegetative tissues. Sequence comparison of the genomic sequence of CnANT in different coconut varieties revealed one single nucleotide polymorphism and one indel in the first exon and first intron, respectively, which differentiate the Tall group of trees from Dwarfs. The indel sequence, which can be considered a simple sequence repeats marker, was successfully used to distinguish the Tall and Dwarf groups as well as to develop a marker system, which may be of value in the identification of parental varieties that are used in coconut breeding programs in Sri Lanka.

Oscheius tipulae as an example of eEF1A gene diversity in nematodes

We characterized four eEF1A genes in the alternative rhabditid nematode model organism Oscheius tipulae. This is twice the copy number of eEF1A genes in C. elegans, C. briggsae, and, probably, many other free-living and parasitic nematodes. The introns show features remarkably different from those of other metazoan eEF1A genes. Most of the introns in the eEF1A genes are specific to O. tipulae and are not shared with any of the other genes described in metazoans. Most of the introns are phase 0 (inserted between two codons), and few are inserted in protosplice sites (introns inserted between the nucleotide sequence A/CAG and G/A). Two of these phase 0 introns are conserved in sequence in two or more of the four eEF1A gene copies, and are inserted in the same position in the genes. Neither of these characteristics has been detected in any of the nematode eEF1A genes characterized to date. The coding sequences were also compared with other eEF1A cDNAs from 11 different nematodes to determine the variability of these genes within the phylum Nematoda. Parsimony and distance trees yielded similar topologies, which were similar to those created using other molecular markers. The presence of more than one copy of the eEF1A gene with nearly identical coding regions makes it difficult to define the orthologous cDNAs. As shown by our data on O. tipulae, careful and extensive examination of intron positions in the eEF1A gene across the phylum is necessary to define their potential for use as valid phylogenetic markers.

Molecular characterization of nitrate reductase gene and its expression in the marine red alga Gracilaria tenuistipitata (Rhodophyta)

The enzyme nitrate reductase (NR) responsible for the conversion of nitrate to nitrite is considered to be the rate-limiting step in nitrogen assimilation. The economically important marine macroalga Gracilaria tenuistipitata presents a circadian oscillation in NR protein content and activity. In order to identify if the regulation of NR in G. tenuistipitata happens at transcriptional levels, the NR cDNA and gene were sequenced and the NR mRNA expression was studied. Analysis of the sequenced gene revealed absence of introns which is unusual for NR genes. The transcriptional profiling revealed a circadian rhythm for NR; furthermore, a rhythm was observed in constant light condition, suggesting a possible regulation by the biological clock at the mRNA levels for NR in G. tenuistipitata.

Effects of the antimicrobial peptide gomesin on the global gene expression profile, virulence and biofilm formation of Xylella fastidiosa

In the xylem vessels of susceptible hosts, such as citrus trees, Xylella fastidiosa forms biofilm-like colonies that can block water transport, which appears to correlate to disease symptoms. Besides aiding host colonization, bacterial biofilms play an important role in resistance against antimicrobial agents, for instance antimicrobial peptides (AMPs). Here, we show that gomesin, a potent AMP from a tarantula spider, modulates X. fastidiosa gene expression profile upon 60 min of treatment with a sublethal concentration. DNA microarray hybridizations revealed that among the upregulated coding sequences, some are related to biofilm production. In addition, we show that the biofilm formed by gomesin-treated bacteria is thicker than that formed by nontreated cells or cells exposed to streptomycin. We have also observed that the treatment of X. fastidiosa with a sublethal concentration of gomesin before inoculation in tobacco plants correlates with a reduction in foliar symptoms, an effect possibly due to the trapping of bacterial cells to fewer xylem vessels, given the enhancement in biofilm production. These results warrant further investigation of how X. fastidiosa would respond to the AMPs produced by citrus endophytes and by the insect vector, leading to a better understanding of the mechanism of action of these molecules on bacterial virulence.

Amplification of multiple copies of mitochondrial Cytochrome b gene fragments in the Australian freshwater crayfish, Cherax destructor Clark (Parastacidae: Decapoda)

Non-coding copies of fragments of the mitochondrial genome translocated to the nucleus or pseudogenes are being found with increasing frequency in a diversity of organisms. As part of a study to evaluate the utility of a range of mitochondrial gene regions for population genetic and systematic studies of the Australian freshwater crayfish, Cherax destructor (the yabby), we report the first detection of Cytochrome b (Cyt b) pseudogenes in crustaceans. We amplified and sequenced fragments of the mitochondrial Cyt b gene from 14 individuals of C. destructor using polymerase chain reaction (PCR) with primers designed from conserved regions of Penaeus monodon and Drosophila melanogaster mitochondrial genomes. The phylogenetic tree produced from the amplified fragments using these primers showed a very different topology to the trees obtained from sequences from three other mitochondrial genes, suggesting one or more nuclear pseudogenes have been amplified. Supporting this conclusion, two highly divergent sequences were isolated from each of two single individuals, and a 2 base pair (bp) deletion in one sequence was observed. There was no evidence to support inadvertent amplification of parasite DNA or contamination of samples from other sources. These results add to other recent observations of pseudogenes suggesting the frequent transfer of mitochondrial DNA (mtDNA) genes to the nucleus and reinforces the necessity of great care in interpreting PCR-generated Cyt b sequences used in population or evolutionary studies in freshwater crayfish and crustaceans more generally.

Phylogenetic relationships among nine scallop species (Bivalvia : Pectinidae) inferred from nucleotide sequences of one mitochondrial and three nuclear gene regions

Current knowledge of the evolutionary relationships among scallop species (Mollusca: Bivalvia: Pectinidae) in the Indo-Pacific region is rather scanty. To enhance the understanding of the relationships within this group, phylogenies of nine species of scallops with the majority from coastal regions of Thailand, were reconstructed by maximum parsimony, maximum likelihood, and Bayesian methods using sequences of the 16S rRNA of the mitochondrial genome, and a fragment containing the ITS1, 5.8S and ITS2 genes of the nuclear DNA. The trees that resulted from the three methods of analysis were topologically identical, however, gained different levels of support at some nodes. Nine species were clustered into two major clades, corresponding to two subfamilies (Pectininae and Chlamydinae) of the three currently recognized subfamilies within Pectinidae. Overall, the relationships reported herein are mostly in accordance with the previous molecular studies that used sequences of the mtDNA cytochrome oxidase subunit I, and the classification system based on microsculpture of shell features and morphological characteristics of juveniles. Levels of divergences were different among genes (i.e., the 5.8S gene showed the lowest levels of nucleotide divergence at all levels, whereas the 16S rRNA showed the highest level of variation within species, and ITS2 gene revealed the highest level of divergence at higher levels).

Low levels of realized seed and pollen gene flow and strong spatial genetic structure in a small, isolated and fragmented population of the tropical tree Copaifera langsdorffii Desf

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Pollen Dispersal Between Isolated Trees in the Brazilian Savannah: A Case Study of the Neotropical Tree Hymenaea stigonocarpa

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Análise do gene E6 de HPV de baixo risco em papilomatose de laringe

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Molecular phylogeny of the genus Saguinus (Platyrrhini, Primates) based on the ND1 mitochondrial gene and implications for conservation

The systematics of the subfamily Callitrichinae (Platyrrhini, Primates), a group of small monkeys from South America and Panama, remains an area of considerable discussion despite many investigations, there being continuing controversy over subgeneric taxonomic classifications based on morphological characters. The purpose of our research was to help elucidate the phylogenetic relationships within the monkey genus Saguinus (Callitrichinae) using a molecular approach to discover whether or not the two different sections containing hairy-faced and bare-faced species are monophyletic, whether Saguinus midas midas and Saguinus bicolor are more closely related than are S. midas midas and Saguinus midas niger, and if Saguinus fuscicollis melanoleucus and Saguinus fuscicollis weddelli really are different species. We sequenced the 957 bp ND1 mitochondrial gene of 21 Saguinus monkeys (belonging to six species and nine morphotypes) and one Cebus monkey (the outgroup) and constructed phylogenetic trees using maximum parsimony, neighbor joining, and maximum likelihood methods. The phylogenetic trees obtained divided the genus Saguinus into two groups, one containing the small-bodied species S. fuscicollis and the other, the large-bodied species S. mystax, S. leucopus, S. oedipus, S. midas, S. bicolor. The most derived taxa, S. midas and S. bicolor, grouped together, while S. fuscicollis melanoleucus and S. f. weddelli showed divergence values that did not support the division of these morphotypes into subspecies. On the other hand, S. midas individuals showed divergence compatible with the existence of three subspecies, two of them with the same morphotype as the subspecies S. midas niger. The results of our study suggest that there is at least one Saguinus subspecies that has not yet been described and that the conservation status of Saguinus species and subspecies should be carefully revised using modern molecular approaches.