Timing and Origins of Local and Distant Metastases in Lung Cancer

INTRODUCTION

Metastasis is the primary cause of lung cancer-related death. Nevertheless, the underlying molecular mechanisms and evolutionary patterns of lung cancer metastases are still elusive.

METHODS

We performed whole-exome sequencing for 40 primary tumors (PTs) and 61 metastases from 47 patients with lung cancer, of which 40 patients had paired PTs and metastases. The PT-metastasis genomic divergence, metastatic drivers, timing of metastatic dissemination, and evolutionary origins were analyzed using appropriate statistical tools and mathematical models.

RESULTS

There were various degrees of genomic heterogeneity when comparing the paired primary and metastatic lesions or comparing metastases of different sites. Multiple metastasis-selected/enriched genetic alterations were found, such as MYC amplification, NKX2-1 amplification, RICTOR amplification, arm 20p gain, and arm 11p loss, and these results were were also featured in a meta-analysis cross-validated using an independent cohort from Memorial Sloan-Kettering Cancer Center database. To elucidate the metastatic seeding time, we applied a metastatic model and found 61.1% of the tumors were late dissemination, in which the metastatic seeding happened approximately 2.74 years before clinical detection. One exception was lymph node metastases whose dissemination time was relatively early. By analyzing the evolutionary origins, we reported that nonlymph node metastases were mainly seeded by the PT (87.5%) rather than the earlier colonized lymph node metastases.

CONCLUSIONS

Our results shed light on the molecular features that potentially drive lung cancer metastases. The distinct temporospatial pattern of disease progression revealed that lung cancer was susceptible to either late dissemination or indolent early lymph node metastases, leaving a potential time window to minimize metastases by early cancer detection.

Functions of intrinsically disordered proteins through evolutionary lenses

Protein sequences are the result of an evolutionary process that involves the balancing act of experimenting with novel mutations and selecting out those that have an undesirable functional outcome. In the case of globular proteins, the function relies on a well-defined conformation, therefore, there is a strong evolutionary pressure to preserve the structure. However, different evolutionary rules might apply for the group of intrinsically disordered regions and proteins (IDR/IDPs) that exist as an ensemble of fluctuating conformations. The function of IDRs can directly originate from their disordered state or arise through different types of molecular recognition processes. There is an amazing variety of ways IDRs can carry out their functions, and this is also reflected in their evolutionary properties. In this chapter we give an overview of the different types of evolutionary behavior of disordered proteins and associated functions in normal and disease settings.

Microbial degradation kinetics and molecular mechanism of 2,6-dichloro-4-nitrophenol by a Cupriavidus strain

2,6-Dichloro-4-nitrophenol (2,6-DCNP) is an emerging chlorinated nitroaromatic pollutant, and its fate in the environment is an important question. However, microorganisms with the ability to utilize 2,6-DCNP have not been reported. In this study, Cupriavidus sp. CNP-8 having been previously reported to degrade various halogenated nitrophenols, was verified to be also capable of degrading 2,6-DCNP. Biodegradation kinetics assay showed that it degraded 2,6-DCNP with the specific growth rate of 0.124 h^-1, half saturation constant of 0.038 mM and inhibition constant of 0.42 mM. Real-time quantitative PCR analyses indicated that the hnp gene cluster was involved in the catabolism of 2,6-DCNP. The hnpA and hnpB gene products were purified to homogeneity by Ni-NTA chromatography. Enzymatic assays showed that HnpAB, a FAD-dependent two-component monooxygenase, converted 2,6-DCNP to 6-chlorohydroxyquinol with a K_m of 3.9 ± 1.4 μM and a k_cat/K_m of 0.12 ± 0.04 μΜ^-1 min^-1. As the oxygenase component encoding gene, hnpA is necessary for CNP-8 to grow on 2,6-DCNP by gene knockout and complementation. The phylogenetic analysis showed that the hnp cluster originated from the cluster involved in the catabolism of chlorophenols rather than nitrophenols. To our knowledge, CNP-8 is the first bacterium with the ability to utilize 2,6-DCNP, and this study fills a gap in the microbial degradation mechanism of this pollutant at the molecular, biochemical and genetic levels. Moreover, strain CNP-8 could degrade three chlorinated nitrophenols rapidly from the synthetic wastewater, indicating its potential in the bioremediation of chlorinated nitrophenols polluted environments.

Morphology and diversity of the forcipules in Strigamia centipedes (Chilopoda, Geophilomorpha)

The morphology of the venomous limbs (forcipules) of 13 species of Strigamia and of six other geophilomorphs was studied with light microscopy, scanning electron microscopy, and, for a subsample, with confocal laser scanning microscopy. In all Strigamia species a well-distinct denticle is present invariantly on the inner side of the terminal article (tarsungulum), in sub-basal position, just proximal to a faint transverse sulcus and a cuticular introflexion that corresponds to the insertion point of a tendon. Strigamia species differ mainly in size and shape of the denticle and thickness of the distal part of the tarsungulum, suggesting some functional diversity in piercing and handling prey. Anatomical evidence supports the hypothesis that the tarsungulum corresponds to two ancestral articles and a denticle at the basis of the tarsungulum originated multiple times within geophilomorphs, however in different positions corresponding to either the ancestral sub-terminal article (in Strigamia, other Geophiloidea and some Schendylidae) or the ancestral terminal article (in the himantariid Thracophilus).

Shigella flexneri serotype 1c derived from serotype 1a by acquisition of gtrIC gene cluster via a bacteriophage

Background

Shigella spp. are the primary causative agents of bacillary dysentery. Since its emergence in the late 1980s, the S. flexneri serotype 1c remains poorly understood, particularly with regard to its origin and genetic evolution. This article provides a molecular insight into this novel serotype and the gtrIC gene cluster that determines its unique immune recognition.

Results

A PCR of the gtrIC cluster showed that serotype 1c isolates from different geographical origins were genetically conserved. An analysis of sequences flanking the gtrIC cluster revealed remnants of a prophage genome, in particular integrase and tRNA^Pro genes. Meanwhile, Southern blot analyses on serotype 1c, 1a and 1b strains indicated that all the tested serotype 1c strains may have had a common origin that has since remained distinct from the closely related 1a and 1b serotypes. The identification of prophage genes upstream of the gtrIC cluster is consistent with the notion of bacteriophage-mediated integration of the gtrIC cluster into a pre-existing serotype.

Conclusions

This is the first study to show that serotype 1c isolates from different geographical origins share an identical pattern of genetic arrangement, suggesting that serotype 1c strains may have originated from a single parental strain. Analysis of the sequence around the gtrIC cluster revealed a new site for the integration of the serotype converting phages of S. flexneri. Understanding the origin of new pathogenic serotypes and the molecular basis of serotype conversion in S. flexneri would provide information for developing cross-reactive Shigella vaccines.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-016-0746-z) contains supplementary material, which is available to authorized users.

Progress in understanding the phylogeny of the Plasmodium vivax lineage

There has been some controversy about the evolutionary origin of Plasmodium vivax, particularly whether it is of Asian or African origin. Recently, a new malaria species which closely related to ape P. vivax was found in chimpanzees, in addition, the host switches of P. vivax from ape to human was confirmed. These findings support the African origin of P. vivax. Previous phylogenetic analyses have shown the position of P. vivax within the Asian primate malaria parasite clade. This suggested an Asian origin of P. vivax. Recent analyses using massive gene data, however, positioned P. vivax after the branching of the African Old World monkey parasite P. gonderi, and before the branching of the common ancestor of Asian primate malaria parasites. This position is consistent with an African origin of P. vivax. We here review the history of phylogenetic analyses on P. vivax, validate previous analyses, and finally present a definitive analysis using currently available data that indicate a tree in which P. vivax is positioned at the base of the Asian primate malaria parasite clade, and thus that is consistent with an African origin of P. vivax.

Evolutionary relationships between seryl-histidine dipeptide and modern serine proteases from the analysis based on mass spectrometry and bioinformatics

Seryl-histidine dipeptide (Ser-His) has been recognized as the shortest peptide with hydrolysis cleavage activity; however, its protein cleavage spectrum has not yet been fully explored. Here, four differently folded proteins were treated with Ser-His, and the digestion products were evaluated with high-resolution mass spectrometry. The cleavage efficiency and cleavage propensity of Ser-His against these protein substrates were calculated at both the primary and secondary sequence levels. The above experiments show that Ser-His cleaves a broad spectrum of substrate proteins of varying secondary structures. Moreover, Ser-His could cleave at all 20 amino acids with different efficiencies according to the protein, which means that Ser-His has the original digestion function of serine proteases. Furthermore, we collected and compared the catalytic sites and cleavage sites of 340 extant serine proteases derived from 17 representative organisms. A consensus motif Ser-[X]-His was identified as the major pattern at the catalytic sites of serine proteases from all of the organisms represented except Danio rerio, which uses Ser-Lys instead. This finding indicates that Ser-His is the core component of the serine protease catalytic site. Moreover, our analysis revealed that the cleavage sites of modern serine proteases have become more specific over the evolutionary history of this family. Based on the above analysis results, it could be found that Ser-His is likely the original serine protease and maybe the evolutionary core of modern serine proteases.

Mitogenomic phylogeny and divergence time estimation of Artemia Leach, 1819 (Branchiopoda: Anostraca) with emphasis on parthenogenetic lineages

The brine shrimp Artemia, a crustacean adapted to the extreme conditions of hypersaline environments, comprises nine regionally distributed sexual species scattered (island-like) over heterogeneous environments and asexual (parthenogenetic) lineages with different ploidies. Such sexual and asexual interaction within the genus raises questions regarding the origin and time of divergence of both sexual species and asexual lineages, including the persistence of the latter over time, a problem not yet clarified using single mitochondrial and nuclear markers. Based on the complete mitochondrial genome of all species and parthenogenetic lineages, this article first describes the mitogenomic characteristics (nucleotide compositions, genome mapping, codon usage, and tRNA secondary structure) of sexual species and asexual types and, secondly, it provides a comprehensive updated phylogenetic analysis. Molecular dating and geographical evidence suggest that the ancestral Artemia taxon originated in ca. 33.97 Mya during the Paleogene Period. The mitogenomic comparisons suggest that the common ancestor of diploid and triploid parthenogenetic lineages (ca. 0.07 Mya) originated from a historical ancestor (ca. 0.61 Mya) in the Late Pleistocene. Additionally, the common ancestor of tetraploid and pentaploid parthenogenetic lineages (ca. 0.05 Mya) diverged from a historical maternal ancestor with A. sinica (ca. 0.96 Mya) in the early Pleistocene. The parthenogenetic lineages do not share a direct ancestor with any sexual species. The Asian clade ancestor diverged more recently (ca. 14.27 Mya, Middle Miocene). The mitogenomic characteristics, maternal phylogenetic tree, and especially divergence time prove that A. monica and A. franciscana are two biological species.

Pathogenic variants in human DNA damage repair genes mostly arose in recent human history

BACKGROUND

Genome stability is maintained by the DNA damage repair (DDR) system composed of multiple DNA repair pathways of hundreds of genes. Germline pathogenic variation (PV) in DDR genes damages function of the affected DDR genes, leading to genome instability and high risk of diseases, in particular, cancer. Knowing evolutionary origin of the PVs in human DDR genes is essential to understand the etiology of human diseases. However, answer to the issue remains largely elusive. In this study, we analyzed evolutionary origin for the PVs in human DDR genes.

METHODS

We identified 169 DDR genes by referring to various databases and identified PVs in the DDR genes of modern humans from ClinVar database. We performed a phylogenetic analysis to analyze the conservation of human DDR PVs in 100 vertebrates through cross-species genomic data comparison using the phyloFit program of the PHAST package and visualized the results using the GraphPad Prism software and the ggplot module. We identified DDR PVs from over 5000 ancient humans developed a database to host the DDR PVs ( https://genemutation.fhs.um.edu.mo/dbDDR-AncientHumans ). Using the PV data, we performed a molecular archeological analysis to compare the DDR PVs between modern humans and ancient humans. We analyzed evolution selection of DDR genes across 20 vertebrates using the CodeML in PAML for phylogenetic analysis.

RESULTS

Our phylogenic analysis ruled out cross-species conservation as the origin of human DDR PVs. Our archeological approach identified rich DDR PVs shared between modern and ancient humans, which were mostly dated within the last 5000 years. We also observed similar pattern of quantitative PV distribution between modern and ancient humans. We further detected a set of ATM, BRCA2 and CHEK2 PVs shared between human and Neanderthals.

CONCLUSIONS

Our study reveals that human DDR PVs mostly arose in recent human history. We propose that human high cancer risk caused by DDR PVs can be a by-product of human evolution.

Evidence for cross-species transmission of human coronavirus OC43 through bioinformatics and modeling infections in porcine intestinal organoids

Cross-species transmission of coronaviruses has been continuously posing a major challenge to public health. Pigs, as the major animal reservoirs for many zoonotic viruses, frequently mediate viral transmission to humans. This study comprehensively mapped the relationship between human and porcine coronaviruses through in-depth bioinformatics analysis. We found that human coronavirus OC43 and porcine coronavirus PHEV share a close phylogenetic relationship, evidenced by high genomic homology, similar codon usage patterns and comparable tertiary structure in spike proteins. Inoculation of infectious OC43 viruses in organoids derived from porcine small and large intestine demonstrated that porcine intestinal organoids (pIOs) are highly susceptible to human coronavirus OC43 infection and support infectious virus production. Using transmission electron microscopy, we visualized OC43 viral particles in both intracellular and extracellular compartments, and observed abnormalities of multiple organelles in infected organoid cells. Robust OC43 infections in pIOs result in a significant reduction of organoids viability and widespread cell death. This study bears essential implications for better understanding the evolutionary origin of human coronavirus OC43, and provides a proof-of-concept for using pIOs as a model to investigate cross-species transmission of human coronavirus.

Large-scale analysis of the ARF and Aux/IAA gene families in 406 horticultural and other plants

The auxin response factor (ARF) and auxin/indole-3-acetic acid (Aux/IAA) family of genes are central components of the auxin signaling pathway and play essential roles in plant growth and development. Their large-scale analysis and evolutionary trajectory of origin are currently not known. Here, we identified the corresponding ARF and Aux/IAA family members and performed a large-scale analysis by scanning 406 plant genomes. The results showed that the ARF and Aux/IAA gene families originated from charophytes. The ARF family sequences were more conserved than the Aux/IAA family sequences. Dispersed duplications were the common expansion mode of ARF and Aux/IAA families in bryophytes, ferns, and gymnosperms; however, whole-genome duplication was the common expansion mode of the ARF and Aux/IAA families in basal angiosperms, magnoliids, monocots, and dicots. Expression and regulatory network analyses revealed that the Arabidopsis thaliana ARF and Aux/IAA families responded to multiple hormone, biotic, and abiotic stresses. The APETALA2 and serum response factor-transcription factor gene families were commonly enriched in the upstream and downstream genes of the ARF and Aux/IAA gene families. Our study provides a comprehensive overview of the evolutionary trajectories, structural functions, expansion mechanisms, expression patterns, and regulatory networks of these two gene families.

Evolutionary origin of germline pathogenic variants in human DNA mismatch repair genes

BACKGROUND

Mismatch repair (MMR) system is evolutionarily conserved for genome stability maintenance. Germline pathogenic variants (PVs) in MMR genes that lead to MMR functional deficiency are associated with high cancer risk. Knowing the evolutionary origin of germline PVs in human MMR genes will facilitate understanding the biological base of MMR deficiency in cancer. However, systematic knowledge is lacking to address the issue. In this study, we performed a comprehensive analysis to know the evolutionary origin of human MMR PVs.

METHODS

We retrieved MMR gene variants from the ClinVar database. The genomes of 100 vertebrates were collected from the UCSC genome browser and ancient human sequencing data were obtained through comprehensive data mining. Cross-species conservation analysis was performed based on the phylogenetic relationship among 100 vertebrates. Rescaled ancient sequencing data were used to perform variant calling for archeological analysis.

RESULTS

Using the phylogenetic approach, we traced the 3369 MMR PVs identified in modern humans in 99 non-human vertebrate genomes but found no evidence for cross-species conservation as the source for human MMR PVs. Using the archeological approach, we searched the human MMR PVs in over 5000 ancient human genomes dated from 45,045 to 100 years before present and identified a group of MMR PVs shared between modern and ancient humans mostly within 10,000 years with similar quantitative patterns.

CONCLUSION

Our study reveals that MMR PVs in modern humans were arisen within the recent human evolutionary history.

Relation Between Obesity and Type 2 Diabetes: Evolutionary Insights, Perspectives and Controversies

PURPOSE OF REVIEW

Since the mid-twentieth century, obesity and its related comorbidities, notably insulin resistance (IR) and type 2 diabetes (T2D), have surged. Nevertheless, their underlying mechanisms remain elusive. Evolutionary medicine (EM) sheds light on these issues by examining how evolutionary processes shape traits and diseases, offering insights for medical practice. This review summarizes the pathogenesis and genetics of obesity-related IR and T2D. Subsequently, delving into their evolutionary connections. Addressing limitations and proposing future research directions aims to enhance our understanding of these conditions, paving the way for improved treatments and prevention strategies.

RECENT FINDINGS

Several evolutionary hypotheses have been proposed to unmask the origin of obesity-related IR and T2D, e.g., the "thrifty genotype" hypothesis suggests that certain "thrifty genes" that helped hunter-gatherer populations efficiently store energy as fat during feast-famine cycles are now maladaptive in our modern obesogenic environment. The "drifty genotype" theory suggests that if thrifty genes were advantageous, they would have spread widely, but proposes genetic drift instead. The "behavioral switch" and "carnivore connection" hypotheses propose insulin resistance as an adaptation for a brain-dependent, low-carbohydrate lifestyle. The thrifty phenotype theory suggests various metabolic outcomes shaped by genes and environment during development. However, the majority of these hypotheses lack experimental validation. Understanding why ancestral advantages now predispose us to diseases may aid in drug development and prevention of disease. EM helps us to understand the evolutionary relation between obesity-related IR and T2D. But still gaps and contradictions persist. Further interdisciplinary research is required to elucidate complete mechanisms.

Evolution, genetic recombination, and phylogeography of goose parvovirus

Goose parvovirus (GPV) has garnered global attention due to its association with severe symptoms in waterfowl. However, the process underlying the global emergence and spread of GPV remains largely elusive. In this study, we illustrated the evolutionary characteristics of GPVs from a global perspective using phylogenetic analysis, recombination analysis, selection pressure analysis, and phylogeographic analysis. Our findings indicate that GPV and muscovy duck parvovirus (MDPV) diverge into two distinct branches. Within GPV, there are two classifications: classical GPV (C-GPV) and novel GPV (N-GPV), each containing three subgroups, underscoring the significant genetic diversity of GPV. Recombination analysis revealed 11 recombination events, suggesting C-GPV, N-GPV, and MDPV co-infections. Further, phylogeographic analysis revealed that China is an important exporter of GPV and that trade might serve as a potential transmission conduit. Nonetheless, a detailed understanding of its geographic transmission dynamics warrants further investigation due to the limited scope of current genomic data in our study. This study offers novel insights into the evolutionary state and spread of GPV, holding promise for informing preventive and containment strategies against GPV infection.

Evolutionary origin and gradual accumulation with plant evolution of the LACS family

BACKGROUND

LACS (long-chain acyl-CoA synthetase) genes are widespread in organisms and have multiple functions in plants, especially in lipid metabolism. However, the origin and evolutionary dynamics of the LACS gene family remain largely unknown.

RESULTS

Here, we identified 1785 LACS genes in the genomes of 166 diverse plant species and identified the clades (I, II, III, IV, V, VI) of six clades for the LACS gene family of green plants through phylogenetic analysis. Based on the evolutionary history of plant lineages, we found differences in the origins of different clades, with Clade IV originating from chlorophytes and representing the origin of LACS genes in green plants. The structural characteristics of different clades indicate that clade IV is relatively independent, while the relationships between clades (I, II, III) and clades (V, VI) are closer. Dispersed duplication (DSD) and transposed duplication (TRD) are the main forces driving the evolution of plant LACS genes. Network clustering analysis further grouped all LACS genes into six main clusters, with genes within each cluster showing significant co-linearity. Ka/Ks results suggest that LACS family genes underwent purifying selection during evolution. We analyzed the phylogenetic relationships and characteristics of six clades of the LACS gene family to explain the origin, evolutionary history, and phylogenetic relationships of different clades and proposed a hypothetical evolutionary model for the LACS family of genes in plants.

CONCLUSIONS

Our research provides genome-wide insights into the evolutionary history of the LACS gene family in green plants. These insights lay an important foundation for comprehensive functional characterization in future research.

Evolutionary history of Plasmodium vivax and Plasmodium simium in the Americas

Malaria is a vector-borne disease caused by protozoan parasites of the genus Plasmodium. Plasmodium vivax is the most prevalent human-infecting species in the Americas. However, the origins of this parasite in this continent are still debated. Similarly, it is now accepted that the existence of Plasmodium simium is explained by a P. vivax transfer from humans to monkey in America. However, many uncertainties still exist concerning the origin of the transfer and whether several transfers occurred. In this review, the most recent studies that addressed these questions using genetic and genomic approaches are presented.