Assessment of the population structure and temporal changes in spatial dynamics and genetic characteristics of the Atlantic bluefin tuna under a fishery independent framework.

As a large and long-lived species with high economic value, restricted spawning areas and short spawning periods, the Atlantic bluefin tuna (BFT; Thunnus thynnus) is particularly susceptible to over-exploitation. Although BFT have been targeted by fisheries in the Mediterranean Sea for thousands of years, it has only been in these last decades that the exploitation rate has reached far beyond sustainable levels. An understanding of the population structure, spatial dynamics, exploitation rates and the environmental variables that affect BFT is crucial for the conservation of the species. The aims of this PhD project were 1) to assess the accuracy of larval identification methods, 2) determine the genetic structure of modern BFT populations, 3) assess the self-recruitment rate in the Gulf of Mexico and Mediterranean spawning areas, 4) estimate the immigration rate of BFT to feeding aggregations from the various spawning areas, and 5) develop tools capable of investigating the temporal stability of population structuring in the Mediterranean Sea. Several weaknesses in modern morphology-based taxonomy including demographic decline of expert taxonomists, flawed identification keys, reluctance of the taxonomic community to embrace advances in digital communications and a general scarcity of modern user-friendly materials are reviewed. Barcoding of scombrid larvae revealed important differences in the accuracy of the taxonomic identifications carried out by different ichthyoplanktologists following morphology-based methods. Using a Genotyping-by-Sequencing a panel of 95 SNPs was developed and used to characterize the population structuring of BFT and composition of adult feeding aggregations. Using novel molecular techniques, DNA was extracted from bluefin tuna vertebrae excavated from late iron age, ancient roman settlements Byzantine-era Constantinople and a 20th century collection. A second panel of 96 SNPs was developed to genotype historical and modern samples in order to elucidate changes in population structuring and allele frequencies of loci associated with selective traits.

Caratterizzazione genomica di virus influenzali suini in Italia mediante next generation sequencing

I sottotipi H1N1, H1N2 e H3N2 di influenza A virus sono largamente diffusi nella popolazione suina di tutto il mondo. Nel presente lavoro è stato sviluppato un protocollo di sequenziamento di c.d. nuova generazione, su piattaforma Ion Torrent PGM, idoneo per l’analisi di tutti i virus influenzali suini (SIV). Per valutare l’evoluzione molecolare dei SIV italiani, sono stati sequenziati ed analizzati mediante analisi genomica e filogenetica un totale di sessantadue ceppi di SIV appartenenti ai sottotipi H1N1, H1N2 e H3N2, isolati in Italia dal 1998 al 2014. Sono stati evidenziati in sei campioni due fenomeni di riassortimento: tutti i SIV H1N2 esaminati presentavano una neuraminidasi di derivazione umana, diversa da quella dei SIV H1N2 circolanti in Europa, inoltre l’emoagglutinina (HA) di due isolati H1N2 era originata dal riassortimento con un SIV H1N1 avian-like. L’analisi molecolare dell’HA ha permesso di rivelare un’inserzione di due amminoacidi in quattro SIV H1N1 pandemici e una delezione di due aminoacidi in quattro SIV H1N2, entrambe a livello del sito di legame con il recettore cellulare. E’ stata inoltre evidenziata un’elevata omologia di un SIV H1N1 con ceppi europei isolati negli anni ’80, suggerendo la possibile origine vaccinale di questo virus. E’ stato possibile, in aggiunta, applicare il nuovo protocollo sviluppato per sequenziare un virus influenzale aviare altamente patogeno trasmesso all’uomo, direttamente da campione biologico. La diversità genetica nei SIV esaminati in questo studio conferma l’importanza di un continuo monitoraggio della costellazione genomica dei virus influenzali nella popolazione suina.

Indagine sulla presenza di Helicobacter pullorum in allevamenti avicoli italiani

From September 2005 to December 2006, in order to define the prevalence of Helicobacter pullorum in broiler chickens, laying hens and turkey, a total of 365 caecum contents of animals reared in 76 different farms were collected at the slaughterhouse. A caecum content of a ostrich was also sampled. In addition, with the aim of investigating the occurrence of H. pullorum in humans, 151 faeces were collected at the Sant’Orsola-Malpighi University Hospital of Bologna from patients suffering of gastroenteritis. A modified Steele–McDermott membrane filter method was used. Gram-negative curved rod bacteria were preliminary identified as H. pullorum by a PCR assay based on 16S rRNA, then subjected to a RFLP-PCR assay to distinguish between H. pullorum and H. canadensis. One isolate from each farm was randomly selected for phenotypic characterization by biochemical methods and 1D SDSPAGE analysis of whole cell proteins profiles. Minimum Inhibitory Concentration (MIC) for seven different antibiotics were also determined by agar dilution method. Moreover, to examine the intraspecific genomic variability, two strains isolated from 17 different farms were submitted to genotyping by Pulse-Field Gel Electrophoresis (PFGE). In order to assess the molecular basis of fluorquinolone resistance in H. pullorum, gyrA of H. pullorum CIP 104787T was sequenced and nucleotide sequences of the Quinolone Resistance Determining Region (QRDR) of a total of 18 poultry isolates, with different MIC values for ciprofloxacin and nalidixic acid, were compared. According to the PCR and PCR-RFLP results, 306 out of 366 animals examined were positive for H. pullorum (83,6%) and 96,1% of farms resulted infected. All positive samples showed a high number of colonies (>50) phenotipically consistent with H. pullorum on the first isolation media, which suggests that this microrganism, when present, colonizes the poultry caecum at an elevate load. No human sample resulted positive for H. pullorum. The 1D SDS-PAGE whole protein profile analysis showed high similarity among the 74 isolates tested and with the type strain H. pullorum CIP 104787T. Regarding the MIC values, a monomodal distribution was found for ampicillin, chloramphenicol, gentamicin and nalidixic acid, whereas a bimodal trend was noticed for erythromycin, ciprofloxacin and tetracycline (indicating an acquired resistance for these antibiotics). Applying the breakpoints indicated by the CSLI, we may assume that all the H. pullorum tested are sensitive only to gentamicin. The intraspecific genomic variability observed in this study confirm that this species don’t have a clonal population structure, as motioned by other autors. The 2490 bp gyrA gene of H. pullorum CIP104787T with an Open Reading Frame (ORF) encoding a polypeptide of 829 amino acids was for the first time sequenced and characterized. All ciprofloxacin resistant poultry isolates showed ACA®ATA (Thr®Ile) substitution at codon 84 of gyrA corresponding to codons of gyrA 86, 87 and 83 of the Campylobacter jejuni, H. pylori and Escherichia coli, respectively. This substitution was functionally confirmed to be associated with the ciprofloxacin resistant phenotype of poultry isolates. This is the first report of isolation of H. pullorum in turkey and in ostrich, indicating that poultry species are the reservoir of this potential zoonotic microorganisms. In order to understand the potential role as food-borne human pathogen of H. pullorum, further studies must be carried on.

Characterization and cloning of chemically-induced mutants in barley (Hordeum vulgare L.)

Induced mutagenesis has been exploited for crop improvement and for investigating gene function and regulation. To unravel molecular mechanisms of stress resilience, we applied state-of-the-art genomics-based gene cloning methods to barley mutant lines showing altered root and shoot architecture and disease lesion mimic phenotypes. With a novel method that we named complementation by sequencing, we cloned NEC3, the causal gene for an orange-spotted disease lesion mimic phenotype. NEC3 belongs to the CYP71P1 gene family and it is involved in serotonin biosynthesis. By comparative phylogenetic analysis we showed that CYP71P1 emerged early in angiosperm evolution but was lost in some lineages including Arabidopsis thaliana. By BSA-Seq, we cloned the gene whose mutation increased leaf width, and we showed that the gene corresponded to the previously cloned BROADLEAF1. By BSA coupled to WGS sequencing, we cloned EGT1 and EGT2, two genes that regulate root gravitropic set point angle. EGT1 encodes a Tubby-like F-box protein and EGT2 encodes a Sterile Alpha Motive protein; EGT2 is phylogenetically related to AtSAM5 in Arabidopsis and to WEEP in peach where it regulates branch angle. Both EGT1 and EGT2 are conserved in wheat. We hypothesized that both participate to an anti-gravitropic offset mechanism since their disruption causes mutant roots to grow along the gravity vector. By the MutMap+ method, we cloned the causal gene of a short and semi-rigid root mutant and found that it encodes for an endoglucanase and is the ortholog of OsGLU3 in rice whose mutant has the same phenotype, suggesting that the gene is conserved in barley and rice. The mutants and the corresponding genes which were cloned in this work are involved in the response to stress and can potentially contribute to crop adaptation.

Identification of novel genetic alterations in pediatric cytogenetically normal acute myeloid leukemia by next-generation sequencing

Pediatric acute myeloid leukemia (AML) is a molecularly heterogeneous disease that arises from genetic alterations in pathways that regulate self-renewal and myeloid differentiation. While the majority of patients carry recurrent chromosomal translocations, almost 20% of childhood AML do not show any recognizable cytogenetic alteration and are defined as cytogenetically normal (CN)-AML. CN-AML patients have always showed a great variability in response to therapy and overall outcome, underlining the presence of unknown genetic changes, not detectable by conventional analyses, but relevant for pathogenesis, and outcome of AML. The development of novel genome-wide techniques such as next-generation sequencing, have tremendously improved our ability to interrogate the cancer genome. Based on this background, the aim of this research study was to investigate the mutational landscape of pediatric CN-AML patients negative for all the currently known somatic mutations reported in AML through whole-transcriptome sequencing (RNA-seq). RNA-seq performed on diagnostic leukemic blasts from 19 pediatric CN-AML cases revealed a considerable incidence of cryptic chromosomal rearrangements, with the identification of 21 putative fusion genes. Several of the fusion genes that were identified in this study are recurrent and might have a prognostic and/or therapeutic relevance. A paradigm of that is the CBFA2T3-GLIS2 fusion, which has been demonstrated to be a common alteration in pediatric CN-AML, predicting poor outcome. Important findings have been also obtained in the identification of novel therapeutic targets. On one side, the identification of NUP98-JARID1A fusion suggests the use of disulfiram; on the other, here we describe alteration-activating tyrosine kinases, providing functional data supporting the use of tyrosine kinase inhibitors to specifically inhibit leukemia cells. This study provides new insights in the knowledge of genetic alterations underlying pediatric AML, defines novel prognostic markers and putative therapeutic targets, and prospectively ensures a correct risk stratification and risk-adapted therapy also for the “all-neg” AML subgroup.

Genetic characterization, population history and evolutionary medicine perspective in two native south american populations: Yanesha and Wichi

Two Amerindian populations from the Peruvian Amazon (Yanesha) and from rural lowlands of the Argentinean Gran Chaco (Wichi) were analyzed. They represent two case study of the South American genetic variability. The Yanesha represent a model of population isolated for long-time in the Amazon rainforest, characterized by environmental and altitudinal stratifications. The Wichi represent a model of population living in an area recently colonized by European populations (the Criollos are the population of the admixed descendents), whose aim is to depict the native ancestral gene pool and the degree of admixture, in relation to the very high prevalence of Chagas disease. The methods used for the genotyping are common, concerning the Y chromosome markers (male lineage) and the mitochondrial markers (maternal lineage). The determination of the phylogeographic diagnostic polymorphisms was carried out by the classical techniques of PCR, restriction enzymes, sequencing and specific mini-sequencing. New method for the detection of the protozoa Trypanosoma cruzi was developed by means of the nested PCR. The main results show patterns of genetic stratification in Yanesha forest communities, referable to different migrations at different times, estimated by Bayesian analyses. In particular Yanesha were considered as a population of transition between the Amazon basin and the Andean Cordillera, evaluating the potential migration routes and the separation of clusters of community in relation to different genetic bio-ancestry. As the Wichi, the gene pool analyzed appears clearly differentiated by the admixed sympatric Criollos, due to strict social practices (deeply analyzed with the support of cultural anthropological tools) that have preserved the native identity at a diachronic level. A pattern of distribution of the seropositivity in relation to the different phylogenetic lineages (the adaptation in evolutionary terms) does not appear, neither Amerindian nor European, but in relation to environmental and living conditions of the two distinct subpopulations.

Novel tools for conservation genetics in marine fish: population structure and evolution of the European hake (Merluccius merluccius) inferred by SNP variation and applications to traceability

The research presented in my PhD thesis is part of a wider European project, FishPopTrace, focused on traceability of fish populations and products. My work was aimed at developing and analyzing novel genetic tools for a widely distributed marine fish species, the European hake (Merluccius merluccius), in order to investigate population genetic structure and explore potential applications to traceability scenarios. A total of 395 SNPs (Single Nucleotide Polymorphisms) were discovered from a massive collection of Expressed Sequence Tags, obtained by high-throughput sequencing, and validated on 19 geographic samples from Atlantic and Mediterranean. Genome-scan approaches were applied to identify polymorphisms on genes potentially under divergent selection (outlier SNPs), showing higher genetic differentiation among populations respect to the average observed across loci. Comparative analysis on population structure were carried out on putative neutral and outlier loci at wide (Atlantic and Mediterranean samples) and regional (samples within each basin) spatial scales, to disentangle the effects of demographic and adaptive evolutionary forces on European hake populations genetic structure. Results demonstrated the potential of outlier loci to unveil fine scale genetic structure, possibly identifying locally adapted populations, despite the weak signal showed from putative neutral SNPs. The application of outlier SNPs within the framework of fishery resources management was also explored. A minimum panel of SNP markers showing maximum discriminatory power was selected and applied to a traceability scenario aiming at identifying the basin (and hence the stock) of origin, Atlantic or Mediterranean, of individual fish. This case study illustrates how molecular analytical technologies have operational potential in real-world contexts, and more specifically, potential to support fisheries control and enforcement and fish and fish product traceability.

Wide-scale population genomics of Atlantic bluefin tuna (Thunnus thynnus) inferred by novel high-throughput technology

My PhD project was focused on Atlantic bluefin tuna, Thunnus thynnus, a fishery resource overexploited in the last decades. For a better management of stocks, it was necessary to improve scientific knowledge of this species and to develop novel tools to avoid collapse of this important commercial resource. To do this, we used new high throughput sequencing technologies, as Next Generation Sequencing (NGS), and markers linked to expressed genes, as SNPs (Single Nucleotide Polymorphisms). In this work we applied a combined approach: transcriptomic resources were used to build cDNA libreries from mRNA isolated by muscle, and genomic resources allowed to create a reference backbone for this species lacking of reference genome. All cDNA reads, obtained from mRNA, were mapped against this genome and, employing several bioinformatics tools and different restricted parameters, we achieved a set of contigs to detect SNPs. Once a final panel of 384 SNPs was developed, following the selection criteria, it was genotyped in 960 individuals of Atlantic bluefin tuna, including all size/age classes, from larvae to adults, collected from the entire range of the species. The analysis of obtained data was aimed to evaluate the genetic diversity and the population structure of Thunnus thynnus. We detect a low but significant signal of genetic differentiation among spawning samples, that can suggest the presence of three genetically separate reproduction areas. The adult samples resulted instead genetically undifferentiated between them and from the spawning populations, indicating a presence of panmictic population of adult bluefin tuna in the Mediterranean Sea, without different meta populations.

Next-Generation Sequencing-Based Mutations Scanning Strategy of the BCR-ABL Kinase Domain in Patients with PhiladelPhia-Chromosome Positive Leukemias Treated with Tyrosine Kinase Inhibitors

In chronic myeloid leukemia and Philadelphia-positive acute lymphoblastic leukemia patients resistant to tyrosine kinase inhibitors (TKIs), BCR-ABL kinase domain mutation status is an essential component of the therapeutic decision algorithm. The recent development of Ultra-Deep Sequencing approach (UDS) has opened the way to a more accurate characterization of the mutant clones surviving TKIs conjugating assay sensitivity and throughput. We decided to set-up and validated an UDS-based for BCR-ABL KD mutation screening in order to i) resolve qualitatively and quantitatively the complexity and the clonal structure of mutated populations surviving TKIs, ii) study the dynamic of expansion of mutated clones in relation to TKIs therapy, iii) assess whether UDS may allow more sensitive detection of emerging clones, harboring critical 2GTKIs-resistant mutations predicting for an impending relapse, earlier than SS. UDS was performed on a Roche GS Junior instrument, according to an amplicon sequencing design and protocol set up and validated in the framework of the IRON-II (Interlaboratory Robustness of Next-Generation Sequencing) International consortium.Samples from CML and Ph+ ALL patients who had developed resistance to one or multiple TKIs and collected at regular time-points during treatment were selected for this study. Our results indicate the technical feasibility, accuracy and robustness of our UDS-based BCR-ABL KD mutation screening approach. UDS was found to provide a more accurate picture of BCR-ABL KD mutation status, both in terms of presence/absence of mutations and in terms of clonal complexity and showed that BCR-ABL KD mutations detected by SS are only the “tip of iceberg”. In addition UDS may reliably pick 2GTKIs-resistant mutations earlier than SS in a significantly greater proportion of patients.The enhanced sensitivity as well as the possibility to identify low level mutations point the UDS-based approach as an ideal alternative to conventional sequencing for BCR-ABL KD mutation screening in TKIs-resistant Ph+ leukemia patients

Analysis of the pig genome for the identification of genomic regions affecting production traits

The aim of this work was to identify markers associated with production traits in the pig genome using different approaches. We focused the attention on Italian Large White pig breed using Genome Wide Association Studies (GWAS) and applying a selective genotyping approach to increase the power of the analyses. Furthermore, we searched the pig genome using Next Generation Sequencing (NSG) Ion Torrent Technology to combine selective genotyping approach and deep sequencing for SNP discovery. Other two studies were carried on with a different approach. Allele frequency changes for SNPs affecting candidate genes and at Genome Wide level were analysed to identify selection signatures driven by selection program during the last 20 years. This approach confirmed that a great number of markers may affect production traits and that they are captured by the classical selection programs. GWAS revealed 123 significant or suggestively significant SNP associated with Back Fat Thickenss and 229 associated with Average Daily Gain. 16 Copy Number Variant Regions resulted more frequent in lean or fat pigs and showed that different copies of those region could have a limited impact on fat. These often appear to be involved in food intake and behavior, beside affecting genes involved in metabolic pathways and their expression. By combining NGS sequencing with selective genotyping approach, new variants where discovered and at least 54 are worth to be analysed in association studies. The study of groups of pigs undergone to stringent selection showed that allele frequency of some loci can drastically change if they are close to traits that are interesting for selection schemes. These approaches could be, in future, integrated in genomic selection plans.

Whole Exome Sequencing widens the genetic landscape of Hereditary Optic Neuropathies

Hereditary optic neuropathies (HON) are a genetic cause of visual impairment characterized by degeneration of retinal ganglion cells. The majority of HON are caused by pathogenic variants in mtDNA genes and in gene OPA1. However, several other genes can cause optic atrophy and can only be identified by high throughput genetic analysis. Whole Exome Sequencing (WES) is becoming the primary choice in rare disease molecular diagnosis, being both cost effective and informative. We performed WES on a cohort of 106 cases, of which 74 isolated ON patients (ON) and 32 syndromic ON patients (sON). The total diagnostic yield amounts to 27%, slightly higher for syndromic ON (31%) than for isolated ON (26%). The majority of genes found are related to mitochondrial function and already reported for harbouring HON pathogenic variants: ACO2, AFG3L2, C19orf12, DNAJC30, FDXR, MECR, MTFMT, NDUFAF2, NDUFB11, NDUFV2, OPA1, PDSS1, SDHA, SSBP1, and WFS1. Among these OPA1, ACO2, and WFS1 were confirmed as the most relevant genetic causes of ON. Moreover, several genes were identified, especially in sON patients, with direct impairment of non-mitochondrial molecular pathways: from autophagy and ubiquitin system (LYST, SNF8, WDR45, UCHL1), to neural cells development and function (KIF1A, GFAP, EPHB2, CACNA1A, CACNA1F), but also vitamin metabolism (SLC52A2, BTD), cilia structure (USH2A), and nuclear pore shuttling (NUTF2). Functional validation on yeast model was performed for pathogenic variants detected in MECR, MTFMT, SDHA, and UCHL1 genes. For SDHA and UCHL1 also muscle biopsy and fibroblast cell lines from patients were analysed, pointing to possible pathogenic mechanisms that will be investigated in further studies. In conclusion, WES proved to be an efficient tool when applied to our ON cohort, for both common disease-genes identification and novel genes discovery. It is therefore recommended to consider WES in ON molecular diagnostic pipeline, as for other rare genetic diseases.

Artisanal food productions of animal origin: exploring food safety in the age of Whole Genome Sequencing

The artisanal food chain is enriched by a wide diversity of local food productions with delightful organoleptic characteristics and valuable nutritional properties. Despite their increasing worldwide popularity and appeal, several food safety challenges are addressed in artisanal facilities context suffering from less standardized processing conditions. In such scenario, recent advances in molecular typing and genomic surveillance (e.g., Whole Genome Sequencing [WGS]) represent an unprecedent solution capable of inferring sources of contamination as well as contributing to food safety along the artisanal food continuum. The overall objective of this PhD thesis was to explore potential microbial hazards among different artisanal food productions of animal origins (dairy and meat-derived) typical of the food culture and heritage landscape belonging to Mediterranean countries. Three different studies were then carried out, specifically focussing on: 1) compare the seasonal variability of microbiological quality and potential occurrence of microbial hazards in two batches of Italian artisanal fermented dairy and meat productions; 2) Investigate genetic relationships as well as virulome and resistome of foodborne pathogens isolated within dairy and meat-derived productions located in Italy, Spain, Portugal and Morocco; 3) investigate the population structure, virulome, resistome and mobilome of Klebsiella spp. isolates collected from study 1, including an extended range of public sequences.

Characterization of the molecular mechanisms mediated by the insulin-like growth factor-2 mRNA-binding proteins and Ephrin receptor A2 (EPHA2) signaling in the malignancy of CIC-DUX4 sarcomas

Ewing sarcoma (EWS) and CIC-DUX4 sarcoma (CDS) are pediatric fusion gene-driven tumors of mesenchymal origin characterized by an extremely stable genome and limited clinical solutions. Post-transcriptional regulatory mechanisms are crucial for understanding the development of this class of tumors. RNA binding proteins (RBPs) play a crucial role in the aggressiveness of these tumors. Numerous RBP families are dysregulated in cancer, including IGF2BPs. Among these, IGF2BP3 is a negative prognostic factor in EWS because it promotes cell growth, chemoresistence, and induces the metastatic process. Based on preliminary RNA sequencing data from clinical samples of EWS vs CDS patients, three major axes that are more expressed in CDS have been identified, two of which are dissected in this PhD work. The first involves the transcription factor HMGA2, IGF2BP2-3, and IGF2; the other involves the ephrin receptor system, particularly EphA2. EphA2 is involved in numerous cellular functions during embryonic stages, and its increased expression in adult tissues is often associated with pathological conditions. In tumors, its role is controversial because it can be associated with both pro- and anti-tumoral mechanisms. In EWS, it has been shown to play a role in promoting cell migration and neoangiogenesis. Our study has confirmed that the HMGA2/IGF2BPs/IGF2 axis contributes to CDS malignancy, and Akt hyperactivation has a strong impact on migration. Using loss/gain of function models for EphA2, we confirmed that it is a substrate of Akt, and Akt hyperactivation in CDS triggers ligand-independent activation of EphA2 through phosphorylation of S897. Moreover, the combination of Trabectedin and NVP/BEZ235 partially inhibits Akt/mTOR activation, resulting in reduced tumor growth in vivo. Inhibition of EphA2 through ALWII 41_27 significantly reduces migration in vitro. The project aim is the identification of target molecules in CDS that can distinguish it from EWS and thus develop new targeted therapeutic strategies.

The RAS-Lung project: implementing blood and tissue genotyping in KRAS-Positive non-small cell lung cancer treatment

Background: The frontline management of non-oncogene addicted non-small cell lung cancer (NSCLC) involves immunotherapy (ICI) alone or combined with chemotherapy (CT-ICI). As therapeutic options expand, refining NSCLC genotyping gains paramount importance. The dynamic landscape of KRAS-positive NSCLC presents a spectrum of treatment options, including ICI, targeted therapy, and combination strategies currently under investigation. Methods: The two-year RASLUNG project, featuring both retrospective and prospective cohorts, aimed to analyze the predictive and prognostic impact of KRAS mutations on tumor tissue and circulating DNA (ctDNA). Secondary objectives included assessing the roles of co-mutations and longitudinal changes in KRAS mutant copies concerning treatment response and survival outcomes. An external validation study confirmed the prognostic or predictive significance of co-mutations. Results: In the prospective cohort (n=24), patients with liver metastases exhibited significantly elevated ctDNA levels(p=0.01), while those with >3 metastatic sites showed increased Allele Frequency (AF) (P=0.002). Median overall survival (OS) was 7.5 months, progression-free survival (PFS) was 4.0 months, and the objective response rate (ORR) was 33.3%. Higher AF correlated with an increased risk of death (HR 1.04, p = 0.03), though not progression. Notably, a reduction in plasma DNA levels was significantly associated with objective response(p=0.01). In the retrospective cohort, KRAS and STK11 mutations co-occurred in 14/21 patients (p=0.053). STK11 mutations were independently detrimental to OS (HR 1.97, p=0.025) after adjusting for various factors. KRAS tissue AF did not correlate with OS or PFS. Within the validation dataset, STK11 mutations were significantly associated with an increased risk of death in univariate (HR 2.01, p<0.001) and multivariate models (HR 1.66, p=0.001) after adjustments. Conclusion: The RAS-Lung Project, employing innovative genotyping techniques, underscores the significance of comprehensive NSCLC genotyping. Tailored next-generation sequencing (NGS) and ctDNA monitoring may offer potential benefits in navigating the evolving landscape of KRAS-positive NSCLC treatment.