Two new HLA-DQB1*03:02:01 variants with substitutions in intron 2: HLA-DQB1*03:02:01:19 and -DQB1*03:02:01:21.

Two different single nucleotide substitutions in intron 2 give rise to novel HLA-DQB1*03:02:01 alleles.

Organellar-genome analyses from the lycophyte genus Isoetes L. show one of the highest frequencies of RNA editing in land plants.

RNA editing is a post-transcriptional process that challenges the central dogma of molecular biology by modifying RNA sequences, introducing nucleotide changes at specific sites, and generating functional diversity beyond the genomic code, especially when it concerns organellar transcripts. In plants, this phenomenon is widespread, but its extent varies significantly among species and organellar genomes. Among land plants, the heterosporous lycophytes (i.e., Isoetes and Selaginella) stand out for their exceptionally high numbers of RNA-editing sites, despite their morphological stasis and ancient lineage. In this study, we explore the complete set of organellar protein-coding genes in the aquatic plant group Isoetes, providing a detailed analysis of RNA editing in both the mitochondrial and plastid genomes. Our findings reveal a remarkable abundance of RNA editing, particularly in the mitochondrial genome, with thousands of editing sites identified. Interestingly, the majority of these edits result in non-silent substitutions, suggesting a role in fine-tuning protein structure and function. Furthermore, we observe a consistent trend of increased hydrophobicity in membrane-bound proteins, supporting the notion that RNA editing may confer a selective advantage by preserving gene functionality in Isoetes. The conservation of highly edited RNA sequences over millions of years underscores the evolutionary significance of RNA editing. Additionally, the study sheds light on the dynamic nature of RNA editing, with shared editing sites reflecting common ancestry whereas exclusive edits matching more recent radiation events within the genus. This work advances our understanding of the intricate interplay between RNA editing, adaptation, and evolution in land plants and highlights the unique genomic features of Isoetes, providing a foundation for further investigations into the functional consequences of RNA editing in this enigmatic plant lineage.

In silico prospecting of the mtDNA of Macrobrachium amazonicum from transcriptome data.

BACKGROUND: Macrobrachium amazonicum is a freshwater prawn widely distributed in South America that is undergoing speciation, so the denomination "M. amazonicum complex" is used for it. The mitochondrial cytochrome c oxidase subunit I (COI) gene has been used to elucidate this speciation, but heteroplasmies and pseudogenes have been recorded, making separation difficult. Obtaining genes from cDNA (RNA) rather than genomic DNA is an effective tool to mitigate those two types of occurrences. The aim of this study was to assemble in silico the mitochondrial DNA (mtDNA) of the Amazonian coastal population of M. amazonicum inhabiting the state of Pará. RESULTS: Sequences were obtained from the prawn's transcriptome using the de novo approach. Six libraries of cDNA from the androgen gland, hepatopancreas, and muscle tissue were used. The mtDNA of M. amazonicum was 14,960 bp in length. It contained 13 protein-coding genes, 21 complete transfer RNAs, and the 12S and 16S subunits of ribosomal RNA. All regions were found on the light strand except tRNAGln, which was on the heavy strand. The control region (D-loop) was not recovered, making for a gap of 793 bp. The cladogram showed the formation of the well-defined Macrobrachium clade, with high support value in the established branches (91-100). The three-dimensional spatial conformation of the mtDNA-encoded proteins showed that most of them were mainly composed of major α-helices that typically shows in those proteins inserted in the membrane (mitochondrial). CONCLUSIONS: It was possible to assemble a large part of the mitochondrial genome of M. amazonicum in silico using data from other genomes deposited in GenBank and to validate it through the similarities between its COI and 16S genes and those from animals of the same region deposited in GenBank. Depositing the M. amazonicum mtDNA sequences in GenBank may help solve the taxonomic problems recorded for the species, in addition to providing complete sequences of candidate coding genes for use as biomarkers in ecological studies.

Solving the unsolved genetic epilepsies: Current and future perspectives.

Many patients with epilepsy undergo exome or genome sequencing as part of a diagnostic workup; however, many remain genetically unsolved. There are various factors that account for negative results in exome/genome sequencing for patients with epilepsy: (1) the underlying cause is not genetic; (2) there is a complex polygenic explanation; (3) the illness is monogenic but the causative gene remains to be linked to a human disorder; (4) family segregation with reduced penetrance; (5) somatic mosaicism or the complexity of, for example, a structural rearrangement; or (6) limited knowledge or diagnostic tools that hinder the proper classification of a variant, resulting in its designation as a variant of unknown significance. The objective of this review is to outline some of the diagnostic options that lie beyond the exome/genome, and that might become clinically relevant within the foreseeable future. These options include: (1) re-analysis of older exome/genome data as knowledge increases or symptoms change; (2) looking for somatic mosaicism or long-read sequencing to detect low-complexity repeat variants or specific structural variants missed by traditional exome/genome sequencing; (3) exploration of the non-coding genome including disruption of topologically associated domains, long range non-coding RNA, or other regulatory elements; and finally (4) transcriptomics, DNA methylation signatures, and metabolomics as complementary diagnostic methods that may be used in the assessment of variants of unknown significance. Some of these tools are currently not integrated into standard diagnostic workup. However, it is reasonable to expect that they will become increasingly available and improve current diagnostic capabilities, thereby enabling precision diagnosis in patients who are currently undiagnosed.

Comparative mitogenomic and evolutionary analysis of Lycaenidae (Insecta: Lepidoptera): Potential association with high-altitude adaptation.

Harsh environments (e.g., hypoxia and cold temperatures) of the Qinghai-Tibetan Plateau have a substantial influence on adaptive evolution in various species. Some species in Lycaenidae, a large and widely distributed family of butterflies, are adapted to the Qinghai-Tibetan Plateau. Here, we sequenced four mitogenomes of two lycaenid species in the Qinghai-Tibetan Plateau and performed a detailed comparative mitogenomic analysis including nine other lycaenid mitogenomes (nine species) to explore the molecular basis of high-altitude adaptation. Based on mitogenomic data, Bayesian inference, and maximum likelihood methods, we recovered a lycaenid phylogeny of [Curetinae + (Aphnaeinae + (Lycaeninae + (Theclinae + Polyommatinae)))]. The gene content, gene arrangement, base composition, codon usage, and transfer RNA genes (sequence and structure) were highly conserved within Lycaenidae. TrnS1 not only lacked the dihydrouridine arm but also showed anticodon and copy number diversity. The ratios of non-synonymous substitutions to synonymous substitutions of 13 protein-coding genes (PCGs) were less than 1.0, indicating that all PCGs evolved under purifying selection. However, signals of positive selection were detected in cox1 in the two Qinghai-Tibetan Plateau lycaenid species, indicating that this gene may be associated with high-altitude adaptation. Three large non-coding regions, i.e., rrnS-trnM (control region), trnQ-nad2, and trnS2-nad1, were found in the mitogenomes of all lycaenid species. Conserved motifs in three non-coding regions (trnE-trnF, trnS1-trnE, and trnP-nad6) and long sequences in two non-coding regions (nad6-cob and cob-trnS2) were detected in the Qinghai-Tibetan Plateau lycaenid species, suggesting that these non-coding regions were involved in high-altitude adaptation. In addition to the characterization of Lycaenidae mitogenomes, this study highlights the importance of both PCGs and non-coding regions in high-altitude adaptation.

Mapping the Constrained Coding Regions in the Human Genome to Their Corresponding Proteins.

Constrained Coding Regions (CCRs) in the human genome have been derived from DNA sequencing data of large cohorts of healthy control populations, available in the Genome Aggregation Database (gnomAD) [1]. They identify regions depleted of protein-changing variants and thus identify segments of the genome that have been constrained during human evolution. By mapping these DNA-defined regions from genomic coordinates onto the corresponding protein positions and combining this information with protein annotations, we have explored the distribution of CCRs and compared their co-occurrence with different protein functional features, previously annotated at the amino acid level in public databases. As expected, our results reveal that functional amino acids involved in interactions with DNA/RNA, protein-protein contacts and catalytic sites are the protein features most likely to be highly constrained for variation in the control population. More surprisingly, we also found that linear motifs, linear interacting peptides (LIPs), disorder-order transitions upon binding with other protein partners and liquid-liquid phase separating (LLPS) regions are also strongly associated with high constraint for variability. We also compared intra-species constraints in the human CCRs with inter-species conservation and functional residues to explore how such CCRs may contribute to the analysis of protein variants. As has been previously observed, CCRs are only weakly correlated with conservation, suggesting that intraspecies constraints complement interspecies conservation and can provide more information to interpret variant effects.

Looking beyond the cytogenetics in haematological malignancies: decoding the role of tandem repeats in DNA repair genes.

BACKGROUND: In cancer research, one of the most significant findings was to characterize the DNA repair deficiency as carcinogenic. Amongst the various repair mechanisms, mismatch repair (MMR) and direct reversal of DNA damage systems are designated as multilevel safeguards in the human genome. Defects in these elevate the rate of mutations and results in dire consequences like cancer. Of the several molecular signatures in human genome, tandem repeats (TRs) appear at various frequencies in the exonic, intronic, and regulatory regions of the DNA. Hypervariability among these repeats in the coding and non-coding regions of the genes is well characterized for solid tumours, but its significance in haematologic malignancies remains to be explored. The purpose of our study was to elucidate the role of nucleotide repeat instability in the coding and non-coding regions of 10 different repair genes in myeloid and lymphoid cell lines compared to the control samples. METHODS AND RESULTS: We selected MMR deficient extensively studied microsatellite instable colorectal cancer (HCT116), and MMR proficient breast cancer (MCF-7) cells along with underemphasized haematologic cancer cell lines to decipher the hypermutability of tandem repeats. A statistically significant TR variation was observed for MSH2 and MSH6 genes in 4 and 3 of the 6 cell lines respectively. KG1 (AML) and Daudi (Burkitt's lymphoma) were found to have compromised DNA repair competency with highly unstable nucleotide repeats. CONCLUSION: Taken together, the results suggest that mutable TRs in intronic and non-intronic regions of repair genes in blood cancer might have a tumorigenic role. Since this is a pilot study on cell lines, high throughput research in large cohorts can be undertaken to reveal novel diagnostic markers for unexplained blood cancer patients with normal karyotypes or otherwise with karyotypic defects.

Identification of Seven Additional Genome Segments of Grapevine-Associated Jivivirus 1.

Jiviruses are a group of recently described viruses characterized with a tripartite genome and having affinities with Virgaviridae (RNA1 and 2) and Flaviviridae (RNA3). Using a combination of high-throughput sequencing, datamining and RT-PCR approaches, we demonstrate here that in grapevine samples infected by grapevine-associated jivivirus 1 (GaJV-1) up to 7 additional molecules can be consistently detected with conserved 5' and 3' non-coding regions in common with the three previously identified GaJV-1 genomic RNAs. RNA4, RNA5, RNA6, RNA7, RNA8 and RNA10, together with a recombinant RNArec7-8, are all members of a family sharing a previously non recognized conserved protein domain, while RNA9 is part of a distinct family characterized by another conserved motif. Datamining of pecan (Carya illinoinensis) public transcriptomic data allowed the identification of two further jiviviruses and the identification of supplementary genomic RNAs with homologies to those of GaJV-1. Taken together, these results reshape our vision of the divided genome of jiviviruses and raise novel questions about the function(s) of the proteins encoded by jiviviruses supplementary RNAs.

Divergent variant patterns among 19 patients with Rubinstein-Taybi syndrome uncovered by comprehensive genetic analysis including whole genome sequencing.

Rubinstein-Taybi syndrome (RSTS) is characterized by dysmorphic facial features, broad thumbs, and intellectual disability. CREB-binding protein (CREBBP) or E1A-binding protein P300 (EP300) are causative genes. To elucidate the underlying genetic and genomic architecture related to the RSTS phenotype, we performed comprehensive genetic analysis targeting CREBBP and/or EP300 in 22 clinically diagnosed patients. During the 11-year study period, we used several analysis methods including high-resolution melting, array-based comparative genomic hybridization, panel-based exome sequencing, whole exome sequencing, and whole genome sequencing (WGS). We identified the causative variants in 19 patients (86.3%), but they were variable and complex, so we must combine multiple analysis methods. Notably, we found genetic alterations in the non-coding regions of two patients (10.5%, 2/19): scattered deletions including a partial 5'-untranslated region of CREBBP in one patient (all coding exons were intact), and a deep 229-bp intronic deletion in another patient, resulting in a splicing error. Furthermore, we identified rare clinical findings: two patients with an EP300 variant showed abnormal development of the neural tube, and one patient with a CREBBP variant had anorectal atresia with a cloaca. Our findings expand the allelic heterogeneity of RSTS, underscore the utility of comprehensive genetic analysis, and suggest that WGS may be a practical diagnostic strategy.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa.

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome's size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

Heritability and De Novo Mutations in Oesophageal Atresia and Tracheoesophageal Fistula Aetiology.

Tracheoesophageal Fistula (TOF) is a congenital anomaly for which the cause is unknown in the majority of patients. OA/TOF is a variable feature in many (often mono-) genetic syndromes. Research using animal models targeting genes involved in candidate pathways often result in tracheoesophageal phenotypes. However, there is limited overlap in the genes implicated by animal models and those found in OA/TOF-related syndromic anomalies. Knowledge on affected pathways in animal models is accumulating, but our understanding on these pathways in patients lags behind. If an affected pathway is associated with both animals and patients, the mechanisms linking the genetic mutation, affected cell types or cellular defect, and the phenotype are often not well understood. The locus heterogeneity and the uncertainty of the exact heritability of OA/TOF results in a relative low diagnostic yield. OA/TOF is a sporadic finding with a low familial recurrence rate. As parents are usually unaffected, de novo dominant mutations seems to be a plausible explanation. The survival rates of patients born with OA/TOF have increased substantially and these patients start families; thus, the detection and a proper interpretation of these dominant inherited pathogenic variants are of great importance for these patients and for our understanding of OA/TOF aetiology.

Episodes of Rapid Recovery of the Functional Activity of the ras85D Gene in the Evolutionary History of Phylogenetically Distant Drosophila Species.

As assemblies of genomes of new species with varying degrees of relationship appear, it becomes obvious that structural rearrangements of the genome, such as inversions, translocations, and transposon movements, are an essential and often the main source of evolutionary variation. In this regard, the following questions arise. How conserved are the regulatory regions of genes? Do they have a common evolutionary origin? And how and at what rate is the functional activity of genes restored during structural changes in the promoter region? In this article, we analyze the evolutionary history of the formation of the regulatory region of the ras85D gene in different lineages of the genus Drosophila, as well as the participation of mobile elements in structural rearrangements and in the replacement of specific areas of the promoter region with those of independent evolutionary origin. In the process, we substantiate hypotheses about the selection of promoter elements from a number of frequently repeated motifs with different degrees of degeneracy in the ancestral sequence, as well as about the restoration of the minimum required set of regulatory sequences using a conversion mechanism or similar.

Rules are made to be broken: a "simple" model organism reveals the complexity of gene regulation.

Global methods for assaying translation have greatly improved our understanding of the protein-coding capacity of the genome. In particular, it is now possible to perform genome-wide and condition-specific identification of translation initiation sites through modified ribosome profiling methods that selectively capture initiating ribosomes. Here we discuss our recent study applying such an approach to meiotic and mitotic timepoints in the simple eukaryote, budding yeast, as an example of the surprising diversity of protein products-many of which are non-canonical-that can be revealed by such methods. We also highlight several key challenges in studying non-canonical protein isoforms that have precluded their prior systematic discovery. A growing body of work supports expanded use of empirical protein-coding region identification, which can help relieve some of the limitations and biases inherent to traditional genome annotation approaches. Our study also argues for the adoption of less static views of gene identity and a broader framework for considering the translational capacity of the genome.

Epigenetic Targeting of Mcl-1 Is Synthetically Lethal with Bcl-xL/Bcl-2 Inhibition in Model Systems of Glioblastoma.

Apoptotic resistance remains a hallmark of glioblastoma (GBM), the most common primary brain tumor in adults, and a better understanding of this process may result in more efficient treatments. By utilizing chromatin immunoprecipitation with next-generation sequencing (CHIP-seq), we discovered that GBMs harbor a super enhancer around the Mcl-1 locus, a gene that has been known to confer cell death resistance in GBM. We utilized THZ1, a known super-enhancer blocker, and BH3-mimetics, including ABT263, WEHI-539, and ABT199. Combined treatment with BH3-mimetics and THZ1 led to synergistic growth reduction in GBM models. Reduction in cellular viability was accompanied by significant cell death induction with features of apoptosis, including disruption of mitochondrial membrane potential followed by activation of caspases. Mechanistically, THZ1 elicited a profound disruption of the Mcl-1 enhancer region, leading to a sustained suppression of Mcl-1 transcript and protein levels, respectively. Mechanism experiments suggest involvement of Mcl-1 in the cell death elicited by the combination treatment. Finally, the combination treatment of ABT263 and THZ1 resulted in enhanced growth reduction of tumors without induction of detectable toxicity in two patient-derived xenograft models of GBM in vivo. Taken together, these findings suggest that combined epigenetic targeting of Mcl-1 along with Bcl-2/Bcl-xL is potentially therapeutically feasible.

Identification of High Molecular Variation Loci in Complete Chloroplast Genomes of Mammillaria (Cactaceae, Caryophyllales).

In plants, partial DNA sequences of chloroplasts have been widely used in evolutionary studies. However, the Cactaceae family (1500-1800 species) lacks molecular markers that allow a phylogenetic resolution between species and genera. In order to identify sequences with high variation levels, we compared previously reported complete chloroplast genomes of seven species of Mammillaria. We identified repeated sequences (RSs) and two types of DNA variation: short sequence repeats (SSRs) and divergent homologous loci. The species with the highest number of RSs was M. solisioides (256), whereas M. pectinifera contained the highest amount of SSRs (84). In contrast, M. zephyranthoides contained the lowest number (35) of both RSs and SSRs. In addition, five of the SSRs were found in the seven species, but only three of them showed variation. A total of 180 homologous loci were identified among the seven species. Out of these, 20 loci showed a molecular variation of 5% to 31%, and 12 had a length within the range of 150 to 1000 bp. We conclude that the high levels of variation at the reported loci represent valuable knowledge that may help to resolve phylogenetic relationships and that may potentially be convenient as molecular markers for population genetics and phylogeographic studies.

Codon usage and modular interactions between messenger RNA coding regions and small RNAs in Escherichia coli.

BACKGROUND: Small RNAs (sRNAs) are key regulators of gene expression in bacteria. In addition to modulating translation initiation, sRNAs can interact with mRNA coding regions to regulate mRNA stability and translation efficiency, enhancing or impeding progression of the ribosome along the mRNA. Since most amino acids are decoded by more than one codon (synonymous) we asked as to whether there is a codon bias in the interaction of sRNAs with coding regions of mRNAs. Therefore, we explored whether there are differences in codon usage or tRNA availability according to whether an mRNA is regulated by sRNAs or not. We also explored these parameters in the coding interaction regions in mRNAs. We focused our analysis on sRNAs that regulate multiple mRNAs. RESULTS: We found differences in codon adaptation index and tRNA adaptation index between sRNA-regulated and non-sRNA-regulated mRNAs. Interestingly, the sRNA-mRNA interacting regions tended to be enriched in unpreferred codons decoded by scarce tRNAs. We also found that sRNAs with multiple targets often contained modular segments capable of recognizing conserved motifs among these mRNAs. CONCLUSIONS: Our results show that sRNAs in E. coli tend to recognize mRNA coding regions in which the ribosome is predicted to advance at low speeds. Identified motifs in interacting regions are conserved among mRNAs that are recognized by the same sRNA.

Genomic non-coding regions reveal hidden patterns of mumps virus circulation in Spain, 2005 to 2015.

BackgroundSince mumps vaccination was introduced in 1981 in Spain, the incidence of the disease has dropped significantly. However, cyclic epidemic waves and outbreaks still occur, despite high vaccination coverage. The World Health Organization (WHO) recommends genotyping to trace the pattern of mumps virus (MuV) circulation. Genotype H was predominant in Spain, but was replaced in 2005 by genotype G which has subsequently remained dominant. Of the small hydrophobic protein gene sequences, 78% are identical and belong to the MuVi/ Sheffield.GBR.1.05/[G]-variant. Aim: Our study aimed to investigate whether the circulation of MuV strains in Spain was continuous after the emergence of genotype G in 2005. Method: We obtained 46 samples from Spanish patients infected with MuVi/Sheffield.GBR.1.05/[G] during two epidemic waves and analysed them using new molecular markers based on genomic non-coding regions (NCRs) that discriminate subvariants of this virus strain. Results: Phylogenetic analyses of the nucleoprotein-phosphoprotein and matrix protein-fusion protein NCR indicated strain replacement after a drop in incidence in 2009, which had not been detectable by SH sequencing. Clustering of sequences from patients epidemiologically linked in the same outbreak suggests a potential use for these NCRs in outbreak characterisation. Conclusion: We suggest to consider their use in conjunction with the SH gene in the future WHO recommendations for MuV epidemiological surveillance.

Functional assessment of human enhancer activities using whole-genome STARR-sequencing.

BACKGROUND: Genome-wide quantification of enhancer activity in the human genome has proven to be a challenging problem. Recent efforts have led to the development of powerful tools for enhancer quantification. However, because of genome size and complexity, these tools have yet to be applied to the whole human genome. RESULTS:  In the current study, we use a human prostate cancer cell line, LNCaP as a model to perform whole human genome STARR-seq (WHG-STARR-seq) to reliably obtain an assessment of enhancer activity. This approach builds upon previously developed STARR-seq in the fly genome and CapSTARR-seq techniques in targeted human genomic regions. With an improved library preparation strategy, our approach greatly increases the library complexity per unit of starting material, which makes it feasible and cost-effective to explore the landscape of regulatory activity in the much larger human genome. In addition to our ability to identify active, accessible enhancers located in open chromatin regions, we can also detect sequences with the potential for enhancer activity that are located in inaccessible, closed chromatin regions. When treated with the histone deacetylase inhibitor, Trichostatin A, genes nearby this latter class of enhancers are up-regulated, demonstrating the potential for endogenous functionality of these regulatory elements. CONCLUSION: WHG-STARR-seq provides an improved approach to current pipelines for analysis of high complexity genomes to gain a better understanding of the intricacies of transcriptional regulation.

Significant differences in terms of codon usage bias between bacteriophage early and late genes: a comparative genomics analysis.

BACKGROUND: Viruses undergo extensive evolutionary selection for efficient replication which effects, among others, their codon distribution. In the current study, we aimed at understanding the way evolution shapes the codon distribution in early vs. late viral genes in terms of their expression during different stages in the viral replication cycle. To this end we analyzed 14 bacteriophages and 11 human viruses with available information about the expression phases of their genes. RESULTS: We demonstrated evidence of selection for distinct composition of synonymous codons in early and late viral genes in 50% of the analyzed bacteriophages. Among others, this phenomenon may be related to the time specific adaptation of the viral genes to the translation efficiency factors involved at different bacteriophage developmental stages. Specifically, we showed that the differences in codon composition in different temporal gene groups cannot be explained only by phylogenetic proximities between the analyzed bacteriophages, and can be partially explained by differences in the adaptation to the host tRNA pool, nucleotide bias, GC content and more. In contrast, no difference in temporal regulation of synonymous codon usage was observed in human viruses, possibly because of a stronger selection pressure due to a larger effective population size in bacteriophages and their bacterial hosts. CONCLUSIONS: The codon distribution in large fractions of bacteriophage genomes tend to be different in early and late genes. This phenomenon seems to be related to various aspects of the viral life cycle, and to various intracellular processes. We believe that the reported results should contribute towards better understanding of viral evolution and may promote the development of relevant procedures in synthetic virology.

Integral Phylogenomic Approach over Ilex L. Species from Southern South America.

The use of molecular markers with inadequate variation levels has resulted in poorly resolved phylogenetic relationships within Ilex. Focusing on southern South American and Asian species, we aimed at contributing informative plastid markers. Also, we intended to gain insights into the nature of morphological and physiological characters used to identify species. We obtained the chloroplast genomes of I.paraguariensis and I. dumosa, and combined these with all the congeneric plastomes currently available to accomplish interspecific comparisons and multilocus analyses. We selected seven introns and nine IGSs as variable non-coding markers that were used in phylogenomic analyses. Eight extra IGSs were proposed as candidate markers. Southern South American species formed one lineage, except for I. paraguariensis, I. dumosa and I. argentina, which occupied intermediate positions among sampled taxa; Euroasiatic species formed two lineages. Some concordant relationships were retrieved from nuclear sequence data. We also conducted integral analyses, involving a supernetwork of molecular data, and a simultaneous analysis of quantitative and qualitative morphological and phytochemical characters, together with molecular data. The total evidence tree was used to study the evolution of non-molecular data, evidencing fifteen non-ambiguous synapomorphic character states and consolidating the relationships among southern South American species. More South American representatives should be incorporated to elucidate their origin.