Complexity of avian evolution revealed by family-level genomes.

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

The water lily genome and the early evolution of flowering plants.

Water lilies belong to the angiosperm order Nymphaeales. Amborellales, Nymphaeales and Austrobaileyales together form the so-called ANA-grade of angiosperms, which are extant representatives of lineages that diverged the earliest from the lineage leading to the extant mesangiosperms1-3. Here we report the 409-megabase genome sequence of the blue-petal water lily (Nymphaea colorata). Our phylogenomic analyses support Amborellales and Nymphaeales as successive sister lineages to all other extant angiosperms. The N. colorata genome and 19 other water lily transcriptomes reveal a Nymphaealean whole-genome duplication event, which is shared by Nymphaeaceae and possibly Cabombaceae. Among the genes retained from this whole-genome duplication are homologues of genes that regulate flowering transition and flower development. The broad expression of homologues of floral ABCE genes in N. colorata might support a similarly broadly active ancestral ABCE model of floral organ determination in early angiosperms. Water lilies have evolved attractive floral scents and colours, which are features shared with mesangiosperms, and we identified their putative biosynthetic genes in N. colorata. The chemical compounds and biosynthetic genes behind floral scents suggest that they have evolved in parallel to those in mesangiosperms. Because of its unique phylogenetic position, the N. colorata genome sheds light on the early evolution of angiosperms.

Dating in the Dark: Elevated Substitution Rates in Cave Cockroaches (Blattodea: Nocticolidae) Have Negative Impacts on Molecular Date Estimates.

Rates of nucleotide substitution vary substantially across the Tree of Life, with potentially confounding effects on phylogenetic and evolutionary analyses. A large acceleration in mitochondrial substitution rate occurs in the cockroach family Nocticolidae, which predominantly inhabit subterranean environments. To evaluate the impacts of this among-lineage rate heterogeneity on estimates of phylogenetic relationships and evolutionary timescales, we analyzed nuclear ultraconserved elements (UCEs) and mitochondrial genomes from nocticolids and other cockroaches. Substitution rates were substantially elevated in nocticolid lineages compared with other cockroaches, especially in mitochondrial protein-coding genes. This disparity in evolutionary rates is likely to have led to different evolutionary relationships being supported by phylogenetic analyses of mitochondrial genomes and UCE loci. Furthermore, Bayesian dating analyses using relaxed-clock models inferred much deeper divergence times compared with a flexible local clock. Our phylogenetic analysis of UCEs, which is the first genome-scale study to include all 13 major cockroach families, unites Corydiidae and Nocticolidae and places Anaplectidae as the sister lineage to the rest of Blattoidea. We uncover an extraordinary level of genetic divergence in Nocticolidae, including two highly distinct clades that separated ~115 million years ago despite both containing representatives of the genus Nocticola. The results of our study highlight the potential impacts of high among-lineage rate variation on estimates of phylogenetic relationships and evolutionary timescales.

Pharmacological Improvement of Cystic Fibrosis Transmembrane Conductance Regulator Function Rescues Airway Epithelial Homeostasis and Host Defense in Children with Cystic Fibrosis.

Rationale: Pharmacological improvement of cystic fibrosis transmembrane conductance regulator (CFTR) function with elexacaftor/tezacaftor/ivacaftor (ETI) provides unprecedented improvements in lung function and other clinical outcomes in patients with cystic fibrosis (CF). However, ETI effects on impaired mucosal homeostasis and host defense at the molecular and cellular levels in the airways of patients with CF remain unknown. Objectives: To investigate effects of ETI on the transcriptome of nasal epithelial and immune cells from children with CF at the single-cell level. Methods: Nasal swabs from 13 children with CF and at least one F508del allele aged 6 to 11 years were collected at baseline and 3 months after initiation of ETI, subjected to single-cell RNA sequencing, and compared with swabs from 12 age-matched healthy children. Measurements and Main Results: Proportions of CFTR-positive cells were decreased in epithelial basal, club, and goblet cells, but not in ionocytes, from children with CF at baseline and were restored by ETI therapy to nearly healthy levels. Single-cell transcriptomics revealed an impaired IFN signaling and reduced expression of major histocompatibility complex classes I and II encoding genes in epithelial cells of children with CF at baseline, which was partially restored by ETI. In addition, ETI therapy markedly reduced the inflammatory phenotype of immune cells, particularly of neutrophils and macrophages. Conclusions: Pharmacological improvement of CFTR function improves innate mucosal immunity and reduces immune cell inflammatory responses in the upper airways of children with CF at the single-cell level, highlighting the potential to restore epithelial homeostasis and host defense in CF airways by early initiation of ETI therapy.

The future of cystic fibrosis treatment: from disease mechanisms to novel therapeutic approaches.

With the 2019 breakthrough in the development of highly effective modulator therapy providing unprecedented clinical benefits for over 90% of patients with cystic fibrosis who are genetically eligible for treatment, this rare disease has become a front runner of transformative molecular therapy. This success is based on fundamental research, which led to the identification of the disease-causing CFTR gene and our subsequent understanding of the disease mechanisms underlying the pathogenesis of cystic fibrosis, working together with a continuously evolving clinical research and drug development pipeline. In this Series paper, we focus on advances since 2018, and remaining knowledge gaps in our understanding of the molecular mechanisms of CFTR dysfunction in the airway epithelium and their links to mucus dysfunction, impaired host defences, airway infection, and chronic inflammation of the lungs of people with cystic fibrosis. We review progress in (and the remaining obstacles to) pharmacological approaches to rescue CFTR function, and novel strategies for improved symptomatic therapies for cystic fibrosis, including how these might be applicable to common lung diseases, such as bronchiectasis and chronic obstructive pulmonary disease. Finally, we discuss the promise of genetic therapies and gene editing approaches to restore CFTR function in the lungs of all patients with cystic fibrosis independent of their CFTR genotype, and the unprecedented opportunities to transform cystic fibrosis from a fatal disease to a treatable and potentially curable one.

Impact of elexacaftor/tezacaftor/ivacaftor therapy on lung clearance index and magnetic resonance imaging in children with cystic fibrosis and one or two F508del alleles.

BACKGROUND: We recently demonstrated that elexacaftor/tezacaftor/ivacaftor (ETI) improves the lung clearance index (LCI) and abnormalities in lung morphology detected by magnetic resonance imaging (MRI) in adolescent and adult patients with cystic fibrosis (CF). However, real-world data on the effect of ETI on these sensitive outcomes of lung structure and function in school-age children with CF have not been reported. The aim of this study was therefore to examine the effect of ETI on the LCI and the lung MRI score in children aged 6-11 years with CF and one or two F508del alleles. METHODS: This prospective, observational, multicentre, post-approval study assessed the longitudinal LCI up to 12 months and the lung MRI score before and 3 months after initiation of ETI. RESULTS: A total of 107 children with CF including 40 heterozygous for F508del and a minimal function mutation (F/MF) and 67 homozygous for F508del (F/F) were enrolled in this study. Treatment with ETI improved the median (interquartile range (IQR)) LCI in F/MF (-1.0 (-2.0- -0.1); p<0.01) and F/F children (-0.8 (-1.9- -0.2); p<0.001) from 3 months onwards. Further, ETI improved the median (IQR) MRI global score in F/MF (-4.0 (-9.0-0.0); p<0.01) and F/F children (-3.5 (-7.3- -0.8); p<0.001). CONCLUSIONS: ETI improves early abnormalities in lung ventilation and morphology in school-age children with CF and at least one F508del allele in a real-world setting. Our results support early initiation of ETI to reduce or even prevent lung disease progression in school-age children with CF.

Pervasive relaxed selection in termite genomes.

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Modulation of TRPV4 protects against degeneration induced by sustained loading and promotes matrix synthesis in the intervertebral disc.

While it is well known that mechanical signals can either promote or disrupt intervertebral disc (IVD) homeostasis, the molecular mechanisms for transducing mechanical stimuli are not fully understood. The transient receptor potential vanilloid 4 (TRPV4) ion channel activated in isolated IVD cells initiates extracellular matrix (ECM) gene expression, while TRPV4 ablation reduces cytokine production in response to circumferential stretching. However, the role of TRPV4 on ECM maintenance during tissue-level mechanical loading remains unknown. Using an organ culture model, we modulated TRPV4 function over both short- (hours) and long-term (days) and evaluated the IVDs' response. Activating TRPV4 with the agonist GSK101 resulted in a Ca2+ flux propagating across the cells within the IVD. Nuclear factor (NF)-κB signaling in the IVD peaked at 6 h following TRPV4 activation that subsequently resulted in higher interleukin (IL)-6 production at 7 days. These cellular responses were concomitant with the accumulation of glycosaminoglycans and increased hydration in the nucleus pulposus that culminated in higher stiffness of the IVD. Sustained compressive loading of the IVD resulted in elevated NF-κB activity, IL-6 and vascular endothelial growth factor A (VEGFA) production, and degenerative changes to the ECM. TRPV4 inhibition using GSK205 during loading mitigated the changes in inflammatory cytokines, protected against IVD degeneration, but could not prevent ECM disorganization due to mechanical damage in the annulus fibrosus. These results indicate TRPV4 plays an important role in both short- and long-term adaptations of the IVD to mechanical loading. The modulation of TRPV4 may be a possible therapeutic for preventing load-induced IVD degeneration.

Evaluating the Accuracy of Methods for Detecting Correlated Rates of Molecular and Morphological Evolution.

Determining the link between genomic and phenotypic change is a fundamental goal in evolutionary biology. Insights into this link can be gained by using a phylogenetic approach to test for correlations between rates of molecular and morphological evolution. However, there has been persistent uncertainty about the relationship between these rates, partly because conflicting results have been obtained using various methods that have not been examined in detail. We carried out a simulation study to evaluate the performance of 5 statistical methods for detecting correlated rates of evolution. Our simulations explored the evolution of molecular sequences and morphological characters under a range of conditions. Of the methods tested, Bayesian relaxed-clock estimation of branch rates was able to detect correlated rates of evolution correctly in the largest number of cases. This was followed by correlations of root-to-tip distances, Bayesian model selection, independent sister-pairs contrasts, and likelihood-based model selection. As expected, the power to detect correlated rates increased with the amount of data, both in terms of tree size and number of morphological characters. Likewise, greater among-lineage rate variation in the data led to improved performance of all 5 methods, particularly for Bayesian relaxed-clock analysis when the rate model was mismatched. We then applied these methods to a data set from flowering plants and did not find evidence of a correlation in evolutionary rates between genomic data and morphological characters. The results of our study have practical implications for phylogenetic analyses of combined molecular and morphological data sets, and highlight the conditions under which the links between genomic and phenotypic rates of evolution can be evaluated quantitatively.

Impacts of Taxon-Sampling Schemes on Bayesian Tip Dating Under the Fossilized Birth-Death Process.

Evolutionary timescales can be inferred by molecular-clock analyses of genetic data and fossil evidence. Bayesian phylogenetic methods such as tip dating provide a powerful framework for inferring evolutionary timescales, but the most widely used priors for tree topologies and node times often assume that present-day taxa have been sampled randomly or exhaustively. In practice, taxon sampling is often carried out so as to include representatives of major lineages, such as orders or families. We examined the impacts of different densities of diversified sampling on Bayesian tip dating on unresolved fossilized birth-death (FBD) trees, in which fossil taxa are topologically constrained but their exact placements are averaged out. We used synthetic data generated by simulations of nucleotide sequence evolution, fossil occurrences, and diversified taxon sampling. Our analyses under the diversified-sampling FBD process show that increasing taxon-sampling density does not necessarily improve divergence-time estimates. However, when informative priors were specified for the root age or when tree topologies were fixed to those used for simulation, the performance of tip dating on unresolved FBD trees maintains its accuracy and precision or improves with taxon-sampling density. By exploring three situations in which models are mismatched, we find that including all relevant fossils, without pruning off those that are incompatible with the diversified-sampling FBD process, can lead to underestimation of divergence times. Our reanalysis of a eutherian mammal data set confirms some of the findings from our simulation study, and reveals the complexity of diversified taxon sampling in phylogenomic data sets. In highlighting the interplay of taxon-sampling density and other factors, the results of our study have practical implications for using Bayesian tip dating to infer evolutionary timescales across the Tree of Life. [Bayesian tip dating; eutherian mammals; fossilized birth-death process; phylogenomics; taxon sampling.].

Treatment and Long-Term Safety Outcomes of Peptide Receptor Radionuclide Therapy for Metastatic Neuroendocrine Tumours: An Asian Experience.

PURPOSE: Peptide receptor radionuclide therapy (PRRT) is a targeted molecular therapy used to treat neuroendocrine tumours (NETs). It has been shown to be effective and well tolerated in patients with metastatic NETs in several centres in the USA, Europe, and Australia. Tolerability and efficacy data emerging from Asian centres remain few. Epidemiological evidence suggests that there are differences in neuroendocrine neoplasms between the population groups. We aim to describe the treatment and safety outcomes of PRRT in the Asian population. METHODS: One hundred and seven (107) patients with metastatic NETs who had undergone PRRT treatment from January 2012 to March 2019 were included in this retrospective study. The response rates using RECIST 1.1 and qualitative analysis were examined. The overall and progression-free survival curves were also evaluated. RESULTS: The median progression-free survival was 49 months. Response assessment after completion of treatment showed that 33 (37.9%) of 87 patients had partial or complete response. Subgroup analysis comparing high- and low-grade NET showed that there was a significant difference in the time to progression curves. Comparison of the number of cycles and progression-free and overall survival also showed a significant difference. Ten patients (9%) had grade 3 or more haematological toxicities. Four patients (4%) had grade 3/4 hepatobiliary toxicities, although the presence of extensive liver metastases was a confounding factor. None of the patients had grade 3/4 acute kidney injury. CONCLUSION: Our results show that PRRT is safe and effective in the treatment of metastatic NET in the Asian population. There was a significant difference in the progression-free survival curves between low-grade and high-grade NET and in the progression-free and overall survival comparing the number of cycles received.

Outcomes after per-oral endoscopic myotomy for Zenker's diverticula (Z-POEM) and correlation with impedance planimetry (FLIP).

INTRODUCTION: Zenker's diverticulum (ZD) is a false pulsion diverticulum of the cervical esophagus. It is typically found in older adults and manifests with dysphagia. The purpose of this study is to describe our experience with Per-oral endoscopic myotomy for Zenker's (Z-POEM) and intraoperative impedance planimetry (FLIP). METHODS: We performed a single institution retrospective review of patients undergoing Z-POEM in a prospective database between 2014 and 2022. Upper esophageal sphincter (UES) distensibility index (DI, mm2/mmHg) was measured by FLIP before and after myotomy. The primary outcome was clinical success. Secondary outcomes included technical failure, adverse events, and quality of life as assessed by the gastroesophageal health-related quality of life (GERD-HRQL), reflux severity index (RSI), and dysphagia score. A statistical analysis of DI was done with the paired t-test (p < 0.05). RESULTS: Fifty-four patients underwent Z-POEM, with FLIP measurements available in 30 cases. We achieved technical success and clinical success in 54/54 (100%) patients and 46/54 patients (85%), respectively. Three patients (6%) experienced contained leaks. Three patients were readmitted: one for aforementioned contained leak, one for dysphagia, and one post-operative pneumonia. Three patients with residual dysphagia underwent additional endoscopic procedures, all of whom had diverticula > 4 cm. Following myotomy, mean DI increased by 2.0 ± 1.7 mm2/mmHg (p < 0.001). In those with good clinical success, change in DI averaged + 1.6 ± 1.1 mm2/mmHg. Significant improvement was found in RSI and GERD-HRQL scores, but not dysphagia score. CONCLUSION: Z-POEM is a safe and feasible for treatment of ZD. We saw zero cases of intraoperative abandonment. We propose that large diverticula (> 4 cm) are a risk factor for poor outcomes and may require additional endoscopic procedures. An improvement in DI is expected after myotomy, however, the ideal range is still not known.

Genomic analysis of the brassica pathogen turnip mosaic potyvirus reveals its spread along the former trade routes of the Silk Road.

Plant pathogens have agricultural impacts on a global scale and resolving the timing and route of their spread can aid crop protection and inform control strategies. However, the evolutionary and phylogeographic history of plant pathogens in Eurasia remains largely unknown because of the difficulties in sampling across such a large landmass. Here, we show that turnip mosaic potyvirus (TuMV), a significant pathogen of brassica crops, spread from west to east across Eurasia from about the 17th century CE. We used a Bayesian phylogenetic approach to analyze 579 whole genome sequences and up to 713 partial sequences of TuMV, including 122 previously unknown genome sequences from isolates that we collected over the past five decades. Our phylogeographic and molecular clock analyses showed that TuMV isolates of the Asian-Brassica/Raphanus (BR) and basal-BR groups and world-Brassica3 (B3) subgroup spread from the center of emergence to the rest of Eurasia in relation to the host plants grown in each country. The migration pathways of TuMV have retraced some of the major historical trade arteries in Eurasia, a network that formed the Silk Road, and the regional variation of the virus is partly characterized by different type patterns of recombinants. Our study presents a complex and detailed picture of the timescale and major transmission routes of an important plant pathogen.

Multiple migrations to the Philippines during the last 50,000 years.

Island Southeast Asia has recently produced several surprises regarding human history, but the region's complex demography remains poorly understood. Here, we report ∼2.3 million genotypes from 1,028 individuals representing 115 indigenous Philippine populations and genome-sequence data from two ∼8,000-y-old individuals from Liangdao in the Taiwan Strait. We show that the Philippine islands were populated by at least five waves of human migration: initially by Northern and Southern Negritos (distantly related to Australian and Papuan groups), followed by Manobo, Sama, Papuan, and Cordilleran-related populations. The ancestors of Cordillerans diverged from indigenous peoples of Taiwan at least ∼8,000 y ago, prior to the arrival of paddy field rice agriculture in the Philippines ∼2,500 y ago, where some of their descendants remain to be the least admixed East Asian groups carrying an ancestry shared by all Austronesian-speaking populations. These observations contradict an exclusive "out-of-Taiwan" model of farming-language-people dispersal within the last four millennia for the Philippines and Island Southeast Asia. Sama-related ethnic groups of southwestern Philippines additionally experienced some minimal South Asian gene flow starting ∼1,000 y ago. Lastly, only a few lowlanders, accounting for <1% of all individuals, presented a low level of West Eurasian admixture, indicating a limited genetic legacy of Spanish colonization in the Philippines. Altogether, our findings reveal a multilayered history of the Philippines, which served as a crucial gateway for the movement of people that ultimately changed the genetic landscape of the Asia-Pacific region.

Escape from cell-cell and cell-matrix adhesion dependence underscores disease progression in gastric cancer organoid models.

OBJECTIVE: Cell-cell (CC) and cell-matrix (CM) adhesions are essential for epithelial cell survival, yet dissociation-induced apoptosis is frequently circumvented in malignant cells. DESIGN: We explored CC and CM dependence in 58 gastric cancer (GC) organoids by withdrawing either ROCK inhibitor, matrix or both to evaluate their tumorigenic potential in terms of apoptosis resistance, correlation with oncogenic driver mutations and clinical behaviour. We performed mechanistic studies to determine the role of diffuse-type GC drivers: ARHGAP fusions, RHOA and CDH1, in modulating CC (CCi) or CM (CMi) adhesion independence. RESULTS: 97% of the tumour organoids were CMi, 66% were CCi and 52% were resistant to double withdrawal (CCi/CMi), while normal organoids were neither CMi nor CCi. Clinically, the CCi/CMi phenotype was associated with an infiltrative tumour edge and advanced tumour stage. Moreover, the CCi/CMi transcriptome signature was associated with poor patient survival when applied to three public GC datasets. CCi/CMi and CCi phenotypes were enriched in diffuse-type GC organoids, especially in those with oncogenic driver perturbation of RHO signalling via RHOA mutation or ARHGAP fusions. Inducible knockout of ARHGAP fusions in CCi/CMi tumour organoids led to resensitisation to CC/CM dissociation-induced apoptosis, upregulation of focal adhesion and tight junction genes, partial reversion to a more normal cystic phenotype and inhibited xenograft formation. Normal gastric organoids engineered with CDH1 or RHOA mutations became CMi or CCi, respectively. CONCLUSIONS: The CCi/CMi phenotype has a critical role in malignant transformation and tumour progression, offering new mechanistic information on RHO-ROCK pathway inhibition that contributes to GC pathogenicity.

A novel thiol-saccharide mucolytic for the treatment of muco-obstructive lung diseases.

BACKGROUND: Mucin disulfide cross-links mediate pathologic mucus formation in muco-obstructive lung diseases. MUC-031, a novel thiol-modified carbohydrate compound, cleaves disulfides to cause mucolysis. The aim of this study was to determine the mucolytic and therapeutic effects of MUC-031 in sputum from patients with cystic fibrosis (CF) and mice with muco-obstructive lung disease (ßENaC-Tg mice). METHODS: We compared the mucolytic efficacy of MUC-031 and existing mucolytics (N-acetylcysteine (NAC) and recombinant human deoxyribonuclease I (rhDNase)) using rheology to measure the elastic modulus (G') of CF sputum, and we tested effects of MUC-031 on airway mucus plugging, inflammation and survival in ßENaC-Tg mice to determine its mucolytic efficacy in vivo. RESULTS: In CF sputum, compared to the effects of rhDNase and NAC, MUC-031 caused a larger decrease in sputum G', was faster in decreasing sputum G' by 50% and caused mucolysis of a larger proportion of sputum samples within 15 min of drug addition. Compared to vehicle control, three treatments with MUC-031 in 1 day in adult ßENaC-Tg mice decreased airway mucus content (16.8±3.2 versus 7.5±1.2 nL·mm-2, p<0.01) and bronchoalveolar lavage cells (73 833±6930 versus 47 679±7736 cells·mL-1, p<0.05). Twice-daily treatment with MUC-031 for 2 weeks also caused decreases in these outcomes in adult and neonatal ßENaC-Tg mice and reduced mortality from 37% in vehicle-treated ßENaC-Tg neonates to 21% in those treated with MUC-031 (p<0.05). CONCLUSION: MUC-031 is a potent and fast-acting mucolytic that decreases airway mucus plugging, lessens airway inflammation and improves survival in ßENaC-Tg mice. These data provide rationale for human trials of MUC-031 in muco-obstructive lung diseases.

Longitudinal effects of elexacaftor/tezacaftor/ivacaftor on sputum viscoelastic properties, airway infection and inflammation in patients with cystic fibrosis.

BACKGROUND: Recent studies demonstrated that the triple combination cystic fibrosis transmembrane conductance regulator (CFTR) modulator therapy elexacaftor/tezacaftor/ivacaftor (ETI) improves lung function and reduces pulmonary exacerbations in cystic fibrosis (CF) patients with at least one F508del allele. However, effects of ETI on downstream consequences of CFTR dysfunction, i.e. abnormal viscoelastic properties of airway mucus, chronic airway infection and inflammation have not been studied. The aim of this study was to determine the longitudinal effects of ETI on airway mucus rheology, microbiome and inflammation in CF patients with one or two F508del alleles aged ≥12 years throughout the first 12 months of therapy. METHODS: In this prospective observational study, we assessed sputum rheology, the microbiome, inflammation markers and proteome before and 1, 3 and 12 months after initiation of ETI. RESULTS: In total, 79 patients with CF and at least one F508del allele and 10 healthy controls were enrolled in this study. ETI improved the elastic modulus and viscous modulus of CF sputum at 3 and 12 months after initiation (all p<0.01). Furthermore, ETI decreased the relative abundance of Pseudomonas aeruginosa in CF sputum at 3 months and increased the microbiome α-diversity at all time points. In addition, ETI reduced interleukin-8 at 3 months (p<0.05) and free neutrophil elastase activity at all time points (all p<0.001), and shifted the CF sputum proteome towards healthy. CONCLUSIONS: Our data demonstrate that restoration of CFTR function by ETI improves sputum viscoelastic properties, chronic airway infection and inflammation in CF patients with at least one F508del allele over the first 12 months of therapy; however, levels close to healthy were not reached.

Genomic analysis reveals strong population structure in the Giant Sydney Crayfish (Euastacus spinifer (Heller, 1865)).

Australia is home to over 140 species of freshwater crayfish (Decapoda: Parastacidae), representing a centre of diversity for this group in the Southern Hemisphere. Species delimitation in freshwater crayfish is difficult because many species show significant variation in colouration and morphology. This is particularly evident in the genus Euastacus, which exhibits large variations in colour and spination throughout its putative range. To understand this variation, we investigated the genetic diversity, population structure, phylogeny, and evolutionary timescale of the Giant Sydney Crayfish (Euastacus spinifer (Heller, 1865)). Our data set is sampled from over 70 individuals from across the ∼600 km range of the species, and includes a combination of two mitochondrial markers and more than 7000 single-nucleotide polymorphisms (SNPs) from the nuclear genome. Data were also obtained for representatives of the close relative, Euastacus vesper McCormack and Ahyong, 2017. Genomic SNP analyses revealed strong population structure, with multiple distinct populations showing little evidence of gene flow or migration. Phylogenetic analyses of mitochondrial data revealed similar structure between populations. Taken together, our analyses suggest that E. spinifer, as currently understood, represents a species complex, of which E. vesper is a member. Molecular clock estimates place the divergences within this group during the Pleistocene. The isolated and highly fragmented populations identified in our analyses probably represent relict populations of a previously widespread ancestral species. Periodic flooding events during the Pleistocene are likely to have facilitated the movement of these otherwise restricted freshwater crayfish within and between drainage basins, including the Murray-Darling and South East Coast Drainages. We present evidence supporting the recognition of populations in the southern parts of the range of E. spinifer as one or two separate species, which would raise the number of species within the E. spinifer complex to at least three. Our results add to the growing body of evidence that many freshwater crayfish exhibit highly fragmented, range-restricted distributions. In combination with the life-history traits of these species, the restricted distributions exacerbate the threats already placed on freshwater crayfish, which are among the five most endangered animal groups globally.

Excluding Loci With Substitution Saturation Improves Inferences From Phylogenomic Data.

The historical signal in nucleotide sequences becomes eroded over time by substitutions occurring repeatedly at the same sites. This phenomenon, known as substitution saturation, is recognized as one of the primary obstacles to deep-time phylogenetic inference using genome-scale data sets. We present a new test of substitution saturation and demonstrate its performance in simulated and empirical data. For some of the 36 empirical phylogenomic data sets that we examined, we detect substitution saturation in around 50% of loci. We found that saturation tends to be flagged as problematic in loci with highly discordant phylogenetic signals across sites. Within each data set, the loci with smaller numbers of informative sites are more likely to be flagged as containing problematic levels of saturation. The entropy saturation test proposed here is sensitive to high evolutionary rates relative to the evolutionary timeframe, while also being sensitive to several factors known to mislead phylogenetic inference, including short internal branches relative to external branches, short nucleotide sequences, and tree imbalance. Our study demonstrates that excluding loci with substitution saturation can be an effective means of mitigating the negative impact of multiple substitutions on phylogenetic inferences. [Phylogenetic model performance; phylogenomics; substitution model; substitution saturation; test statistics.].

Evolutionary Rate Variation among Lineages in Gene Trees has a Negative Impact on Species-Tree Inference.

Phylogenetic analyses of genomic data provide a powerful means of reconstructing the evolutionary relationships among organisms, yet such analyses are often hindered by conflicting phylogenetic signals among loci. Identifying the signals that are most influential to species-tree estimation can help to inform the choice of data for phylogenomic analysis. We investigated this in an analysis of 30 phylogenomic data sets. For each data set, we examined the association between several branch-length characteristics of gene trees and the distance between these gene trees and the corresponding species trees. We found that the distance of each gene tree to the species tree inferred from the full data set was positively associated with variation in root-to-tip distances and negatively associated with mean branch support. However, no such associations were found for gene-tree length, a measure of the overall substitution rate at each locus. We further explored the usefulness of the best-performing branch-based characteristics for selecting loci for phylogenomic analyses. We found that loci that yield gene trees with high variation in root-to-tip distances have a disproportionately distant signal of tree topology compared with the complete data sets. These results suggest that rate variation across lineages should be taken into consideration when exploring and even selecting loci for phylogenomic analysis.[Branch support; data filtering; nucleotide substitution model; phylogenomics; substitution rate; summary coalescent methods.].