Resampling-based approaches to study variation in morphological modularity

Modularity has been suggested to be connected to evolvability because a higher degree of independence among parts allows them to evolve as separate units. Recently, the Escoufier RV coefficient has been proposed as a measure of the degree of integration between modules in multivariate morphometric datasets. However, it has been shown, using randomly simulated datasets, that the value of the RV coefficient depends on sample size. Also, so far there is no statistical test for the difference in the RV coefficient between a priori defined groups of observations. Here, we (1), using a rarefaction analysis, show that the value of the RV coefficient depends on sample size also in real geometric morphometric datasets; (2) propose a permutation procedure to test for the difference in the RV coefficient between a priori defined groups of observations; (3) show, through simulations, that such a permutation procedure has an appropriate Type I error; (4) suggest that a rarefaction procedure could be used to obtain sample-size-corrected values of the RV coefficient; and (5) propose a nearest-neighbor procedure that could be used when studying the variation of modularity in geographic space. The approaches outlined here, readily extendable to non-morphometric datasets, allow study of the variation in the degree of integration between a priori defined modules. A Java application – that will allow performance of the proposed test using a software with graphical user interface – has also been developed and is available at the Morphometrics at Stony Brook Web page (http://life.bio.sunysb.edu/morph/).

Bat species comparisons based on external morphology: A test of traditional versus geometric morphometric approaches

External morphology is commonly used to identify bats as well as to investigate flight and foraging behavior, typically relying on simple length and area measures or ratios. However, geometric morphometrics is increasingly used in the biological sciences to analyse variation in shape and discriminate among species and populations. Here we compare the ability of traditional versus geometric morphometric methods in discriminating between closely related bat species – in this case European horseshoe bats (Rhinolophidae, Chiroptera) – based on morphology of the wing, body and tail. In addition to comparing morphometric methods, we used geometric morphometrics to detect interspecies differences as shape changes. Geometric morphometrics yielded improved species discrimination relative to traditional methods. The predicted shape for the variation along the between group principal components revealed that the largest differences between species lay in the extent to which the wing reaches in the direction of the head. This strong trend in interspecific shape variation is associated with size, which we interpret as an evolutionary allometry pattern.

Characterization of epitopes for virus-neutralizing monoclonal antibodies to Ross River virus E2 using phage-displayed random peptide libraries

Ross River virus (RRV) is the predominant cause of epidemic polyarthritis in Australia, yet the antigenic determinants are not well defined. We aimed to characterize epitope(s) on RRV-E2 for a panel of monoclonal antibodies (MAbs) that recognize overlapping conformational epitopes on the E2 envelope protein of RRV and that neutralize virus infection of cells in vitro. Phage-displayed random peptide libraries were probed with the MAbs T1E7, NB3C4, and T10C9 using solution-phase and solid-phase biopanning methods. The peptides VSIFPPA and KTAISPT were selected 15 and 6 times, respectively, by all three of the MAbs using solution-phase biopanning. The peptide LRLPPAP was selected 8 times by NB3C4 using solid-phase biopanning; this peptide shares a trio of amino acids with the peptide VSIFPPA. Phage that expressed the peptides VSIFPPA and LRLPPAP were reactive with T1E7 and/or NB3C4, and phage that expressed the peptides VSIFPPA, LRLPPAP, and KTAISPT partially inhibited the reactivity of T1E7 with RRV. The selected peptides resemble regions of RRV-E2 adjacent to sites mutated in neutralization escape variants of RRV derived by culture in the presence of these MAbs (E2 210-219 and 238-245) and an additional region of E2 172-182. Together these sites represent a conformational epitope of E2 that is informative of cellular contact sites on RRV.

Interleukin 17A is an immune marker for chlamydial disease severity and pathogenesis in the koala (Phascolarctos cinereus)

The koala (Phascolarctos cinereus) is an iconic Australian marsupial species that is facing many threats to its survival. Chlamydia pecorum infections are a significant contributor to this ongoing decline. A major limiting factor in our ability to manage and control chlamydial disease in koalas is a limited understanding of the koala’s cell-mediated immune response to infections by this bacterial pathogen. To identify immunological markers associated with chlamydial infection and disease in koalas, we used koala-specific Quantitative Real Time PCR (qrtPCR) assays to profile the cytokine responses of Peripheral Blood Mononuclear Cells (PBMCs) collected from 41 koalas with different stages of chlamydial disease. Target cytokines included the principal Th1 (Interferon gamma; IFNγ), Th2 (Interleukin 10; IL10), and pro-inflammatory cytokines (Tumor Necrosis Factor alpha; TNFα). A novel koala-specific IL17A qrtPCR assay was also developed as part of this study to quantitate the gene expression of this Th17 cytokine in koalas. A statistically significant higher IL17A gene expression was observed in animals with current chlamydial disease compared to animals with asymptomatic chlamydial infection. A modest up-regulation of pro-inflammatory cytokines, such as TNFα and IFNγ, was also observed in these animals with signs of current chlamydial disease. IL10 gene expression was not evident in the majority of animals from both groups. Future longitudinal studies are now required to confirm the role played by cytokines in pathology and/or protection against C. pecorum infection in the koala.

Site and gender specificity of inheritance of bone mineral density

Differences in genetic control of BMD by skeletal sites and genders were examined by complex segregation analysis in 816 members of 147 families with probands with extreme low BMD. Spine BMD correlated more strongly in male-male comparisons and hip BMD in female-female comparisons, consistent with gender- and site-specificity of BMD heritability. Introduction: Evidence from studies in animals and humans suggests that the genetic control of bone mineral density (BMD) may differ at different skeletal sites and between genders. This question has important implications for the design and interpretation of genetic studies of osteoporosis. Methods: We examined the genetic profile of 147 families with 816 individuals recruited through probands with extreme low BMD (T-score < −2.5, Z-score < −2.0). Complex segregation analysis was performed using the Pedigree Analysis Package. BMD was measured by DXA at both lumbar spine (L1-L4) and femoral neck. Results: Complex segregation analysis excluded purely monogenic and environmental models of segregation of lumbar spine and femoral neck BMD in these families. Pure polygenic models were excluded at the lumbar spine when menopausal status was considered as a covariate, but not at the femoral neck. Mendelian models with a residual polygenic component were not excluded. These models were consistent with the presence of a rare Mendelian genotype of prevalence 3–19 %, causing high BMD at the hip and spine in these families, with additional polygenic effects. Total heritability range at the lumbar spine was 61–67 % and at the femoral neck was 44–67 %. Significant differences in correlation of femoral neck and lumbar spine BMD were observed between male and female relative pairs, with male-male comparisons exhibiting stronger lumbar spine BMD correlation than femoral neck, and female-female comparisons having greater femoral neck BMD correlation than lumbar spine. These findings remained true for parent-offspring correlations when menopausal status was taken into account. The recurrence risk ratio for siblings of probands of a Z-score < −2.0 was 5.4 at the lumbar spine and 5.9 at the femoral neck. Conclusions: These findings support gender- and site-specificity of the inheritance of BMD. These results should be considered in the design and interpretation of genetic studies of osteoporosis.

Adventures in plant anatomy – not just a pigment of my imagination

An overview of teaching and research activities in the area of plant anatomy at QUT is provided. The current status of teaching of technical skills in plant anatomy is discussed briefly. Examples of applications of plant anatomy to a diverse range of fields are provided, including the crossover between art and science.

Evidence from molecules and morphology expands Podonomopsis Brundin (Diptera: Chironomidae: Podonominae) to include ‘genus Chile'

The informal taxon ‘genus Chile’ of Brundin, based solely on pupal exuviae of a podonomine Chironomidae, has remained inadequately known for half a century. New collections reveal life associations, and provide molecular data to hypothesise a precise phylogenetic placement in the austral Podonominae. A densely sampled molecular phylogeny based on two nuclear and one mitochondrial DNA markers shows ‘genus Chile’ to be the sister group to Podonomopsis Brundin, 1966. Within Podonomopsis a clade of South American species is sister to all Australian species. We discuss how to rank such a sister group taxon and treat ‘genus Chile’ as a new subgenus Araucanopsis, subg. nov. with the new species, Podonomopsis (Araucanopsis) avelasse, sp. nov. from Chile and Argentina as genotype of the monotypic subgenus. We describe P. (A.) avelasse in all stages and provide an expanded diagnosis and description of Podonomopsis to include Araucanopsis. A dated biogeographic hypothesis (chronogram) infers the most recent common ancestor (tmcra) of expanded Podonomopsis at 95 million years ago (Mya) (68–122 Mya 95% highest posterior density), ‘core’ Podonomopsis at 83 Mya (58–108) and Australian Podonomopsis at 65 Mya (44–87). All dates are before the South America–Australia geological separation through Antarctica, supporting previous conclusions that the taxon distribution is ‘Gondwanan’ in origin. Podonomopsis, even as expanded here, remains unknown from New Zealand or elsewhere on extant Zealandia.

A comparative structural bioinformatics analysis of the insulin receptor family ectodomain based on phylogenetic information

The insulin receptor (IR), the insulin-like growth factor 1 receptor (IGF1R) and the insulin receptor-related receptor (IRR) are covalently-linked homodimers made up of several structural domains. The molecular mechanism of ligand binding to the ectodomain of these receptors and the resulting activation of their tyrosine kinase domain is still not well understood. We have carried out an amino acid residue conservation analysis in order to reconstruct the phylogeny of the IR Family. We have confirmed the location of ligand binding site 1 of the IGF1R and IR. Importantly, we have also predicted the likely location of the insulin binding site 2 on the surface of the fibronectin type III domains of the IR. An evolutionary conserved surface on the second leucine-rich domain that may interact with the ligand could not be detected. We suggest a possible mechanical trigger of the activation of the IR that involves a slight ‘twist’ rotation of the last two fibronectin type III domains in order to face the likely location of insulin. Finally, a strong selective pressure was found amongst the IRR orthologous sequences, suggesting that this orphan receptor has a yet unknown physiological role which may be conserved from amphibians to mammals.

'Mutiny on the Bounty': The genetic history of Norfolk Island reveals extreme gender-biased admixture

Background The Pacific Oceania region was one of the last regions of the world to be settled via human migration. Here we outline a settlement of this region that has given rise to a uniquely admixed population. The current Norfolk Island population has arisen from a small number of founders with mixed Caucasian and Polynesian ancestry, descendants of a famous historical event. The ‘Mutiny on the Bounty’ has been told in history books, songs and the big screen, but recently this story can be portrayed through comprehensive molecular genetics. Written history details betrayal and murder leading to the founding of Pitcairn Island by European mutineers and the Polynesian women who left Tahiti with them. Investigation of detailed genealogical records supports historical accounts. Findings Using genetics, we show distinct maternal Polynesian mitochondrial lineages in the present day population, as well as a European centric Y-chromosome phylogeny. These results comprehensively characterise the unique gender-biased admixture of this genetic isolate and further support the historical records relating to Norfolk Island. Conclusions Our results significantly refine previous population genetic studies investigating Polynesian versus Caucasian diversity in the Norfolk Island population and add information that is beneficial to future disease and gene mapping studies.

The complete mitochondrial genome of the eastern grey kangaroo (Macropus giganteus)

We present the complete mitochondrial genome (accession number: LK995454) of an iconic Australian species, the eastern grey kangaroo (Macropus giganteus). The mitogenomic organization is consistent with other marsupials, encoding 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes, an origin of light strand replication and a control region or Dloop. No repetitive sequences were detected in the control region. The M. giganteus mitogenome exemplifies a combination of tRNA gene order and structural peculiarities that appear to be unique to marsupials. We present a maximum likelihood phylogeny based on complete mitochondrial protein and RNA coding sequences that confirms the phylogenetic position of the grey kangaroo among macropodids.

Geomolecular dating and the origin of placental mammals

In modern evolutionary divergence analysis the role of geological information extends beyond providing a timescale, to informing molecular rate variation across the tree. Here I consider the implications of this development. I use fossil calibrations to test the accuracy of models of molecular rate evolution for placental mammals, and reveal substantial misspecification associated with life history rate correlates. Adding further calibrations to reduce dating errors at specific nodes unfortunately tends to transfer underlying rate errors to adjacent branches. Thus, tight calibration across the tree is vital to buffer against rate model errors. I argue that this must include allowing maximum bounds to be tight when good fossil records permit, otherwise divergences deep in the tree will tend to be inflated by the interaction of rate errors and asymmetric confidence in minimum and maximum bounds. In the case of placental mammals I sought to reduce the potential for transferring calibration and rate model errors across the tree by focusing on well-supported calibrations with appropriately conservative maximum bounds. The resulting divergence estimates are younger than others published recently, and provide the long-anticipated molecular signature for the placental mammal radiation observed in the fossil record near the 66 Ma Cretaceous–Paleogene extinction event.

Comment on “Whole-genome analyses resolve early branches in the tree of life of modern birds”

Jarvis et al. (Research Articles, 12 December 2014, p. 1320) presented molecular clock analyses that suggested that most modern bird orders diverged just after the mass extinction event at the Cretaceous-Paleogene boundary (about 66 million years ago). We demonstrate that this conclusion results from the use of a single inappropriate maximum bound, which effectively precludes the Cretaceous diversification overwhelmingly supported by previous molecular studies.

Avian diversification patterns across the K-Pg boundary: Influence of calibrations, datasets, and model misspecification

Birds represent the most diverse extant tetrapod clade, with ca. 10,000 extant species, and the timing of the crown avian radiation remains hotly debated. The fossil record supports a primarily Cenozoic radiation of crown birds, whereas molecular divergence dating analyses generally imply that this radiation was well underway during the Cretaceous. Furthermore, substantial differences have been noted between published divergence estimates. These have been variously attributed to clock model, calibration regime, and gene type. One underappreciated phenomenon is that disparity between fossil ages and molecular dates tends to be proportionally greater for shallower nodes in the avian Tree of Life. Here, we explore potential drivers of disparity in avian divergence dates through a set of analyses applying various calibration strategies and coding methods to a mitochondrial genome dataset and an 18-gene nuclear dataset, both sampled across 72 taxa. Our analyses support the occurrence of two deep divergences (i.e., the Palaeognathae/Neognathae split and the Galloanserae/Neoaves split) well within the Cretaceous, followed by a rapid radiation of Neoaves near the K-Pg boundary. However, 95% highest posterior density intervals for most basal divergences in Neoaves cross the boundary, and we emphasize that, barring unreasonably strict prior distributions, distinguishing between a rapid Early Paleocene radiation and a Late Cretaceous radiation may be beyond the resolving power of currently favored divergence dating methods. In contrast to recent observations for placental mammals, constraining all divergences within Neoaves to occur in the Cenozoic does not result in unreasonably high inferred substitution rates. Comparisons of nuclear DNA (nDNA) versus mitochondrial DNA (mtDNA) datasets and NT- versus RY-coded mitochondrial data reveal patterns of disparity that are consistent with substitution model misspecifications that result in tree compression/tree extension artifacts, which may explain some discordance between previous divergence estimates based on different sequence types. Comparisons of fully calibrated and nominally calibrated trees support a correlation between body mass and apparent dating error. Overall, our results are consistent with (but do not require) a Paleogene radiation for most major clades of crown birds.