A mitochondrial genome phylogeny of termites (Blattodea: Termitoidae) : robust support for interfamilial relationships and molecular synapomorphies define major clades

Despite their ecological significance as decomposers and their evolutionary significance as the most speciose eusocial insect group outside the Hymenoptera, termite (Blattodea: Termitoidae or Isoptera) evolutionary relationships have yet to be well resolved. Previous morphological and molecular analyses strongly conflict at the family level and are marked by poor support for backbone nodes. A mitochondrial (mt) genome phylogeny of termites was produced to test relationships between the recognised termite families, improve nodal support and test the phylogenetic utility of rare genomic changes found in the termite mt genome. Complete mt genomes were sequenced for 7 of the 9 extant termite families with additional representatives of each of the two most speciose families Rhinotermitidae (3 of 7 subfamilies) and Termitidae (3 of 8 subfamilies). The mt genome of the well supported sister group of termites, the subsocial cockroach Cryptocercus, was also sequenced. A highly supported tree of termite relationships was produced by all analytical methods and data treatment approaches, however the relationship of the termites + Cryptocercus clade to other cockroach lineages was highly affected by the strong nucleotide compositional bias found in termites relative to other dictyopterans. The phylogeny supports previously proposed suprafamilial termite lineages, the Euisoptera and Neoisoptera, a later derived Kalotermitidae as sister group of the Neoisoptera and a monophyletic clade of dampwood (Stolotermitidae, Archotermopsidae) and harvester termites (Hodotermitidae). In contrast to previous termite phylogenetic studies, nodal supports were very high for family-level relationships within termites. Two rare genomic changes in the mt genome control region were found to be molecular synapomorphies for major clades. An elongated stem-loop structure defined the clade Polyphagidae + (Cryptocercus + termites), and a further series of compensatory base changes in this stem loop is synapomorphic for the Neoisoptera. The complicated repeat structures first identified in Reticulitermes, composed of short (A-type) and long (B-type repeats) defines the clade Heterotermitinae + Termitidae, while the secondary loss of A-type repeats is synapomorphic for the non-macrotermitine Termitidae.

Molecular characterization of "Candidatus Similichlamydia latridicola" gen. nov., sp. nov. (Chlamydiales: "Candidatus Parilichlamydiaceae"), a novel Chlamydia-like epitheliocystis agent in the striped trumpeter, Latris lineata (Forster)

Histological analysis of gill samples taken from individuals of Latris lineata reared in aquaculture in Tasmania, Australia, and those sampled from the wild revealed the presence of epitheliocystis-like basophilic inclusions. Subsequent morphological, in situ hybridization, and molecular analyses were performed to confirm the presence of this disease and discovered a Chlamydia-like organism associated with this condition, and the criteria set by Fredericks and Relman's postulates were used to establish disease causation. Three distinct 16S rRNA genotypes were sequenced from 16 fish, and phylogenetic analyses of the nearly full-length 16S rRNA sequences generated for this bacterial agent indicated that they were nearly identical novel members of the order Chlamydiales. This new taxon formed a well-supported clade with "Candidatus Parilichlamydia carangidicola" from the yellowtail kingfish (Seriola lalandi). On the basis of sequence divergence over the 16S rRNA region relative to all other members of the order Chlamydiales, a new genus and species are proposed here for the Chlamydia-like bacterium from L. lineata, i.e., "Candidatus Similichlamydia latridicola" gen. nov., sp. nov.

One and the same: Integrative taxonomic evidence that Bactrocera invadens (Diptera: Tephritidae) is the same species as the Oriental fruit fly Bactrocera dorsalis

The invasive fruit fly Bactrocera invadens Drew, Tsuruta & White, and the Oriental fruit fly Bactrocera dorsalis (Hendel) are highly destructive horticultural pests of global significance. Bactrocera invadens originates from the Indian subcontinent and has recently invaded all of sub-Saharan Africa, while B. dorsalis principally occurs from the Indian subcontinent towards southern China and South-east Asia. High morphological and genetic similarity has cast doubt over whether B. invadens is a distinct species from B. dorsalis. Addressing this issue within an integrative taxonomic framework, we sampled from across the geographic distribution of both taxa and: (i) analysed morphological variation, including those characters considered diagnostic (scutum colour, length of aedeagus, width of postsutural lateral vittae, wing size, and wing shape); (ii) sequenced four loci (ITS1, ITS2, cox1 and nad4) for phylogenetic inference, and; (iii) generated a cox1 haplotype network to examine population structure. Molecular analyses included the closely related species, Bactrocera kandiensis Drew & Hancock. Scutum colour varies from red-brown to fully black for individuals from Africa and the Indian subcontinent. All individuals east of the Indian subcontinent are black except for a few red-brown individuals from China. The postsutural lateral vittae width of B. invadens is narrower than B. dorsalis from eastern Asia, but the variation is clinal, with subcontinent B. dorsalis populations intermediate in size. Aedeagus length, wing shape and wing size cannot discriminate between the two taxa. Phylogenetic analyses failed to resolve B. invadens from B. dorsalis, but did resolve B. kandiensis. Bactrocera dorsalis and B. invadens shared cox1 haplotypes, yet the haplotype network pattern does not reflect current taxonomy or patterns in thoracic colour. Some individuals of B. dorsalis/B. invadens possessed haplotypes more closely related to B. kandiensis than to conspecifics, suggestive of mitochondrial introgression between these species. The combined evidence fails to support the delimitation of B. dorsalis and B. invadens as separate biological species. Consequently, existing biological data for B. dorsalis may be applied to the invasive population in Africa. Our recommendation, in line with other recent publications, is that B. invadens be synonymized with B. dorsalis.

Refining prognosis and identifying targetable pathways for high-risk endometrial cancer; a TransPORTEC initiative

This study aimed to investigate whether molecular analysis can be used to refine risk assessment, direct adjuvant therapy, and identify actionable alterations in high-risk endometrial cancer. TransPORTEC, an international consortium related to the PORTEC3 trial, was established for translational research in high-risk endometrial cancer. In this explorative study, routine molecular analyses were used to detect prognostic subgroups: p53 immunohistochemistry, microsatellite instability and POLE proofreading mutation. Furthermore, DNA was analyzed for hotspot mutations in 13 additional genes (BRAF, CDKNA2, CTNNB1, FBXW7, FGFR2, FGFR3, FOXL2, HRAS, KRAS, NRAS, PIK3CA, PPP2R1A, and PTEN) and protein expression of ER, PR, PTEN, and ARID1a was analyzed. Rates of distant metastasis, recurrence-free, and overall survival were calculated using the Kaplan-Meier method and log-rank test. In total, samples of 116 high-risk endometrial cancer patients were included: 86 endometrioid; 12 serous; and 18 clear cell. For endometrioid, serous, and clear cell cancers, 5-year recurrence-free survival rates were 68%, 27%, and 50% (P=0.014) and distant metastasis rates 23%, 64%, and 50% (P=0.001), respectively. Four prognostic subgroups were identified: (1) a group of p53-mutant tumors; (2) microsatellite instable tumors; (3) POLE proofreading-mutant tumors; and (4) a group with no specific molecular profile (NSMP). In group 3 (POLE-mutant; n=14) and group 2 (microsatellite instable; n=19) patients, no distant metastasis occurred, compared with 50% distant metastasis rate in group 1 (p53-mutant; n=36) and 39% in group 4 (NSMP; P<0.001). Five-year recurrence-free survival was 93% and 95% for group 3 (POLE-mutant) and group 2 (microsatellite instable) vs 42% (group 1, p53-mutant) and 52% (group 4, NSMP; P<0.001). Targetable FBXW7 and FGFR2 mutations (6%), alterations in the PI3K-AKT pathway (60%) and hormone receptor positivity (45%) were frequently found. In conclusion, molecular analysis of high-risk endometrial cancer identifies four distinct prognostic subgroups, with potential therapeutic implications. High frequencies of targetable alterations were identified and may serve as targets for individualized treatment

Resolving cryptic species complexes of major tephritid pests

An FAO/IAEA Co-ordinated Research Project (CRP) on “Resolution of Cryptic Species Complexes of Tephritid Pests to Overcome Constraints to SIT Application and International Trade” was conducted from 2010 to 2015. As captured in the CRP title, the objective was to undertake targeted research into the systematics and diagnostics of taxonomically challenging fruit fly groups of economic importance. The scientific output was the accurate alignment of biological species with taxonomic names; which led to the applied outcome of assisting FAO and IAEA Member States in overcoming technical constraints to the application of the Sterile Insect Technique (SIT) against pest fruit flies and the facilitation of international agricultural trade. Close to 50 researchers from over 20 countries participated in the CRP, using coordinated, multidisciplinary research to address, within an integrative taxonomic framework, cryptic species complexes of major tephritid pests. The following progress was made for the four complexes selected and studied: Anastrepha fraterculus complex – Eight morphotypes and their geographic and ecological distributions in Latin America were defined. The morphotypes can be considered as distinct biological species on the basis of differences in karyotype, sexual incompatibility, post-mating isolation, cuticular hydrocarbon, pheromone, and molecular analyses. Discriminative taxonomic tools using linear and geometric morphometrics of both adult and larval morphology were developed for this complex. Bactrocera dorsalis complex – Based on genetic, cytogenetic, pheromonal, morphometric, and behavioural data, which showed no or only minor variation between the Asian/African pest fruit flies Bactrocera dorsalis, B. papayae, B. philippinensis and B. invadens, the latter three species were synonymized with B. dorsalis. Of the five target pest taxa studied, only B. dorsalis and B. carambolae remain as scientifically valid names. Molecular and pheromone markers are now available to distinguish B. dorsalis from B. carambolae. Ceratitis FAR Complex (C. fasciventris, C. anonae, C. rosa) – Morphology, morphometry, genetic, genomic, pheromone, cuticular hydrocarbon, ecology, behaviour, and developmental physiology data provide evidence for the existence of five different entities within this fruit fly complex from the African region. These are currently recognised as Ceratitis anonae, C. fasciventris (F1 and F2), C. rosa and a new species related to C. rosa (R2). The biological limits within C. fasciventris (i.e. F1 and F2) are not fully resolved. Microsatellites markers and morphological identification tools for the adult males of the five different FAR entities were developed based on male leg structures. Zeugodacus cucurbitae (formerly Bactrocera (Zeugodacus) cucurbitae) – Genetic variability was studied among melon fly populations throughout its geographic range in Africa and the Asia/Pacific region and found to be limited. Cross-mating studies indicated no incompatibility or sexual isolation. Host preference and genetic studies showed no evidence for the existence of host races. It was concluded that the melon fly does not represent a cryptic species complex, neither with regard to geographic distribution nor to host range. Nevertheless, the higher taxonomic classification under which this species had been placed, by the time the CRP was started, was found to be paraphyletic; as a result the subgenus Zeugodacus was elevated to genus level.

Proteomic identification of putative microRNA394 target genes in Arabidopsis thaliana identifies MLP family members critical for normal development

Expression of the F-Box protein Leaf Curling Responsiveness (LCR) is regulated by microRNA, miR394, and alterations to this interplay in Arabidopsis thaliana produce defects in leaf polarity and shoot apical meristem (SAM) organisation. Although the miR394-LCR node has been documented in Arabidopsis, the identification of proteins targeted by LCR F-box itself has proven problematic. Here, a proteomic analysis of shoot apices from plants with altered LCR levels identified a member of the Major Latex Protein (MLP) family gene as a potential LCR F-box target. Bioinformatic and molecular analyses also suggested that other MLP family members are likely to be targets for this post-translational regulation. Direct interaction between LCR F-Box and MLP423 was validated. Additional MLP members had reduction in protein accumulation, in varying degrees, mediated by LCR F-Box. Transgenic Arabidopsis lines, in which MLP28 expression was reduced through an artificial miRNA technology, displayed severe developmental defects, including changes in leaf patterning and morphology, shoot apex defects, and eventual premature death. These phenotypic characteristics resemble those of Arabidopsis plants modified to over-express LCR. Taken together, the results demonstrate that MLPs are driven to degradation by LCR, and indicate that MLP gene family is target of miR394-LCR regulatory node, representing potential targets for directly post-translational regulation mediated by LCR F-Box. In addition, MLP28 family member is associated with the LCR regulation that is critical for normal Arabidopsis development.

Molecular investigation of the mechanical properties of single actin filaments based on vibration analyses

The mechanical vibration properties of single actin filaments from 50 to 288 nm are investigated by the molecular dynamics simulation in this study. The natural frequencies obtained from the molecular simulations agree with those obtained from the analytical solution of the equivalent Euler–Bernoulli beam model. Through the convergence study of the mechanical properties with respect to the filament length, it was found that the Euler–Bernoulli beam model can only be reliably used when the single actin filament is of the order of hundreds of nanometre scale. This molecular investigation not only provides the evidence for the use of the continuum beam model in characterising the mechanical properties of single actin filaments, but also clarifies the criteria for the effective use of the Euler–Bernoulli beam model.

Molecular genetic analyses of RFLPs for PCR-amplified growth hormone gene, renal kallikrein gene and atrial natriuretic factor gene in essential hypertension

The present study examined polymorphisms of genes that might be involved in the onset of essential hypertension (HT). These included the (i) growth hormone gene (GH1), whose locus has recently been linked to elevated blood pressure (BP) in the stroke-prone SHR, although recent sib-pair analysis of a polymorphism near the human chorionic somatomammotropin gene (a member of the GH cluster) was unable to show linkage with HT; (ii) renal kallikrein gene (KLK1); and (iii) atrial natriuretic factor gene (ANF), where a primary defect in production or activity of kallikrein or ANF could cause NaCl retention and vasoconstriction. Association analyses were conducted to compare restriction fragment length polymorphisms (RFLPs) of each gene in 85 HT and 95 normotensive (NT) Caucasian subjects whose parents had a similar BP status at age ≥50 years. The frequency of the minor allele of (i) a RsaI RFLP in the promoter of GH1, amplified from leukocyte DNA by the polymerase chain reaction, was 0.15 in the HT group and 0.14 in the NT group (χ1=0.34, P=0.55); (ii) a TaqI RFLP for KLK1 was 0.035 in the HT group and 0.015 in the NT group (χ2=1.5, P=0.21); and (iii) a XhoI RFLP for ANF was 0.50 in HTs and 0.46 in NTs (χ2=0.20, P=0.65). Studies of HT pedigrees found one family in which the ANF locus and HT were not linked, owing to an obligate recombinant. The present data thus provide no evidence for involvement of the growth hormone, renal kallikrein, nor ANF gene in the causation of essential hypertension.

Fractal and complex network analyses of protein molecular dynamics

Based on protein molecular dynamics, we investigate the fractal properties of energy, pressure and volume time series using the multifractal detrended fluctuation analysis (MF-DFA) and the topological and fractal properties of their converted horizontal visibility graphs (HVGs). The energy parameters of protein dynamics we considered are bonded potential, angle potential, dihedral potential, improper potential, kinetic energy, Van der Waals potential, electrostatic potential, total energy and potential energy. The shape of the h(q)h(q) curves from MF-DFA indicates that these time series are multifractal. The numerical values of the exponent h(2)h(2) of MF-DFA show that the series of total energy and potential energy are non-stationary and anti-persistent; the other time series are stationary and persistent apart from series of pressure (with H≈0.5H≈0.5 indicating the absence of long-range correlation). The degree distributions of their converted HVGs show that these networks are exponential. The results of fractal analysis show that fractality exists in these converted HVGs. For each energy, pressure or volume parameter, it is found that the values of h(2)h(2) of MF-DFA on the time series, exponent λλ of the exponential degree distribution and fractal dimension dBdB of their converted HVGs do not change much for different proteins (indicating some universality). We also found that after taking average over all proteins, there is a linear relationship between 〈h(2)〉〈h(2)〉 (from MF-DFA on time series) and 〈dB〉〈dB〉 of the converted HVGs for different energy, pressure and volume.

Exclusion of angiotensinogen gene in molecular basis of human hypertension: Sibpair linkage and association analyses in Australian anglo-caucasians

Linkage with essential hypertension has been claimed for a microsatellite marker near the angiotensinogen gene (AGT; chromosome 1q42), as has association for the AGT variants M235T, G(-6)A and A(-20)C. To more rigorously evaluate AGT as a candidate gene for hypertension we performed sibpair analysis with multiple microsatellite markers surrounding this locus and using more sophisticated analysis programs. We also performed an association study of the AGT variants in unrelated subjects with a strong family history (two affected parents). For the linkage study, single and multiplex polymerase chain reaction (PCRs) and automated genescan analysis were conducted on DNA from 175 Australian Anglo-Celtic Caucasian hypertensives for the following markers: D1S2880-(2.1 cM)-D1S213-(2.8 cM)-D1S251-(6.5 cM)-AGT-(2.0 cM) -D1S235. Statistical evaluation of genotype data by nonparametric methods resulted in the following scores: Single-point analysis - SPLINK, P > 0.18; APM method, P > 0.25; ASPEX, MLOD < 0.28; SIB-PAIR, P > 0. 24; Multipoint analysis - MAPMAKER/SIBS, MLOD < 0.24; GENEHUNTER, P > 0.35. Exclusion scores of Lod -4.1 to -5.1 were obtained for these markers using MAPMAKER/SIBS for a lambda(s) of 1.6. The association study of G(-6)A, A(-20)C and M235T variants in 111 hypertensives with strong family history and 190 normotensives with no family history showed significant linkage disequilibrium between particular haplotypes, but we could find no association with hypertension. The present study therefore excludes AGT in the etiology of hypertension, at least in the population of Australian Anglo-Celtic Caucasians studied.

Molecular analysis of the myosin gene family in Arabidopsis thaliana

Myosin is believed to act as the molecular motor for many actin-based motility processes in eukaryotes. It is becoming apparent that a single species may possess multiple myosin isoforms, and at least seven distinct classes of myosin have been identified from studies of animals, fungi, and protozoans. The complexity of the myosin heavy-chain gene family in higher plants was investigated by isolating and characterizing myosin genomic and cDNA clones from Arabidopsis thaliana. Six myosin-like genes were identified from three polymerase chain reaction (PCR) products (PCR1, PCR11, PCR43) and three cDNA clones (ATM2, MYA2, MYA3). Sequence comparisons of the deduced head domains suggest that these myosins are members of two major classes. Analysis of the overall structure of the ATM2 and MYA2 myosins shows that they are similar to the previously-identified ATM1 and MYA1 myosins, respectively. The MYA3 appears to possess a novel tail domain, with five IQ repeats, a six-member imperfect repeat, and a segment of unique sequence. Northern blot analyses indicate that some of the Arabidopsis myosin genes are preferentially expressed in different plant organs. Combined with previous studies, these results show that the Arabidopsis genome contains at least eight myosin-like genes representing two distinct classes.

Genetic and gene expression analyses of the polycystic ovary syndrome candidate gene fibrillin-3 and other fibrillin family members in human ovaries

Several studies have demonstrated an association between polycystic ovary syndrome (PCOS) and the dinucleotide repeat microsatellite marker D19S884, which is located in intron 55 of the fibrillin-3 (FBN3) gene. Fibrillins, including FBN1 and 2, interact with latent transforming growth factor (TGF)-β-binding proteins (LTBP) and thereby control the bioactivity of TGFβs. TGFβs stimulate fibroblast replication and collagen production. The PCOS ovarian phenotype includes increased stromal collagen and expansion of the ovarian cortex, features feasibly influenced by abnormal fibrillin expression. To examine a possible role of fibrillins in PCOS, particularly FBN3, we undertook tagging and functional single nucleotide polymorphism (SNP) analysis (32 SNPs including 10 that generate non-synonymous amino acid changes) using DNA from 173 PCOS patients and 194 controls. No SNP showed a significant association with PCOS and alleles of most SNPs showed almost identical population frequencies between PCOS and control subjects. No significant differences were observed for microsatellite D19S884. In human PCO stroma/cortex (n = 4) and non-PCO ovarian stroma (n = 9), follicles (n = 3) and corpora lutea (n = 3) and in human ovarian cancer cell lines (KGN, SKOV-3, OVCAR-3, OVCAR-5), FBN1 mRNA levels were approximately 100 times greater than FBN2 and 200–1000-fold greater than FBN3. Expression of LTBP-1 mRNA was 3-fold greater than LTBP-2. We conclude that FBN3 appears to have little involvement in PCOS but cannot rule out that other markers in the region of chromosome 19p13.2 are associated with PCOS or that FBN3 expression occurs in other organs and that this may be influencing the PCOS phenotype.

Time dependency of molecular rate estimates and systematic overestimation of recent divergence times

Studies of molecular evolutionary rates have yielded a wide range of rate estimates for various genes and taxa. Recent studies based on population-level and pedigree data have produced remarkably high estimates of mutation rate, which strongly contrast with substitution rates inferred in phylogenetic (species-level) studies. Using Bayesian analysis with a relaxed-clock model, we estimated rates for three groups of mitochondrial data: avian protein-coding genes, primate protein-coding genes, and primate d-loop sequences. In all three cases, we found a measurable transition between the high, short-term (<1–2 Myr) mutation rate and the low, long-term substitution rate. The relationship between the age of the calibration and the rate of change can be described by a vertically translated exponential decay curve, which may be used for correcting molecular date estimates. The phylogenetic substitution rates in mitochondria are approximately 0.5% per million years for avian protein-coding sequences and 1.5% per million years for primate protein-coding and d-loop sequences. Further analyses showed that purifying selection offers the most convincing explanation for the observed relationship between the estimated rate and the depth of the calibration. We rule out the possibility that it is a spurious result arising from sequence errors, and find it unlikely that the apparent decline in rates over time is caused by mutational saturation. Using a rate curve estimated from the d-loop data, several dates for last common ancestors were calculated: modern humans and Neandertals (354 ka; 222–705 ka), Neandertals (108 ka; 70–156 ka), and modern humans (76 ka; 47–110 ka). If the rate curve for a particular taxonomic group can be accurately estimated, it can be a useful tool for correcting divergence date estimates by taking the rate decay into account. Our results show that it is invalid to extrapolate molecular rates of change across different evolutionary timescales, which has important consequences for studies of populations, domestication, conservation genetics, and human evolution.

Branch-length estimation bias misleads molecular dating for a vertebrate mitochondrial phylogeny

Despite recent methodological advances in inferring the time-scale of biological evolution from molecular data, the fundamental question of whether our substitution models are sufficiently well specified to accurately estimate branch-lengths has received little attention. I examine this implicit assumption of all molecular dating methods, on a vertebrate mitochondrial protein-coding dataset. Comparison with analyses in which the data are RY-coded (AG → R; CT → Y) suggests that even rates-across-sites maximum likelihood greatly under-compensates for multiple substitutions among the standard (ACGT) NT-coded data, which has been subject to greater phylogenetic signal erosion. Accordingly, the fossil record indicates that branch-lengths inferred from the NT-coded data translate into divergence time overestimates when calibrated from deeper in the tree. Intriguingly, RY-coding led to the opposite result. The underlying NT and RY substitution model misspecifications likely relate respectively to “hidden” rate heterogeneity and changes in substitution processes across the tree, for which I provide simulated examples. Given the magnitude of the inferred molecular dating errors, branch-length estimation biases may partly explain current conflicts with some palaeontological dating estimates.