Nucleotide Repeats Research Articles

The molecular biology of FMRP: new insights into fragile X syndrome.

Fragile X mental retardation protein (FMRP) is the product of the fragile X mental retardation 1 gene (FMR1), a gene that - when epigenetically inactivated by a triplet nucleotide repeat expansion - causes the neurodevelopmental disorder fragile X syndrome (FXS). FMRP is a widely expressed RNA-binding protein with activity that is essential for proper synaptic plasticity and architecture, aspects of neural function that are known to go awry in FXS. Although the neurophysiology of FXS has been described in remarkable detail, research focusing on the molecular biology of FMRP has only scratched the surface. For more than two decades, FMRP has been well established as a translational repressor; however, recent whole transcriptome and translatome analyses in mouse and human models of FXS have shown that FMRP is involved in the regulation of nearly all aspects of gene expression. The emerging mechanistic details of the mechanisms by which FMRP regulates gene expression may offer ways to design new therapies for FXS.

Mapping Internally and Environmentally Stimulated RNA Foci Formation

Formation of aberrant membraneless cellular compartments, RNA foci, in the cell nucleus is a common pathological feature of nucleotide repeat expansion disorders including Huntington's disease, myotonic dystrophy, and amyotrophic lateral sclerosis (ALS). Nucleotide repeat expansions are defined by consecutive repeats of short sequences in the coding or non-coding region of an RNA transcript. Above a critical length these repeat RNAs form RNA foci which mislocalize to nuclear speckles (another membraneless compartment) instead of being exported to the cytoplasm for decay.

Development of EST-Molecular Markers from RNA Sequencing for Genetic Management and Identification of Growth Traits in Potato Grouper (Epinephelus tukula)

Simple SummaryThe potato grouper is a novel aquaculture species in Taiwan. Due to the lack of genetic information concerning this species, we have developed molecular markers based on transcriptome sequencing and further characterized their association with gene diversity and growth traits of this species. Ultimately, these markers could be utilized as accurate and efficient tools for genetic management and marker-assisted selection of potato grouper with distinct growth traits.The accuracy and efficiency of marker-assisted selection (MAS) has been proven for economically critical aquaculture species. The potato grouper (Epinephelus tukula), a novel cultured grouper species in Taiwan, shows large potential in aquaculture because of its fast growth rate among other groupers. Because of the lack of genetic information for the potato grouper, the first transcriptome and expressed sequence tag (EST)-derived simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were developed. Initially, the transcriptome was obtained from seven cDNA libraries by using the Illumina platform. De novo transcriptome of the potato grouper yielded 51.34 Gb and 111,490 unigenes. The EST-derived SSR and SNP markers were applied in genetic management, in parentage analysis, and to discover the functional markers of economic traits. The F1 juveniles were identified as siblings from one pair of parents (80 broodstocks). Fast- and slow-growth individuals were analyzed using functional molecular markers and through their association with growth performance. The results revealed that two SNPs were correlated with growth traits. The transcriptome database obtained in this study and its derived SSR and SNP markers may be applied not only for MAS but also to maintain functional gene diversity in the novel cultured grouper.

Label-free non-invasive classification of rice seeds using optical coherence tomography assisted with deep neural network

Identification of the seed varieties is essential in the quality control and high yield crop growth. The existing methods of varietal identification rely primarily on visual examination and DNA fingerprinting. Although the pattern of DNA fingerprinting allows precise classification of seed varieties but fraught with challenges such as low rate of polymorphism amongst closely related species, destructive method of analysis and a huge cost involved in identification of robust markers such as simple sequence repeat (SSR) and single nucleotide polymorphisms. Here, we propose a fast, non-contact and non-invasive technique, deep learning assisted optical coherence tomography (OCT) for subsurface imaging in order to distinguish different seed varieties. The volumetric dataset of, (a) four rice varieties (PUSA Basmati 1, PUSA 1509, PUSA 44 and IR 64) and, (b) seven morphologically similar seeds of rice landrace Pokkali was acquired using OCT technique. A feedforward deep neural network is implemented for deep feature extraction and to classify the OCT images into their relevant classes. The proposed method provides the classification accuracy of 89.6% for the dataset of total 158,421 OCT images and 82.5% in classifying the dataset of total 56,301 OCT images collected from Pokkali seeds. The current technique can accurately classify seed varieties irrespective of the morphological similarities and can be adopted for the removal of varietal duplication and assessment of the purity of the seeds.

Drug-resistant Mycobacterium tuberculosis and its genotypes isolated from an outbreak in western Thailand.

Multidrug-resistant TB (MDR-TB) outbreaks have occurred in the Thamaka district, Kanchanaburi province in Thailand. Seventy-two isolates, which included 7% mono-, 30.6% MDR and extensively drug-resistant TB (XDR-TB), were genotyped by spoligotyping, mycobacterial interspersed repetitive unit-variable-number tandem repeat (MIRU-VNTR) and single nucleotide polymorphism genotyping, and their drug resistance was analysed. The spoligotyping results showed that Beijing spoligo-international type (SIT)1 was predominant (n=38; 52.8%) while the remaining were non-Beijing sublineages (n=34). The MIRU-VNTR analysis showed that Beijing isolates, most of which belonged to the modern type (n=37), formed 5 clusters and 13 individual patterns. In katG, only mutation Ser315Thr was identified. In rpoB, Ser531Leu was predominant, except for His526Arg and Leu533Pro, which were found in two isolates. A cluster of 14 Beijing strains contained these common mutations and shared the MIRU-VNTR genotype with isolates in the Thamaka district that had spread previously. Two U SIT523 isolates contained the mutations A1400G in rrs and Asp94Gly in gyrA genes, indicating a spread of XDR-TB. Most mutations were associated with drug resistance and the specific MDR Beijing and XDR-TB in U SIT523 isolates remain. This genotyping is a key tool for tracking TB transmission in the Thamaka district of Thailand.

CitSATdb: Genome-Wide Simple Sequence Repeat (SSR) Marker Database of Citrus Species for Germplasm Characterization and Crop Improvement.

Microsatellites or simple sequence repeats (SSRs) are popular co-dominant markers that play an important role in crop improvement. To enhance genomic resources in general horticulture, we identified SSRs in the genomes of eight citrus species and characterized their frequency and distribution in different genomic regions. Citrus is the world’s most widely cultivated fruit crop. We have implemented a microsatellite database, citSATdb, having the highest number (~1,296,500) of putative SSR markers from the genus Citrus, represented by eight species. The database is based on a three-tier approach using MySQL, PHP, and Apache. The markers can be searched using multiple search parameters including chromosome/scaffold number(s), motif types, repeat nucleotides (1–6), SSR length, patterns of repeat motifs and chromosome/scaffold location. The cross-species transferability of selected markers can be checked using e-PCR. Further, the markers can be visualized using the Jbrowse feature. These markers can be used for distinctness, uniformity, and stability (DUS) tests of variety identification, marker-assisted selection (MAS), gene discovery, QTL mapping, and germplasm characterization. citSATdb represents a comprehensive source of markers for developing/implementing new approaches for molecular breeding, required to enhance Citrus productivity. The potential polymorphic SSR markers identified by cross-species transferability could be used for genetic diversity and population distinction in other species.

Targeting Telomeres: Molecular Dynamics and Free Energy Simulation of Gold-Carbene Binding to DNA

DNA sequences in regulatory regions and in telomers at the ends of chromosomes frequently contain tandem repeats of guanine nucleotides that can form stacked structures stabilized by Hoogsten pairing and centrally bound monovalent cations. The replication and elongation of telomeres requires the disruption of these G-quadruplex structures. Hence, drug molecules such as gold (Au)-carbene that stabilize G-quadruplexes may also interfere with the elongation of telomeres and, in turn, could be used to control cell replication and growth. To better understand the molecular mechanism of Au-carbene binding to G-quadruplexes, we employed molecular dynamics simulations and free energy simulations. Whereas very restricted mobility of two Au-carbene ligands was found upon binding as a doublet to one side of the G-quadruplex, much larger translational and orientational mobility was observed for a single Au-carbene binding at the second G-quadruplex surface. Comparative simulations on duplex DNA in the presence of Au-carbene ligands indicates a preference for the minor groove and weaker unspecific and more salt-dependent binding than to the G-quadruplex surface. Analysis of energetic contributions reveals a dominance of nonpolar and van der Waals interactions to drive binding. The simulations can also be helpful for proposing possible modifications that could improve Au-carbene affinity and specificity for G-quadruplex binding.

Antisense Transcription across Nucleotide Repeat Expansions in Neurodegenerative and Neuromuscular Diseases: Progress and Mysteries.

Unstable repeat expansions and insertions cause more than 30 neurodegenerative and neuromuscular diseases. Remarkably, bidirectional transcription of repeat expansions has been identified in at least 14 of these diseases. More remarkably, a growing number of studies has been showing that both sense and antisense repeat RNAs are able to dysregulate important cellular pathways, contributing together to the observed clinical phenotype. Notably, antisense repeat RNAs from spinocerebellar ataxia type 7, myotonic dystrophy type 1, Huntington’s disease and frontotemporal dementia/amyotrophic lateral sclerosis associated genes have been implicated in transcriptional regulation of sense gene expression, acting either at a transcriptional or posttranscriptional level. The recent evidence that antisense repeat RNAs could modulate gene expression broadens our understanding of the pathogenic pathways and adds more complexity to the development of therapeutic strategies for these disorders. In this review, we cover the amazing progress made in the understanding of the pathogenic mechanisms associated with repeat expansion neurodegenerative and neuromuscular diseases with a focus on the impact of antisense repeat transcription in the development of efficient therapies.

Genomic characteristics and profile of microsatellite primers for Acanthogobius ommaturus by genome survey sequencing.

Acanthogobius ommaturus is one of the suitable species to study the genetic mechanism of adaptive evolution, but there are few reports on its genetics. In the present study, the genomic survey sequencing method was used to analyze the genome characters of A. ommaturus. A total of 50.50 G high-quality sequence data were obtained in the present study. From the 19-mer distribution frequency, the estimated genome size was 928.01 Mb. The calculated sequence repeat rate was about 38.31%, the heterozygosity was approximately 0.17%, and the GC% content was approximately 40.88%. Moreover, 475,724 simple sequence repeats (SSRs) were identified. Among them, dinucleotide repeats were the most (53.70% of the total SSRs), followed by tri- (35.36%), hexa- (4.59%), tetra- (4.57%) and penta- (1.77%) nucleotide repeats type. This is the first genome-wide feature of this species to be reported.

Genome survey sequencing of Atractylodes lancea and identification of its SSR markers.

Atractylodes lancea (Thunb.) DC. is a traditional Chinese medicine rich in sesquiterpenes that has been widely used in China and Japan for the treatment of viral infections. Despite its important pharmacological value, genomic information regarding A. lancea is currently unavailable. In the present study, the whole genome sequence of A. lancea was obtained using an Illumina sequencing platform. The results revealed an estimated genome size for A. lancea of 4,159.24 Mb, with 2.28% heterozygosity, and a repeat rate of 89.2%, all of which indicate a highly heterozygous genome. Based on the genomic data of A. lancea, 27,582 simple sequence repeat (SSR) markers were identified. The differences in representation among nucleotide repeat types were large, e.g., the mononucleotide repeat type was the most abundant (54.74%) while the pentanucleotide repeats were the least abundant (0.10%), and sequence motifs GA/TC (31.17%) and TTC/GAA (7.23%) were the most abundant among the dinucleotide and trinucleotide repeat motifs, respectively. A total of 93,434 genes matched known genes in common databases including 48,493 genes in the Gene Ontology (GO) database and 34,929 genes in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. This is the first report to sequence and characterize the whole genome of A. lancea and will provide a theoretical basis and reference for further genome-wide deep sequencing and SSR molecular marker development of A. lancea.

Mathematical Connection between Short Telomere Induced Senescence Calculation and Mortality Rate Data.

The last 20 years have seen a surge in scientific activity and promising results in the study of aging and longevity. Many researchers have focused on telomeres, which are composed of a series of TTAGGG repeat nucleotide sequences at the ends of each chromosome. Measurements of the length of these telomere strands show that they decrease in length with increasing age, leading many authors to propose that when the length of these telomere strands decreases sufficiently, the cells enter into a state of replicative senescence, eventually leading to disease and death. These ideas are supported by evidence that short telomere length is correlated with increased mortality. In this paper, we extend this idea to make an actual calculation of the predicted mortality rate caused by short telomere length induced senescence (STLIS). We derive a simple equation for the mathematical relationship between telomere length and mortality rate. Using only three parameters based on telomere length measurement data of Canadians, we have calculated both the magnitude and the age dependence of the mortality rate for both men and women. We show that these calculated data are in good quantitative agreement with the actual number of Canadians that die. This agreement demonstrates the quantitative correlation between the mortality calculated by the STLIS model and the mortality of the major diseases of aging (e.g., cardiovascular disease, many cancers and diabetes mellitus), which dominate human mortality. This result represents significant progress in our understanding of the factors behind the cause of aging.

ISSR-PCR approach as a means of studying genetic variation in moth bean (Vigna aconitifolia (Jacq.) Maréchal)

Abstract Moth bean (Vigna aconitifolia (Jacq.) Marechal) is considered as one of the underutilised food legume vis-a-vis the Indian sub-continent. The seeds of this legume are a rich source of soluble proteins for nutrition. PCR-based marker characterizations for assessing genetic diversity are almost negligible in this crop. Present investigation, therefore, employed inter simple sequence repeat (ISSR) markers to assess the genetic diversity amongst 25 accessions. The di and penta nucleotide repeat markers showed higher polymorphism on amplification than those observed using trinucleotide markers. ISSR marker UBC-811 detected 36.15% polymorphism while UBC-879 detected 61.03% with an average of 9.7% on each primer. Based on dendrogram analysis, a logical similarity between molecular data and geographical origin was corroborated. The accessions cultivated in western coast of India were closely related than north-western and arid regions of India. The accessions like Maru, Jadiya and Jawala can be recommended as parents for further breeding programmes to procure additional variation complexes in moth bean genetic stocks.

The Spillover of African Swine Fever in Western Poland Revealed Its Estimated Origin on the Basis of O174L, K145R, MGF 505-5R and IGR I73R/I329L Genomic Sequences.

The African swine fever epidemic occurred in Poland at the beginning of 2014 and, up to date, the disease has been spreading mainly in the eastern part of the country. Unexpectedly, in November 2019 an infected wild boar case was confirmed in Lubuskie voivodship in western Poland. During the following weeks, several dozen African swine fever virus (ASFV)-positive animals were notified in the neighboring area, causing severe concern regarding further spread of the disease to the mostly pig-dense region in Poland, namely, Wielkopolskie voivodship. Moreover, almost a year after, several infected wild boar cases were confirmed for the first time in Germany, just beyond the Polish border, sending out a shock wave through the global pig market. The whole genome sequence of ASFV, isolated from the first case of ASF in western Poland, and three selected viruses from other affected areas, revealed the tandem repeat and single nucleotide polymorphism (SNP) variations in reference to the Georgia 2007/1 strain. These data, supported by the conventional sequencing of selected genomic regions from a total of 154 virus samples isolated between 2017 and 2020 in Poland, shed a new light on pathogen epidemiology. The sequence variations within the O174L gene detected in this study showed that cases identified in western Poland might be originating from the so-called southern Warsaw cluster. Moreover, the viruses originating from the northern Warsaw cluster do not possess single nucleotide polymorphism (SNP) mutations within the K145R and MGF 505-5R genes, which are specific to all of the other Polish ASFV strains. These results led to a conclusion of their distinct origin. Supporting these results, the nucleotide sequencing of I73R/I329L intergenic region revealed its new, previously undescribed variant, called IGR IV, with an additional three tandem repeats of 10 nucleotides in comparison to the reference sequence of the Georgia 2007/1 strain.

Genotyping of Russian isolates of fungal pathogen Trichophyton rubrum, based on simple sequence repeat and single nucleotide polymorphism.

The Trichophyton rubrum species group consists of prevalent causative agents of human skin, nail and hair infections, including Trubrum sensu stricto and Tviolaceum, as well as other less well-established or debatable taxa like Tsoudanense, Tkuryangei and Tmegninii. Our previous study provided limited evidence in favour of the existence of two genetic lineages in the Russian Trubrum sensu stricto population. We aimed to study the genetic structure of the Russian population of Trubrum and to identify factors shaping this structure. We analysed the polymorphism of 12 simple sequence repeat (SSR or microsatellite) markers and single nucleotide polymorphism in the TERG_02941 protein-coding gene in 70 Trubrum isolates and performed a phylogenomic reconstruction. All three types of data provided conclusive evidence that the population consists of two genetic lineages. Clustering, performed by means of microsatellite length polymorphism analysis, was strongly dependent on the number of nucleotide repeats in the 5'-area of the fructose-1,6-bisphosphate aldolase gene. Analysis of molecular variance (AMOVA) on the basis of SSR typing data indicated that 22%-48% of the variability was among groups within Trubrum. There was no clear connection of population structure with types of infection, places of geographic origin, aldolase gene expression or urease activity. Our results suggest that the Russian population of Trubrum consists of two cosmopolitan genetic lineages.

Kinship analysis on single cells after whole genome amplification

Short Tandem Repeat (STR-) and Single Nucleotide Polymorphism (SNP-) genotyping have been extensively studied within forensic kinship analysis. Nevertheless, no results have been reported on kinship analysis after whole genome amplification (WGA) of single cells. This WGA step is a necessary procedure in several applications, such as cell-based non-invasive prenatal testing (cbNIPT) and pre-implantation genetic diagnosis (PGD). In cbNIPT, all putative fetal cells must be discriminated from maternal cells after enrichment from whole blood. This study investigates the efficacy and evidential value of STR- and SNP-genotyping methods for the discrimination of 24 single cells after WGA, within three families. Formaldehyde-fixed and unfixed cells are assessed in offspring-parent duos and offspring-mother-father trios. Results demonstrate that both genotyping methods can be used in all tested conditions and scenarios with 100% sensitivity and 100% specificity, with a similar evidential value for fixed and unfixed cells. Moreover, sequence-based SNP-genotyping results in a higher evidential value than length-based STR-genotyping after WGA, which is not observed using high-quality offspring bulk DNA samples. Finally, it is also demonstrated that the availability of the DNA genotypes of both parents strongly increases the evidential value of the results.

Cross genera transferability of microsatellite markers from other members of family Bignoniaceae to Tecomella undulata (Sm.) Seem

Transferability of SSR markers from related species/genera is a rapid and cost-effective method to enhance genomic database of genetically neglected crops. Monotypic genus Tecomella undulata, a pharmacologically important endangered timber tree, found in hyper arid regions is one such genetic resource limited species, wherein neither SSR identification nor cross transferability studies have been initiated. Transferability of 69 cross-genera SSR primers selected from other members of family bignoniaceae (Incarvillea sinensis, I. mairei, Jacaranda copaia, Tabebuia aurea and Arrabidaea chica) showed 40.58% (28 primers) transferability amplifying 1–10 amplicons ranging in the size from 80to 600 bp in 24 accessions of T. undulata. Per cent marker transferability varied from 66.67% (A. chica) to 38.89% (J. copaia) among the genera tested while no transferable markers could be identified from I. sinensis. Transferability of simple di-/tri-nucleotide repeat (34.2–37.5%) and complex nucleotide sequence (60.0%) based SSR motifs were higher as compared to penta-/hexa-nucleotide based repeat motifs. Within the 28 transferable markers, 26 (92.86%) were polymorphic while 2 (7.14%) were monomorphic. Average number of alleles per primer was 5.39 while average polymorphic alleles per locus were 5.11. The PIC values of SSR markers ranged from 0.23 to 0.85 with an average of 0.59. Efficiency of the identified cross species SSR markers was tested for diversity analysis using a set of 24 trees (8 each of yellow, orange and red flower colour morphotypes). The 24 samples clustered into six groups with similarity coefficient ranging from 0.29 to 0.67, while six accessions formed distinct out groups. Distribution of accessions in the SSR dendrogram showed no correlation with flower colour. More diversity (99%) was observed within populations compared to that existing between populations (1%). The present study can enrich the genomic background of T. undulata by identifying suitable polymorphic markers from confamiliar species that can be applied for genetic diversity studies, mapping of QTLs and cultivar identification as well.

The RNA exosome complex degrades expanded hexanucleotide repeat RNA in C9orf72 FTLD/ALS

Nucleotide repeat expansions in the C9orf72 gene cause frontotemporal lobar degeneration (FTLD) and amyotrophic lateral sclerosis (ALS). Transcribed repeat RNA accumulates within RNA foci and is also translated into toxic dipeptide repeat proteins (DPR). The mechanism of repeat RNA accumulation, however, remains unclear. The RNA exosome complex is a multimeric ribonuclease involved in degradation of defective RNA. Here, we uncover the RNA exosome as a major degradation complex for pathogenic C9orf72‐derived repeat RNA. Knockdown of EXOSC10, the catalytic subunit of the complex, enhanced repeat RNA and DPR protein expression levels. RNA degradation assays confirmed that EXOSC10 can degrade both sense and antisense repeats. Furthermore, EXOSC10 reduction increased RNA foci and repeat transcripts in patient‐derived cells. Cells expressing toxic poly‐GR or poly‐PR proteins accumulate a subset of small nucleolar RNA precursors, which are physiological substrates of EXOSC10, as well as excessive repeat RNA, indicating that arginine‐rich DPR proteins impair the intrinsic activity of EXOSC10. Collectively, arginine‐rich DPR‐mediated impairment of EXOSC10 and the RNA exosome complex compromises repeat RNA metabolism and may thus exacerbate C9orf72‐FTLD/ALS pathologies in a vicious cycle.

Uniparental disomy of chromosome 21: A statistical approach and application in paternity tests

Considering the overall frequency of paternity investigation cases including mutational events, there is a real possibility that at least a fraction of all inconsistencies reported in paternity cases are caused not by polymerase slippage mutations, but to chromosomic abnormalities, as Uniparental Disomy (UPD). We report here the investigation of a trio paternity case (mother, child and alleged father), with observed inconsistencies that can alternatively be explained by occurrence of maternal uniparental isodisomy of chromosome 21 (miUPD21). A total of 350 short tandem repeat (STR) and single nucleotide polymorphism (SNP) markers were tested, statistically suggesting true biological linkage within the trio. Additionally, we propose miUPD21 explains, with significantly greater probability, the occurrence of detected inconsistencies, when compared to alternative hypothesis of multiple and simultaneous slippage mutations. Similar cases could have their statistical conclusions improved or even altered by including unusual chromosomal segregation patterns in the hypothesis formulation, as well as in mathematical calculations. Such reports of allelic inconsistencies being explained by chromosomal alterations are common in clinical genetics, and such situations might have impact on forensic investigation.

Relatively semi-conservative replication and a folded slippage model for short tandem repeats

BackgroundThe ubiquitous presence of short tandem repeats (STRs) in virtually all genomes implicates their functional relevance, while a widely-accepted definition of STR is yet to be established. Previous studies majorly focus on relatively longer STRs, while shorter repeats were generally excluded. Herein, we have adopted a more generous criteria to define shorter repeats, which has led to the definition of a much larger number of STRs that lack prior analysis. Using this definition, we analyzed the short repeats in 55 randomly selected segments in 55 randomly selected genomic sequences from a fairly wide range of species covering animals, plants, fungi, protozoa, bacteria, archaea and viruses.ResultsOur analysis reveals a high percentage of short repeats in all 55 randomly selected segments, indicating that the universal presence of high-content short repeats could be a common characteristic of genomes across all biological kingdoms. Therefore, it is reasonable to assume a mechanism for continuous production of repeats that can make the replicating process relatively semi-conservative. We have proposed a folded replication slippage model that considers the geometric space of nucleotides and hydrogen bond stability to explain the mechanism more explicitly, with improving the existing straight-line slippage model. The folded slippage model can explain the expansion and contraction of mono- to hexa- nucleotide repeats with proper folding angles. Analysis of external forces in the folding template strands also suggests that expansion exists more commonly than contraction in the short tandem repeats.ConclusionThe folded replication slippage model provides a reasonable explanation for the continuous occurrences of simple sequence repeats in genomes. This model also contributes to the explanation of STR-to-genome evolution and is an alternative model that complements semi-conservative replication.

A Forensic Genomics Approach for the Identification of Sister Marija Crucifiksa Kozulić.

Sister Marija Krucifiksa Kozulić (1852–1922) was a Croatian nun who is in consideration for beatification by the Vatican, which is facilitated by the identification of her 20th-century remains. Sister Marija was buried in a tomb in Rijeka, Croatia, along with other nuns including her biological sister, Tereza Kozulić (1861–1933). When the remains were exhumed in 2011, they were found in a deteriorated state and commingled with several other sets of remains. Thus, mitochondrial genome sequencing of the long bones was performed to sort the remains by mitochondrial haplotype. Two similar but unique haplotypes belonging to haplogroup H1bu were identified, and samples from these bones were subjected to autosomal short tandem repeat (STR) and single nucleotide polymorphism (SNP) sequencing. Although only partial profiles were obtained, the data were sufficient for kinship analysis with the profile of a paternal niece of Sister Marija (Fides Kozulić). The data indicate that it is 574,195-fold more likely that the two sets of skeletal remains represent 2nd-degree relatives of Fides than sisters who are unrelated to Fides. Although it is impossible to discern which set of remains belongs to Marija and which belongs to Tereza, forensic genomics methods have enabled identification of the sisters.