Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores

Darwin's paradigm holds that the diversity of present-day organisms has arisen via a process of genetic descent with modification, as on a bifurcating tree. Evidence is accumulating that genes are sometimes transferred not along lineages but rather across lineages. To the extent that this is so, Darwin's paradigm can apply only imperfectly to genomes, potentially complicating or perhaps undermining attempts to reconstruct historical relationships among genomes (i.e., a genome tree). Whether most genes in a genome have arisen via treelike (vertical) descent or by lateral transfer across lineages can be tested if enough complete genome sequences are used. We define a phylogenetically discordant sequence (PDS) as an open reading frame (ORF) that exhibits patterns of similarity relationships statistically distinguishable from those of most other ORFs in the same genome. PDSs represent between 6.0 and 16.8% (mean, 10.8%) of the analyzable ORFs in the genomes of 28 bacteria, eight archaea, and one eukaryote (Saccharomyces cerevisiae). In this study we developed and assessed a distance-based approach, based on mean pairwise sequence similarity, for generating genome trees. Exclusion of PDSs improved bootstrap support for basal nodes but altered few topological features, indicating that there is little systematic bias among PDSs. Many but not all features of the genome tree from which PDSs were excluded are consistent with the 16S rRNA tree.

Distributional profiles of homologous open reading frames among bacterial phyla: implications for vertical and lateral transmission

If open reading frames (ORFs) have been transmitted primarily by vertical descent, the distributional profile of orthologues of each ORF should be congruent with the organismal tree or a subtree thereof. Distributional patterns not reconciled parsimoniously with tree-like descent and loss are prima facie evidence of lateral gene transfer. Herein, a rigorous criterion for recognizing ORF distributions is described and implemented; it does not require the inference of phylogenetic trees, nor does it assume any specific tree. Because lineage-specific differences in rates of sequence change can also generate unexpected distributional patterns, rate artefacts, were controlled for by requiring pairwise matches between ORFs to exceed a rigorous inclusion threshold, but absence of a match was assessed against a more-permissive exclusion threshold. Applying this dual-threshold criterion to cross-domain and cross-phylum distributional patterns for ORFs in 23 bacterial genomes, a relative abundance of ORFs was observed that find a match in exactly seven other bacterial phyla; 94-99% of these ORFs also find matches among the Archaea and/or Eukarya. In the larger (and some smaller) bacterial genomes, ORFs that find matches in exactly one other bacterial phylum are also relatively abundant, but fewer of these have non-bacterial homologues; most of their matches within the Bacteria are to the Proteobacteria and/or Firmicutes, which cannot be sister lineages to all bacteria. ORFs that are neither distributed universally among the Bacteria, nor necessarily shared with topologically adjacent lineages, are preferentially enriched in large bacterial genomes.

Candidatus "Scalindua brodae", spec. nov., Candidatus "Scalindua wagneri", spec. nov., two new species of anaerobic ammonium oxidizing bacteria

Anaerobic ammonium oxidation (anammox) is both a promising process in wastewater treatment and a long overlooked microbial physiology that can contribute significantly to biological nitrogen cycling in the world's oceans. Anammox is mediated by a monophyletic group of bacteria that branches deeply in the Planctomycetales. Here we describe a new genus and species of anaerobic ammonium oxidizing planctomycetes, discovered in a wastewater treatment plant (wwtp) treating landfill leachate in Pitsea, UK. The biomass from this wwtp showed high anammox activity (5.0 +/- 0.5 nmol/mg protein/min) and produced hydrazine from hydroxylamine, one of the unique features of anammox bacteria. Eight new planctomycete 16S rRNA gene sequences were present in the 16S rRNA gene clone library generated from the biomass. Four of these were affiliated to known anammox 16S rRNA gene sequences, but branched much closer to the root of the planctomycete line of descent. Fluorescence in situ hybridization (FISH) with oligonucleotide probes specific for these new sequences showed that two species (belonging to the same genus) together made up > 99% of the planctomycete population which constituted 20% of the total microbial community. The identification of these organisms as typical anammox bacteria was confirmed with electron microscopy and lipid analysis. The new species, provisionally named Candidatus Scalindua brodae and Scalindua wagneri considerably extend the biodiversity of the anammox lineage on the 16S rRNA gene level, but otherwise resemble known anammox bacteria. Simultaneously, another new species of the same genus, Candidatus Scalindua sorokinii, was detected in the water column of the Black Sea, making this genus the most widespread of all anammox bacteria described so far.

Effect of alcohol on iron storage diseases of the liver

The consumption of excess alcohol in patients with liver iron storage diseases, in particular the iron-overload disease hereditary haemochromatosis (HH), has important clinical consequences. HH, a common genetic disorder amongst people of European descent, results in a slow, progressive accumulation of excess hepatic iron. If left untreated, the condition may lead to fibrosis, cirrhosis and primary hepatocellular carcinoma. The consumption of excess alcohol remains an important cause of hepatic cirrhosis and alcohol consumption itself may lead to altered iron homeostasis. Both alcohol and iron independently have been shown to result in increased oxidative stress causing lipid peroxidation and tissue damage. Therefore, the added effects of both toxins may exacerbate the pathogenesis of disease and impose an increased risk of cirrhosis. This review discusses the concomitant effects of alcohol and iron on the pathogenesis of liver disease. We also discuss the implications of co-existent alcohol and iron in end-stage liver disease.

Prospects for whole genome linkage disequilibrium mapping in domestic dog breeds

Linkage disequilibrium (LD) mapping is commonly used as a fine mapping tool in human genome mapping and has been used with some success for initial disease gene isolation in certain isolated inbred human populations. An understanding of the population history of domestic dog breeds suggests that LID mapping could be routinely utilized in this species for initial genome-wide scans. Such an approach offers significant advantages over traditional linkage analysis. Here, we demonstrate, using canine copper toxicosis in the Bedlington terrier as the model, that LID mapping could be reasonably expected to be a useful strategy in low-resolution, genome-wide scans in pure-bred dogs. Significant LID was demonstrated over distances up to 33.3 cM. It is very unlikely, for a number of reasons discussed, that this result could be extrapolated to the rest of the genome. It is, however, consistent with the expectation given the population structure of canine breeds and, in this breed at least, with the hypothesis that it may be possible to utilize LID in a genome-wide scan. In this study, LD mapping confirmed the location of the copper toxicosis in Bedlington terrier gene (CT-BT) and was able to do so in a population that was refractory to traditional linkage analysis.

Genetic covariation of pelvic organ and elbow mobility in twins and their sisters

A range of environmental risk factors, with childbirth the most notable, have been associated with the development of pelvic organ prolapse and urinary incontinence. However, indications of genetic influence (positive family histories, ethnic differences) have prompted research into the heritability of measures of pelvic organ descent and joint mobility, which have also been associated with prolapse and incontinence. Genes appear to influence about half of the variation in these measures and, furthermore, the pelvic organ measures are associated with elbow hyperextension at a phenotypic level (r approximate to .2). We examined these measures in young, nulligravid women to determine if their association is due to a common genetic source. Data were collected from 178 Caucasian female co-twins and non-twin sisters, 50 of whom returned to be retested, which allowed reliability to be estimated and unreliable variance to be isolated in the multivariate analyses. Structural equation modeling was used to estimate genetic associations between latent elbow and bladder mobility factors for which heritabilities were estimated to be 0.80 and 0.64 respectively. The association between these factors appeared to be mediated by common genes (genetic r = .48, non-shared environmental r = -.06), with genes influencing latent elbow mobility accounting for 14% of the variation in latent bladder mobility. We speculate that genes influencing connective tissue structure may underlie this association.

Using identity-by-de scent information in affected sib pairs to increase the efficiency of genetic association studies

T he aim of this study was to determine whether identity-by-descent (IBD) information for affected sib pairs (ASPs) can be used to select a sample of cases for a genetic case-control study which will provide more power for detecting association with loci in a known linkage region. By modeling the expected frequency of the disease allele in ASPs showing IBD sharing of 0, 1, or 2 alleles, and considering additive, recessive, and dominant disease models, we show that cases selected from IBD 2 families are best for this purpose, followed by those selected from IBD 1 families; least useful are cases selected from IBD 0 families.

Simulating the effects of dominance and epistasis on selection response in the CIMMYT wheat breeding program using QuCim

Plant breeders use many different breeding methods to develop superior cultivars. However, it is difficult, cumbersome, and expensive to evaluate the performance of a breeding method or to compare the efficiencies of different breeding methods within an ongoing breeding program. To facilitate comparisons, we developed a QU-GENE module called QuCim that can simulate a large number of breeding strategies for self-pollinated species. The wheat breeding strategy Selected Bulk used by CIMMYT's wheat breeding program was defined in QuCim as an example of how this is done. This selection method was simulated in QuCim to investigate the effects of deviations from the additive genetic model, in the form of dominance and epistasis, on selection outcomes. The simulation results indicate that the partial dominance model does not greatly influence genetic advance compared with the pure additive model. Genetic advance in genetic systems with overdominance and epistasis are slower than when gene effects are purely additive or partially dominant. The additive gene effect is an appropriate indicator of the change in gene frequency following selection when epistasis is absent. In the absence of epistasis, the additive variance decreases rapidly with selection. However, after several cycles of selection it remains relatively fixed when epistasis is present. The variance from partial dominance is relatively small and therefore hard to detect by the covariance among half sibs and the covariance among full sibs. The dominance variance from the overdominance model can be identified successfully, but it does not change significantly, which confirms that overdominance cannot be utilized by an inbred breeding program. QuCim is an effective tool to compare selection strategies and to validate some theories in quantitative genetics.

Petrology and geochemistry of early cretaceous bimodal continental flood volcanism of the NW Etendeka, Namibia. Part 2: Characteristics and petrogenesis of the high-Ti latite and high-Ti and low-Ti voluminous quartz latite eruptives

As a result of their relative concentration towards the respective Atlantic margins, the silicic eruptives of the Parana (Brazil)-Etendeka large igneous province are disproportionately abundant in the Etendeka of Namibia. The NW Etendeka silicic units, dated at similar to132 Ma, occupy the upper stratigraphic levels of the volcanic sequences, restricted to the coastal zone, and comprise three latites and five quartz latites (QL). The large-volume Fria QL is the only low-Ti type. Its trace element and isotopic signatures indicate massive crustal input. The remaining NW Etendeka silicic units are enigmatic high-Ti types, geochemically different from low-Ti types. They exhibit chemical affinities with the temporally overlapping Khumib high-Ti basalt (see Ewart et al. Part 1) and high crystallization temperatures (greater than or equal to980 to 1120degreesC) inferred from augite and pigeonite phenocrysts, both consistent with their evolution from a mafic source. Geochemically, the high-Ti units define three groups, thought genetically related. We test whether these represent independent liquid lines of descent from a common high-Ti mafic parent. Although the recognition of latites reduces the apparent silica gap, difficulty is encountered in fractional crystallization models by the large volumes of two QL units. Numerical modelling does, however, support large-scale open-system fractional crystallization, assimilation of silicic to basaltic materials, and magma mixing, but cannot entirely exclude partial melting processes within the temporally active extensional environment. The fractional crystallization and mixing signatures add to the complexity of these enigmatic and controversial silicic magmas. The existence, however, of temporally and spatially overlapping high-Ti basalts is, in our view, not coincidental and the high-Ti character of the silicic magmas ultimately reflects a mantle signature.

Prevalence of rectocele in young nulliparous women

Background: It is generally assumed that fascial defects in the rectovaginal septum are the result of childbirth. However, rectoceles do occur in women who have never delivered vaginally. Aims: To determine the incidence of rectocele in a cohort of asymptomatic, young nulliparous women. Methods: Observational cohort study on 178 nulliparous caucasian women (aged 18-24) recruited for a twin study of pelvic floor dysfunction. All women were interviewed and examined by translabial ultrasound, supine and after voiding. In 52 women, 3D imaging was obtained and 171 datasets were complete and available for analysis. Ultrasound findings were reviewed for rectovaginal septal integrity by an assessor blinded against interview and demographic data for rectovaginal septal integrity. Results: A discontinuity of the anterior rectal wall with extrusion of rectal mucosa or contents (depth of ! 10 mm) was observed in 21/171 (12%). The depth of this herniation ranged from 10 to 25 mm and was filled with stool (n = 10) or rectal mucosa (n = 11). Defects were associated with a higher BMI (P = 0.049), with the complaint of constipation (P = 0.049) and non-significantly with straining at stool (P = 0.09). Descent of the ampulla to beyond the level of the symphysis pubis without fascial defect, that is, significant perineal relaxation, was observed in 23/171 (13%). Conclusions: Twelve percent of 171 young nulligravid caucasian women showed a defect of the rectovaginal septum. Associations were observed with higher body mass index and a history of constipation. It is hypothesised that in some women defects of the rectovaginal septum and perineal hypermobility may be congenital in nature.

Genetic time-series analysis identifies a major QTL for in vivo alcohol metabolism not predicted by in vitro studies of structural protein polymorphism at the ADH1B or ADH1C loci

After ingestion of a standardized dose of ethanol, alcohol concentrations were assessed, over 3.5 hours from blood (six readings) and breath (10 readings) in a sample of 412 MZ and DZ twins who took part in an Alcohol Challenge Twin Study (ACTS). Nearly all participants were subsequently genotyped on two polymorphic SNPs in the ADH1B and ADH1C loci known to affect in vitro ADH activity. In the DZ pairs, 14 microsatellite markers covering a 20.5 cM region on chromosome 4 that includes the ADH gene family were assessed, Variation in the timed series of autocorrelated blood and breath alcohol readings was studied using a bivariate simplex design. The contribution of a quantitative trait locus (QTL) or QTL's linked to the ADH region was estimated via a mixture of likelihoods weighted by identity-by-descent probabilities. The effects of allelic substitution at the ADH1B and ADH1C loci were estimated in the means part of the model simultaneously with the effects sex and age. There was a major contribution to variance in alcohol metabolism due to a QTL which accounted for about 64% of the additive genetic covariation common to both blood and breath alcohol readings at the first time point. No effects of the ADH1B*47His or ADH1C*349Ile alleles on in vivo metabolism were observed, although these have been shown to have major effects in vitro. This implies that there is a major determinant of variation for in vivo alcohol metabolism in the ADH region that is not accounted for by these polymorphisms. Earlier analyses of these data suggested that alcohol metabolism is related to drinking behavior and imply that this QTL may be protective against alcohol dependence.

Population-based continuous optimization, probabilistic modelling and mean shift

Evolutionary algorithms perform optimization using a population of sample solution points. An interesting development has been to view population-based optimization as the process of evolving an explicit, probabilistic model of the search space. This paper investigates a formal basis for continuous, population-based optimization in terms of a stochastic gradient descent on the Kullback-Leibler divergence between the model probability density and the objective function, represented as an unknown density of assumed form. This leads to an update rule that is related and compared with previous theoretical work, a continuous version of the population-based incremental learning algorithm, and the generalized mean shift clustering framework. Experimental results are presented that demonstrate the dynamics of the new algorithm on a set of simple test problems.

Coupled thermal, structural and vibrational analysis of a hypersonic engine for flight test

This paper describes a relatively simple and quick method for implementing aerodynamic heating models into a finite element code for non-linear transient thermal-structural and thermal-structural-vibrational analyses of a Mach 10 generic HyShot scramjet engine. The thermal-structural-vibrational response of the engine was studied for the descent trajectory from 60 to 26 km. Aerodynamic heating fluxes, as a function of spatial position and time for varying trajectory points, were implemented in the transient heat analysis. Additionally, the combined effect of varying dynamic pressure and thermal loads with altitude was considered. This aero-thermal-structural analysis capability was used to assess the temperature distribution, engine geometry distortion and yielding of the structural material due to aerodynamic heating during the descent trajectory, and for optimising the wall thickness, nose radius of leading edge, etc. of the engine intake. A structural vibration analysis was also performed following the aero-thermal-structural analysis to determine the changes in natural frequencies of the structural vibration modes that occur at the various temperatures associated with the descent trajectory. This analysis provides a unique and relatively simple design strategy for predicting and mitigating the thermal-structural-vibrational response of hypersonic engines. (C) 2006 Elsevier SAS. All rights reserved.

Solution structure and characterization of the LGR8 receptor binding surface of insulin-like peptide 3

Insulin-like peptide 3 (INSL3), a member of the relaxin peptide family, is produced in testicular Leydig cells and ovarian thecal cells. Gene knock-out experiments have identified a key biological role in initiating testes descent during fetal development. Additionally, INSL3 has an important function in mediating male and female germ cell function. These actions are elicited via its recently identified receptor, LGR8, a member of the leucine-rich repeat-containing G-protein- coupled receptor family. To identify the structural features that are responsible for the interaction of INSL3 with its receptor, its solution structure was determined by NMR spectroscopy together with in vitro assays of a series of B-chain alanine-substituted analogs. Synthetic human INSL3 was found to adopt a characteristic relaxin/ insulin-like fold in solution but is a highly dynamic molecule. The four termini of this two-chain peptide are disordered, and additional conformational exchange is evident in the molecular core. Alanine-substituted analogs were used to identify the key residues of INSL3 that are responsible for the interaction with the ectodomain of LGR8. These include Arg(B16) and Val(B19), with His(B12) and Arg(B20) playing a secondary role, as evident from the synergistic effect on the activity in double and triple mutants involving these residues. Together, these amino acids combine with the previously identified critical residue, Trp(B27), to form the receptor binding surface. The current results provide clear direction for the design of novel specific agonists and antagonists of this receptor.

HLA and Genomewide Allele Sharing in Dizygotic Twins

Gametic selection during fertilization or the effects of specific genotypes on the viability of embryos may cause a skewed transmission of chromosomes to surviving offspring. A recent analysis of transmission distortion in humans reported significant excess sharing among full siblings. Dizygotic (DZ) twin pairs are a special case of the simultaneous survival of two genotypes, and there have been reports of DZ pairs with excess allele sharing around the HLA locus, a candidate locus for embryo survival. We performed an allele-sharing study of 1,592 DZ twin pairs from two independent Australian cohorts, of which 1,561 pairs were informative for linkage on chromosome 6. We also analyzed allele sharing in 336 DZ twin pairs from The Netherlands. We found no evidence of excess allele sharing, either at the HLA locus or in the rest of the genome. In contrast, we found evidence of a small but significant (P = .003 for the Australian sample) genomewide deficit in the proportion of two alleles shared identical by descent among DZ twin pairs. We reconciled conflicting evidence in the literature for excess genomewide allele sharing by performing a simulation study that shows how undetected genotyping errors can lead to an apparent deficit or excess of allele sharing among sibling pairs, dependent on whether parental genotypes are known. Our results imply that gene-mapping studies based on affected sibling pairs that include DZ pairs will not suffer from false-positive results due to loci involved in embryo survival.