Updated Brazilian STR allele frequency data using over 100,000 individuals: An analysis of CSF1PO, D3S1358, D5S818, D7S820, D8S1179, D13S317, D16S539, D18S51, D21S11, FGA, Penta D, Penta E, TH01, TPOX and vWA loci

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.

Forensic characterization of Brazilian regional populations through massive parallel sequencing of 124 SNPs included in HID ion Ampliseq Identity Panel

Use of Massive Parallel Sequencing (MPS) techniques has been investigated by forensic community aiming introduction of such methods in routine forensic casework analyses. Interesting features presented by MPS include high-throughput, ability to simultaneous genotyping of significant number of samples and forensic markers, workflow automation, among others. Emergence of single nucleotide polymorphism (SNP) as forensic relevant markers was facilitated in this process, since concurrent typing of larger marker sets is necessary for obtaining same levels of individual discrimination provided by other marker categories. In this context, HID Ion Ampliseq Identity Panel is a commercial solution with forensic purposes comprising simultaneous analysis of 90 highly informative autosomal SNPs and 34 Y -chromosome superior clade SNPs for male lineage haplotyping. SNP typing can be obtained with smaller amplicons, and this panel was designed for efficient processing of critical or challenging forensic samples. In this work, a sample of 432 individuals from all five Brazilian geopolitical regions was evaluated with this panel, in order to access feasibility of this panel use in a national basis. Results obtained for all five regions, including forensic parameters, show that this marker set can be efficiently employed for Brazilian nationals in human identification or kinship determination applications, due to high levels of genetic discriminative information content displayed by Brazilians. Interpopulation comparison studies were executed among Brazilian regional populations and 26 worldwide populations, in order to access genetic stratification occurrence. Some levels of population structure were identified, and impact on database design was discussed. Y-chromosome haplotyping of Brazilian samples revealed high levels of European ancestry in Brazilian male lineages, and utility of haplotyping in real forensic casework is addressed. Finally, genotyping and sequencing efficiency with this panel were addressed, as an effort to appraise the adequacy of this panel use in Brazilian national forensic demands.

Evaluating DNA evidence in a genetically complex population

In forensic genetics, the likelihood ratio (LR), measuring the value of DNA profile evidence, is computed from a database of allele frequencies. Here, we address the choice of database and adjustments for population structure and sample size in the context of Brazil. The Brazilian population underwent a complex process of colonization, migration and mating, which created an admixed genetic composition that makes it difficult to obtain an appropriate database for a given case. National databases are now available, as well as databases for many Brazilian states. However, those databases are not statistically random samples, and state boundaries may not accurately reflect the sub-structuring of genetic diversity. We compared the LR calculated using the relevant state-specific database with the statistics calculated when a national database and when international databases were used. We evaluated two methods of adjustment for population structure, due to Wright [13] and Balding and Nichols [14]. We also considered two adjustments for database sample size: the Balding size bias correction [15] and a minimum allele frequency [16]. Our results show that the use of a national database with the Balding and Nichols adjustment and θ = 0.002 generated lower LR values than did the state-specific database in more than 50% of the profiles simulated using the state-based allele frequencies, while θ = 0.01 produced lower LRs for more than 90% of the profiles. We conclude that the utilization of a national database for Brazilian cases can be justified in association with the appropriate adjustment for population structure.

Genetic analysis of Southern Brazil subjects using the PowerSeq™ AUTO/Y system for short tandem repeat sequencing

With the advent of Next-Generation Sequencing technology, sequencing of short tandem repeats (STRs) allows for a more detailed analysis when compared to size-based fragment methods (capillary electrophoresis-CE). The implementation of high-throughput sequencing can help uncover deeper genetic diversities of different populations. Subjects from the South region of Brazil present a particular and more homogeneous ancestry background when compared to other regions of the country. Both autosomal and Y- STRs have been analyzed in these individuals; however, all analyses published to date encompass data from CE-based fragment analysis. In this study, a genetic analysis of 59 individuals from Southern Brazil was performed on STR sequences. Forensically relevant STRs were PCR-enriched using a prototype of the PowerSeq™ AUTO/Y system (Promega Corp.). Next-generation sequencing was performed on an Illumina MiSeq instrument. The raw data (FASTQ files) were processed using a custom designed sequence processing tool, Altius. Isoalleles, which are sequence-based allelic variants that do not differ in length, were observed in nine autosomal and in six Y- STRs from the core global forensic marker set. The number of distinctive alleles based on sequence was higher when compared to those based on length, 37.3% higher in autosomal STRs and 13.8% higher in Y-STRs. The most polymorphic autosomal locus was D12S391, which presented 38 different sequence-based alleles. Among the loci in the Y chromosome, DYS389II presented the highest number of isoalleles. In comparison to CE analysis, Observed and Expected Heterozygosity, Polymorphic Information Content (PIC) and Genetic Diversity also presented higher values when the alleles were analyzed based on their sequence. For autosomal loci, Polymorphic Information Content (PIC) was 2.6% higher for sequence-based data. Diversity was 9.3% and 6.5% higher for autosomal and Y markers, respectively. In the analysis of the repeat structures for the STR loci, a new allele variant was found for allele 18 in the vWA locus. The STR flanking regions were also further investigated and sixteen variations were observed at nine autosomal STR loci and one Y-STR locus. The results obtained in this study demonstrate the importance of genetic analysis based on sequencing and highlight the diversity of the South Brazilian population when characterized by STR sequencing.

A new panel of SNPs to assess thyroid carcinoma risk: a pilot study in a Brazilian admixture population

BACKGROUND

Thyroid cancer is a common malignant disease of the endocrine system with increasing incidence rates over the last few decades. In this study, we sought to analyze the possible association of 45 single nucleotide polymorphisms (SNPs) with thyroid cancer in a population from Rio Grande do Norte, Brazil.

METHODS

Based on histological analysis by a pathologist, 80 normal thyroid specimens of tissue adjacent to thyroid tumors were obtained from the biobank at the Laboratory of Pathology of Liga Norte Riograndense Contra o Câncer, Natal, RN. Patient samples were then genotyped using the MassARRAY platform (Sequenon, Inc) followed by statistical analysis employing the SNPassoc package in R program. The genotypic frequencies of all 45 SNPs obtained from the International HapMap Project database and based on data from the ancestral populations of European and African origin were used to compose the control study group.

RESULTS

In our study, the following 9 SNPs showed significant differences in their frequency when comparing the study and control groups: rs3744962, rs258107, rs1461855, rs4075022, rs9943744, rs4075570, rs2356508, rs17485896, and rs2651339. Furthermore, the SNPs rs374492 C/T and rs258107 C/T were associated with a relative risk for thyroid carcinoma of 3.78 (p = 6.27 × 10e^-5) and 2.91 (p = 8.27 × 10e^-5), respectively, after Bonferroni's correction for multiple comparisons.

CONCLUSIONS

These nine polymorphisms could be potential biomarkers of predisposition to thyroid carcinoma in the population from Rio Grande do Norte. However, complementary studies including a control group with samples obtained from healthy subjects in Rio Grande do Norte state, should be conducted to confirm these results.

Introduction of the Python script STRinNGS for analysis of STR regions in FASTQ or BAM files and expansion of the Danish STR sequence database to 11 STRs

This work introduces the in-house developed Python application STRinNGS for analysis of STR sequence elements in BAM or FASTQ files. STRinNGS identifies sequence reads with STR loci by their flanking sequences, it analyses the STR sequence and the flanking regions, and generates a report with the assigned SNP-STR alleles. The main output file from STRinNGS contains all sequences with read counts above 1% of the total number of reads per locus. STR sequences are automatically named according to the nomenclature used previously and according to the repeat unit definitions in STRBase (http://www.cstl.nist.gov/strbase/). The sequences are named with (1) the locus name, (2) the length of the repeat region divided by the length of the repeat unit, (3) the sequence(s) of the repeat unit(s) followed by the number of repeats and (4) variations in the flanking regions. Lower case letters in the main output file are used to flag sequences with previously unknown variations in the STRs. SNPs in the flanking regions are named by their "rs" numbers and the nucleotides in the SNP position. Data from 207 Danes sequenced with the Ion Torrent™ HID STR 10-plex that amplified nine STRs (CSF1PO, D3S1358, D5S818, D7S820, D8S1179, D16S539, TH01, TPOX, vWA), and Amelogenin was analysed with STRinNGS. Sequencing uncovered five common SNPs near four STRs and revealed 20 new alleles in the 207 Danes. Three short homopolymers in the D8S1179 flanking regions caused frequent sequencing errors. In 29 of 3726 allele calls (0.8%), sequences with homopolymer errors were falsely assigned as true alleles. An in-house developed script in R compensated for these errors by compiling sequence reads that had identical STR sequences and identical nucleotides in the five common SNPs. In the output file from the R script, all SNP-STR haplotype calls were correct. The 207 samples and six additional samples were sequenced for D3S1358, D12S391, and D21S11 using the 454 GS Junior platform in this and a previous work. Overall, next generation sequencing (NGS) of the 11 STRs lowered the mean match probability 386 times and increased the typical paternity indexes (i.e. the geometric mean) for trios and duos 47 and 23 times, respectively, compared to the traditional PCR-CE typing of the same population.

Exploring the applicability of analysing X chromosome STRs in Brazilian admixed population

Kinship and parentage analyses always involve one sample being compared to another sample or a few samples with a specific relationship question in mind. In most cases, the analysis of autosomal STR markers is sufficient to determine the genetic kinship. However, when genetic profiles are reconstructed from supposed relatives, for whom the family configuration available for analysis is deficient, the examination may be inconclusive. This study reports practical examples of actual cases analysing the efficiency of the chromosome X STR (STR-ChrX) markers. Three cases with different degrees of efficiency and impact were selected as follows: the identification of two charred bodies in a traffic accident, in which the family setting available was not complete, and one filiation analysis resulting from rape. This is the first paper reporting the use of the multiplex STR 12 ChrX in actual cases using the software Familias 1.8 and Brazilian regional frequency data. Our study clarifies the complex analysis using this powerful tool for professionals in the forensic science community, for both civil and criminal justice. We also discuss state-of-the-art ChrX STR markers and its implications and applications for legal procedures. The data presented here should be used in other studies of complex cases to improve the progress of the current justice system.

Mutation rate estimation for 15 autosomal STR loci in a large population from Mainland China

STR, short tandem repeats, are well known as a type of powerful genetic marker and widely used in studying human population genetics. Compared with the conventional genetic markers, the mutation rate of STR is higher. Additionally, the mutations of STR loci do not lead to genetic inconsistencies between the genotypes of parents and children; therefore, the analysis of STR mutation is more suited to assess the population mutation. In this study, we focused on 15 autosomal STR loci. DNA samples from a total of 42,416 unrelated healthy individuals (19,037 trios) from the population of Mainland China collected between Jan 2012 and May 2014 were successfully investigated. In our study, the allele frequencies, paternal mutation rates, maternal mutation rates and average mutation rates were detected. Furthermore, we also investigated the relationship between paternal ages, maternal ages, area, the time of pregnancy and average mutation rate. We found that the paternal mutation rate was higher than the maternal mutation rate and the paternal, maternal, and average mutation rates had a positive correlation with paternal age, maternal age and the time of pregnancy respectively. Additionally, the average mutation rate of coastal areas was higher than that of inland areas.

•

15 autosomal STR loci were analyzed in a large scale population in Mainland China.

•

With increases of repeat units the event of loss is more than that of gain.

•

The paternal mutation rate is higher than the maternal mutation rate.

•

There is a positive correlation between mutation rates and age and pregnancy time.

•

The average mutation rates of coastal areas are higher than that of inland areas.

Population data of 16 autosomal STR loci of the Powerplex ESX 17 System in a Brazilian Population from the State of São Paulo

The State of São Paulo is the most populous state in Brazil, including approximately one fifth of the population of the country. In addition to a strong economy, the state has relatively good social indicators when compared with the rest of the country. The capital city, also called São Paulo, is the sixth largest city in the world. Its population is considered the most multicultural and racially mixed in Brazil. Currently, the largest populations in São Paulo are of Italian, Lebanese, Spanish and Japanese origin, and the state has the largest number of Northeasterners outside of the Northeast region. This population structure may lead to a particular genotype frequency. In this context, the formation of a new database containing the allele frequencies of five new genetic markers (D2S441, D10S1248, D22S1045, D1S1656 and D12S391) in a sample population is relevant. The allele frequencies of 16 STR loci, including the five new European Standard Set (ESS) loci, were calculated in a sample of 1088-1098 unrelated individuals, who geographically represent the Capital city.