TEX13B is essential for metabolic reprogramming during germ cell differentiation

STUDY QUESTION

What is the functional significance of Tex13b in male germ cell development and differentiation?

SUMMARY ANSWER

Tex13b regulates male germ cell differentiation by metabolic reprogramming during spermatogenesis.

WHAT IS KNOWN ALREADY

Studies in mice and humans suggest that TEX13B is a transcription factor and is exclusively expressed in germ cells.

STUDY DESIGN, SIZE, DURATION

We sequenced the coding regions of TEX13B in 628 infertile men and 427 ethnically matched fertile control men. Further, to identify the molecular function of Tex13b, we created a Tex13b knockout and conditional overexpression system in GC-1spg (hereafter, GC-1) cells.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Our recent exome sequencing study identified novel candidate genes for male infertility. TEX13B was found to be one of the potential candidates, hence we explored the role of TEX13B in male infertility within a large infertile case-control cohort. We performed functional analyses of Tex13b in a GC-1 cell line using CRISPR-Cas9. We differentially labelled the cell proteins by stable isotope labelling of amino acids in cell culture (SILAC) and performed mass spectrometry-based whole-cell proteomics to identify the differential protein regulation in knockout cells compared to wild-type cells. We found that Tex13b knockout leads to downregulation of the OXPHOS complexes and upregulation of glycolysis genes, which was further validated by western blotting. These results were further confirmed by respirometry analysis in Tex13b knockout cells. Further, we also performed a conditional overexpression of TEX13B in GC-1 cells and studied the expression of OXPHOS complex proteins by western blotting.

MAIN RESULTS AND THE ROLE OF CHANCE

We identified a rare variant, rs775429506 (p.Gly237Glu), exclusively in two non-obstructive-azoospermia (NOA) men, that may genetically predispose these men for infertility. Further, we demonstrated that Tex13b functions in the transcription regulation of OXPHOS complexes.

LARGE SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

We examined the function of Tex13b in GC-1 in vitro by knocking out and conditional overexpression, for understanding the function of Tex13b in germ cells. Unfortunately, this could not be replicated in either an animal model or in patient-derived tissue due to the non-availability of an animal model or patient's testis biopsies.

WIDER IMPLICATIONS OF THE FINDINGS

This study identified that Tex13b plays an important role in male germ cell development and differentiation. The findings of this study would be useful for screening infertile males with spermatogenic failure and counselling them before the implementation of assisted reproduction technique(s).

STUDY FUNDING/COMPETING INTEREST(S)

Funding was provided by the Council of Scientific and Industrial Research (CSIR) under the network project (BSC0101 and MLP0113) and SERB, the Department of Science and Technology, Government of India (J C Bose Fellowship: JCB/2019/000027). The authors do not have any competing interest.

The genetic identity of the Vedda: A language isolate of South Asia

Linguistic data from South Asia identified several language isolates in the subcontinent. The Vedda, an indigenous population of Sri Lanka, are the least studied amongst them. Therefore, to understand the initial peopling of Sri Lanka and the genetic affinity of the Vedda with other populations in Eurasia, we extensively studied the high-resolution autosomal and mitogenomes from the Vedda population of Sri Lanka. Our autosomal analyses suggest a close genetic link of Vedda with the tribal populations of India despite no evidence of close linguistic affinity, thus suggesting a deep genetic link of the Vedda with these populations. The mitogenomic analysis supports this association by pointing to an ancient link with Indian populations. We suggest that the Vedda population is a genetically drifted group with limited gene flow from neighbouring Sinhalese and Sri Lankan Tamil populations. Interestingly, the genetic ancestry sharing of Vedda neglects the isolation-by-distance model. Collectively, the demography of Sri Lanka is unique, where Sinhalese and Sri Lankan Tamil populations excessively admixed, whilst Vedda largely preserved their isolation and deep genetic association with India.

Mitochondria in biology and medicine - 2023

Mitochondria are an indispensable part of the cell that plays a crucial role in regulating various signaling pathways, energy metabolism, cell differentiation, proliferation, and cell death. Since mitochondria have their own genetic material, they differ from their nuclear counterparts, and dysregulation is responsible for a broad spectrum of diseases. Mitochondrial dysfunction is associated with several disorders, including neuro-muscular disorders, cancer, and premature aging, among others. The intricacy of the field is due to the cross-talk between nuclear and mitochondrial genes, which has also improved our knowledge of mitochondrial functions and their pathogenesis. Therefore, interdisciplinary research and communication are crucial for mitochondrial biology and medicine due to the challenges they pose for diagnosis and treatment. The ninth annual conference of the Society for Mitochondria Research and Medicine (SMRM)- India, titled "Mitochondria in Biology and Medicine" was organized at the Centre for DNA Fingerprinting and Diagnostics (CDFD), Hyderabad, India, on June 21-23, 2023. The latest advancements in the field of mitochondrial biology and medicine were discussed at the conference. In this article, we summarize the entire event for the benefit of researchers working in the field of mitochondrial biology and medicine.

The maternal genetic origin and diversity of the extant populations of the Ladakh region in India

Ladakh lies at a strategic location between the Indus River valley and the Hindu Khush Mountains, which makes the "Land of high passes" one of the major routes of movement. Through the years the region has faced multi-layered cultural movements, genetic assimilation and demographic changes. The initial settlement in the years goes back to the early Neolithic age and still continues despite its harsh, unhospitable and cold climate. Previous studies mostly covered the patrilineal markers of the region and an in-depth study lacked to represent the matrilineal ancestry and possible genetic inflow in the region. Hence, our current study first time generated complete mitogenomes of 108 unrelated individuals from Ladakh belonging to three population groups namely, Changpa (n = 38), Brokpa (n = 32) and Monpa (n = 38). In the in-depth analysis, we found that the mitogenome of the three Ladakhi groups are highly diverse in terms of maternal haplogroup distribution carrying lineages specific to East Asia (M9a), Tibbet (A21) and South Asia (M3, M30, U2). In our analysis we found that Changpa and Monpa probably have shared maternal ancestry compared to Brokpa, which is very distinct and also later suffered possible historical Bottleneck. Bayesian evolutionary and Network analysis indicates more ancient maternal lineage of Changpa and Monpa in terms of M9a haplotypes, but they also share some genetic history with Tibeto-Burman speakers in past. These findings conclusively indicate possible matrilineal genetic inflow in Ladakh from three directions, primarily from East Asia or South East Asia during post-glacial, West Eurasia and also from South Asia.

Genetic Affinities and Adaptation of the South-West Coast Populations of India

Evolutionary event has not only altered the genetic structure of human populations but also associated with social and cultural transformation. South Asian populations were the result of migration and admixture of genetically and culturally diverse groups. Most of the genetic studies pointed to large-scale admixture events between Ancestral North Indian (ANI) and Ancestral South Indian (ASI) groups, also additional layers of recent admixture. In the present study, we have analyzed 213 individuals inhabited in South-west coast India with traditional warriors and feudal lord status and historically associated with migratory events from North/North West India and possible admixture with West Eurasian populations, whose genetic links are still missing. Analysis of autosomal Single Nucleotide Polymorphism (SNP) markers suggests that these groups possibly derived their ancestry from some groups of North West India having additional Middle Eastern genetic components. Higher distribution of West Eurasian mitochondrial haplogroups also points to female-mediated admixture. Estimation of Effective Migration Surface (EEMS) analysis indicates Central India and Godavari basin as a crucial transition zone for population migration from North and North West India to South-west coastal India. Selection screen using 3 distinct outlier-based approaches revealed genetic signatures related to Immunity and protection from Viral infections. Thus, our study suggests that the South-west coastal groups with traditional warriors and feudal lords' status are of a distinct lineage compared to Dravidian and Gangetic plain Indo-Europeans and are remnants of very early migrations from North West India following the Godavari basin to Karnataka and Kerala.

Neuromuscular disease genetics in under-represented populations: increasing data diversity

Neuromuscular diseases (NMDs) affect ∼15 million people globally. In high income settings DNA-based diagnosis has transformed care pathways and led to gene-specific therapies. However, most affected families are in low-to-middle income countries (LMICs) with limited access to DNA-based diagnosis. Most (86%) published genetic data is derived from European ancestry. This marked genetic data inequality hampers understanding of genetic diversity and hinders accurate genetic diagnosis in all income settings. We developed a cloud-based transcontinental partnership to build diverse, deeply-phenotyped and genetically characterized cohorts to improve genetic architecture knowledge, and potentially advance diagnosis and clinical management. We connected 18 centres in Brazil, India, South Africa, Turkey, Zambia, Netherlands and the UK. We co-developed a cloud-based data solution and trained 17 international neurology fellows in clinical genomic data interpretation. Single gene and whole exome data were analysed via a bespoke bioinformatics pipeline and reviewed alongside clinical and phenotypic data in global webinars to inform genetic outcome decisions. We recruited 6001 participants in the first 43 months. Initial genetic analyses 'solved' or 'possibly solved' ∼56% probands overall. In-depth genetic data review of the four commonest clinical categories (limb girdle muscular dystrophy, inherited peripheral neuropathies, congenital myopathy/muscular dystrophies and Duchenne/Becker muscular dystrophy) delivered a ∼59% 'solved' and ∼13% 'possibly solved' outcome. Almost 29% of disease causing variants were novel, increasing diverse pathogenic variant knowledge. Unsolved participants represent a new discovery cohort. The dataset provides a large resource from under-represented populations for genetic and translational research. In conclusion, we established a remote transcontinental partnership to assess genetic architecture of NMDs across diverse populations. It supported DNA-based diagnosis, potentially enabling genetic counselling, care pathways and eligibility for gene-specific trials. Similar virtual partnerships could be adopted by other areas of global genomic neurological practice to reduce genetic data inequality and benefit patients globally.

Novel MYBPC3 Mutations in Indian Population with Cardiomyopathies

Background

Mutations in Myosin Binding Protein C (MYBPC3) are one of the most frequent causes of cardiomyopathies in the world, but not much data are available in India.

Methods

We carried out targeted direct sequencing of MYBPC3 in 115 hypertrophic (HCM) and 127 dilated (DCM) cardiomyopathies against 197 ethnically matched healthy controls from India.

Results

We detected 34 single nucleotide variations in MYBPC3, of which 19 were novel. We found a splice site mutation [(IVS6+2T) T>G] and 16 missense mutations in Indian cardiomyopathies [5 in HCM; E258K, T262S, H287L, R408M, V483A: 4 in DCM; T146N, V321L, A392T, E393K and 7 in both HCM and DCM; L104M, V158M, S236G, R272C, T290A, G522E, A626V], but those were absent in 197 normal healthy controls. Interestingly, we found 7 out of 16 missense mutations (V158M, E258K, R272C, A392T, V483A, G522E, and A626V) in MYBPC3 were altering the evolutionarily conserved native amino acids, accounted for 8.7% and 6.3% in HCM and DCM, respectively. The bioinformatic tools predicted that those 7 missense mutations were pathogenic. Moreover, the co-segregation of those 7 mutations in families further confirmed their pathogenicity. Remarkably, we also identified compound mutations within the MYBPC3 gene of 6 cardiomyopathy patients (5%) with more severe disease phenotype; of which, 3 were HCM (2.6%) [(1. K244K + E258K + (IVS6+2T) T>G); (2. L104M + G522E + A626V); (3. P186P + G522E + A626V]; and 3 were DCM (2.4%) [(1. 5'UTR + A392T; 2. V158M+G522E; and 3.V158M + T262T + A626V].

Conclusion

The present comprehensive study on MYBPC3 has revealed both single and compound mutations in MYBPC3 and their association with disease in Indian Population with Cardiomyopathies. Our findings may perhaps help in initiating diagnostic strategies and eventually recognizing the targets for therapeutic interventions.

Deciphering the West Eurasian Genetic Footprints in Ancient South India

Since 2006, Pattanam coastal village of the Ernakulam District in Kerala, India, has witnessed multi-disciplinary archaeological investigations in collaboration with leading research institutions across the world. The results confirm that the Pattanam site could be an integral part of the lost ancient port of Muziris, which, as per the material evidence from Pattanam and its contemporary sites, played an important role in the transoceanic exchanges between 100 BCE (Before Common Era) and 300 CE (Common Era). So far, the material evidence with direct provenance to the maritime exchanges related to ancient cultures of the Mediterranean, West Asian, Red Sea, African, and Asian regions have been identified at Pattanam. However, the genetic evidence supporting the impact of multiple cultures or their admixing is still missing for this important archaeological site of South India. Hence, in the current study, we tried to infer the genetic composition of the skeletal remains excavated from the site in a broader context of South Asian and worldwide maternal affinity. We applied the MassArray-based genotyping approach of mitochondrial makers and observed that ancient samples of Pattanam represent a mixed maternal ancestry pattern of both the West Eurasian ancestry and the South Asian ancestry. We observed a high frequency of West Eurasian haplogroups (T, JT, and HV) and South Asian-specific mitochondrial haplogroups (M2a, M3a, R5, and M6). The findings are consistent with the previously published and ongoing archaeological excavations, in which material remains from over three dozen of sites across the Indian Ocean, Red Sea, and Mediterranean littoral regions have been unearthed. This study confirms that people belonging to multiple cultural and linguistic backgrounds have migrated, probably settled, and eventually died on the South-western coast of India.

Exome sequencing and functional analyses revealed CETN1 variants leads to impaired cell division and male fertility

Human spermatogenesis requires an orchestrated expression of numerous genes in various germ cell subtypes. Therefore, the genetic landscape of male infertility is highly complex. Known genetic factors alone account for at least 15% of male infertility. However, ~40% of infertile men remain undiagnosed and are classified as idiopathic infertile men. We performed exome sequencing in 47 idiopathic infertile men (discovery cohort), followed by replication study (40 variants in 33 genes) in 844 infertile men and 709 controls using Sequenom MassARRAY® based genotyping. We report 17 variants in twelve genes that comprise both previously reported (DNAH8, DNAH17, FISP2 and SPEF2) and novel candidate genes (BRDT, CETN1, CATSPERD, GMCL1, SPATA6, TSSK4, TSKS and ZNF318) for male infertility. The latter have a strong biological nexus to human spermatogenesis and their respective mouse knockouts are concordant with human phenotypes. One candidate gene CETN1, identified in this study, was sequenced in another independent cohort of 840 infertile and 689 fertile men. Further, CETN1 variants were functionally characterized using biophysical and cell biology approaches. We demonstrate that CETN1 variant- p.Met72Thr leads to multipolar cells, fragmented nuclei during mitosis leading to cell death and show significantly perturbed ciliary disassembly dynamics. Whereas CETN1-5' UTR variant; rs367716858 leads to loss of a methylation site and increased reporter gene expression in vitro. We report a total of eight novel candidate genes identified by exome sequencing, which may have diagnostic relevance and can contribute to improved diagnostic workup and clinical management of male infertility.

Signatures of high altitude adaptation in Tibeto-Burman tribes of the Darjeeling Hill Region

OBJECTIVES

The long-term isolation, endogamy practices, and environmental adaptations have shaped the enormous human diversity in India. The genetic and morphological variations in mainland Indians are well studied. However, the data on the Indian Himalayan populations are scattered. Thus, the present study attempts to understand variations in the selected parameter among four Tibeto-Burman speaking ethnic tribal populations from the Darjeeling Hill Region (DHR) in the Eastern Himalaya Biodiversity Hotspot region of India.

METHODS

A total of 178 healthy male individuals (Lepcha 98, Sherpa 31, Bhutia 27, and Tibetans 22) living at an altitudinal range of 1467-2258 m above the sea level were studied for the 10 parameters namely, weight (kg), height (cm), body mass index (BMI) (kg/m² ) systolic and diastolic pressure (mm of Hg), pulse rate (per minute), saturation of peripheral oxygen (SPO₂ ) (%), hemoglobin (g/dl), hematocrit (HCT) (%), and blood glucose (mg/dl). The data was statistically analyzed using analysis of variance and multiple linear regression methods.

RESULTS

Our analysis revealed comparatively lower hemoglobin and HCT levels, and higher systolic and diastolic blood pressure in the Sherpas followed by the Tibetans. This may be reflecting the persistence of high-altitude adaptation signatures even in lowlands. Interestingly, the Tibetans differed significantly from other populations in terms of their higher body weight, height, and BMI.

CONCLUSION

Thus, our study showed the persistence of high altitude signatures in Tibetans and Sherpa inhabited the DHR. Additionally, we also observed significant differences in the anthropometric and physiological parameters among the Tibeto-Burman populations of the DHR.

Genotype-phenotype correlation and natural history study of dysferlinopathy: a single-centre experience from India

Dysferlinopathies are a group of limb-girdle muscular dystrophies causing significant disability in the young population. There is a need for studies on large cohorts to describe the clinical, genotypic and natural history in our subcontinent. To describe and correlate the clinical, genetic profile and natural history of genetically confirmed dysferlinopathies. We analysed a retrospective cohort of patients with dysferlinopathy from a single quaternary care centre in India. A total of 124 patients with dysferlinopathy were included (40 females). Median age at onset and duration of illness were 21 years (range, 13-50) and 48 months (range, 8-288), respectively. The average follow-up period was 60 months (range, 12-288). Fifty-one percent had LGMD pattern of weakness at onset; 23.4% each had Miyoshi and proximo-distal type while isolated hyperCKemia was noted in 1.6%. About 60% were born to consanguineous parents and 26.6% had family history of similar illness. Twenty-three patients (18.6%) lost ambulation at follow-up; the median time to loss of independent ambulation was 120 months (range, 72-264). Single-nucleotide variants (SNVs) constituted 78.2% of patients; INDELs 14.5% and 7.3% had both SNVs and INDELs. Earlier age at onset was noted with SNVs. There was no correlation between the other clinical parameters and ambulatory status with the genotype. Thirty-seven (45.7%) novel pathogenic/likely pathogenic (P/LP) variants were identified out of a total of 81 variations. The c.3191G > A variant was the most recurrent mutation. Our cohort constitutes a clinically and genetically heterogeneous group of dysferlinopathies. There is no significant correlation between the clinico-genetic profile and the ambulatory status.

Generation of a new human induced pluripotent stem cell (hiPSC) line from a South Asian Indian with a MYBPC3Δ25bp variant

Myosin binding protein C3 (MYBPC3) is a thick filament contractile protein that interacts with myosin, titin and actin and regulates cardiac muscle contraction. Genetic variations in the MYBPC3 gene are known causal factors for cardiomyopathy and heart failure. Previously, we identified a recurrent MYBPC3 deletion (25 base pairs) among South Asians associated with cardiomyopathy and heart failure. Here, we generated an induced pluripotent stem cell (iPSC) line using peripheral blood mononuclear cells (PBMC) from an Indian harboring MYBPC3 deletion. This iPSC line displayed embryonic stem cell morphology, expressed pluripotency markers, differentiated into three germ layers and exhibited normal karyotype.

Proteomic and genetic dissection of testis-specific histone 2B in infertile men reveals its contribution to meiosis and sperm motility

OBJECTIVE

To investigate testis-specific histone 2B (TSH2B) and its gene anomalies in infertile men.

DESIGN

Case-control study.

SETTING

Basic science laboratory.

PATIENT(S)

Fertile and infertile men.

INTERVENTION(S)

Not applicable.

MAIN OUTCOME MEASURE(S)

The histone and protamine status of sperm was studied by aniline blue and chromomycin A3 staining, respectively. Testis-specific histone 2B, total H2B, and phosphorylated TSH2B (pTSH2B) were estimated by Western blot analysis. The frequency of genetic polymorphisms and rare variants in H2BC1 was studied by Sanger sequencing. Phosphosites on TSH2B in sperm were identified by reverse-phase high-performance liquid chromatography purification of TSH2B followed by mass spectrometric analysis.

RESULT(S)

Aniline blue and chromomycin A3 staining revealed significantly higher histone retention and low protamine in sperm of infertile men. Sperm TSH2B and total H2B levels were significantly lower in oligozoospermic and oligoasthenozoospermic men (in both groups). The TSH2B levels were comparable in asthenozoospermic men; however, the pTSH2B level was significantly low. The H2BC1 gene sequencing identified 6 variants, of which 2 are rare variants (rs368672899 and rs544942090) and 4 (rs4711096, rs4712959, rs4712960 and rs4712961) are single nucleotide polymorphisms. Minor allele frequency of 5'-untranslated region variant rs4711096 was significantly lower in infertile men (OR = 0.65). The rare nonsynonymous variant, rs368672899, p.Ser5Pro was seen in 1 oligoasthenoteratozoospermic individual. Interestingly, mass spectrometric analysis identified a site on TSH2B to bear a phosphate group in the sperm of fertile men.

CONCLUSION(S)

Our study reveals a defect in the replacement of somatic histones with testis-specific variants in infertile men. Chromatin compaction positively correlates with sperm motility, which is suggestive of its utility in diagnostic semen analysis of infertile individuals. Our observations with TSH2B and its cognate gene in sperm of infertile men indicate an essential role for TSH2B in meiosis and its phosphorylation in sperm motility, respectively.

Ribosomal protein S6 kinase beta-1 gene variants cause hypertrophic cardiomyopathy

BACKGROUND

Hypertrophic cardiomyopathy (HCM) is a genetic heart muscle disease with preserved or increased ejection fraction in the absence of secondary causes. Mutations in the sarcomeric protein-encoding genes predominantly cause HCM. However, relatively little is known about the genetic impact of signalling proteins on HCM.

METHODS AND RESULTS

Here, using exome and targeted sequencing methods, we analysed two independent cohorts comprising 401 Indian patients with HCM and 3521 Indian controls. We identified novel variants in ribosomal protein S6 kinase beta-1 (RPS6KB1 or S6K1) gene in two unrelated Indian families as a potential candidate gene for HCM. The two unrelated HCM families had the same heterozygous missense S6K1 variant (p.G47W). In a replication association study, we identified two S6K1 heterozygotes variants (p.Q49K and p.Y62H) in the UK Biobank cardiomyopathy cohort (n=190) compared with matched controls (n=16 479). These variants are neither detected in region-specific controls nor in the human population genome data. Additionally, we observed an S6K1 variant (p.P445S) in an Arab patient with HCM. Functional consequences were evaluated using representative S6K1 mutated proteins compared with wild type in cellular models. The mutated proteins activated the S6K1 and hyperphosphorylated the rpS6 and ERK1/2 signalling cascades, suggesting a gain-of-function effect.

CONCLUSIONS

Our study demonstrates for the first time that the variants in the S6K1 gene are associated with HCM, and early detection of the S6K1 variant carriers can help to identify family members at risk and subsequent preventive measures. Further screening in patients with HCM with different ethnic populations will establish the specificity and frequency of S6K1 gene variants.

Pinpointing the Geographic Origin of 165-Year-Old Human Skeletal Remains Found in Punjab, India: Evidence From Mitochondrial DNA and Stable Isotope Analysis

In 2014, 157 years after the Sepoy Mutiny of 1857, several unidentified human skeletons were discovered in an abandoned well at Ajnala, Punjab. The most prevailing hypothesis suggested them as Indian soldiers who mutinied during the Indian uprising of 1857. However, there is an intense debate on their geographic affinity. Therefore, to pinpoint their area of origin, we have successfully isolated DNA from cementum-rich material of 50 good-quality random teeth samples and analyzed mtDNA haplogroups. In addition to that, we analyzed 85 individuals for oxygen isotopes (δ¹⁸O values). The mtDNA haplogroup distribution and clustering pattern rejected the local ancestry and indicated their genetic link with the populations living east of Punjab. In addition, the oxygen isotope analysis (δ¹⁸O values) from archaeological skeletal remains corroborated the molecular data and suggested the closest possible geographical affinity of these skeletal remains toward the eastern part of India, largely covering the Gangetic plain region. The data generated from this study are expected to expand our understanding of the ancestry and population affinity of martyr soldiers.

Novel Mutations in β-MYH7 Gene in Indian Patients With Dilated Cardiomyopathy

Background

Methods

Results

Conclusions

COVID-19: Impact on linguistic and genetic isolates of India

The rapid expansion of coronavirus SARS-CoV-2 has impacted various ethnic groups all over the world. The burden of infectious diseases including COVID-19 are generally reported to be higher for the Indigenous people. The historical knowledge have also suggested that the indigenous populations suffer more than the general populations in the pandemic. Recently, it has been reported that the indigenous groups of Brazil have been massively affected by COVID-19. Series of studies have shown that many of the indigenous communities reached at the verge of extinction due to this pandemic. Importantly, South Asia also has several indigenous and smaller communities, that are living in isolation. Till date, despite the two consecutive waves in India, there is no report on the impact of COVID-19 for indigenous tribes. Since smaller populations experiencing drift may have greater risk of such pandemic, we have analysed Runs of Homozygosity (ROH) among South Asian populations and identified several populations with longer homozygous segments. The longer runs of homozygosity at certain genomic regions may increases the susceptibility for COVID-19. Thus, we suggest extreme careful management of this pandemic among isolated populations of South Asia.

Genomic characterization and epidemiology of an emerging SARS-CoV-2 variant in Delhi, India

Delhi, the national capital of India, experienced multiple severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreaks in 2020 and reached population seropositivity of >50% by 2021. During April 2021, the city became overwhelmed by COVID-19 cases and fatalities, as a new variant, B.1.617.2 (Delta), replaced B.1.1.7 (Alpha). A Bayesian model explains the growth advantage of Delta through a combination of increased transmissibility and reduced sensitivity to immune responses generated against earlier variants (median estimates: 1.5-fold greater transmissibility and 20% reduction in sensitivity). Seropositivity of an employee and family cohort increased from 42% to 87.5% between March and July 2021, with 27% reinfections, as judged by increased antibody concentration after a previous decline. The likely high transmissibility and partial evasion of immunity by the Delta variant contributed to an overwhelming surge in Delhi.

Novel FCN2 Variants and Haplotypes are Associated with Rheumatic Heart Disease

Ficolins are pattern recognition molecules that are involved in innate immune defense. Ficonin-2 (FCN2) has a specific affinity for lipoteichoic acid present in the cell wall of Streptococcus pyogenes, an etiological agent for rheumatic heart disease (RHD). We have estimated FCN2 serum levels and analyzed the functional variants of FCN2 in 400 RHD patients, 617 healthy controls, and 581 individuals belonged to various ethnic populations, who are inhabited in various geographical regions of India. Our study revealed that the FCN2 -986A and +6359T alleles were the risk factors for RHD susceptibility (p = 0.0007 for -986G>A; p = 0.0004 for +6359C>T). The haplotype AGGT (p = 0.0024) was observed to be a risk factor for RHD susceptibility, and the haplotype GGAC (p = 0.002) was found to confer protection against RHD. The level of serum FCN2 was significantly higher in controls (p < 0.0001) and in controls with GGAC haplotypes (p < 0.0001). The frequency of the risk alleles -986A and +6359T was found to be more prevalent in Northern and North-Western (Indo-European) India. The protective GGAC haplotype was found more prevalent in Eastern (Tibeto-Burman) and South-Western (Dravidian) India. Alleles -986A and +6359T were in positive correlation with the prevalence of RHD (regression coefficient = 1.84 and 1.94, respectively), whereas GGAC haplotype was in negative correlation with prevalence of RHD (regression coefficient = -1.71). In conclusion, we found that low level of serum ficolin-2 is significantly associated with RHD. Further, FCN2 -986A and +6359T alleles and AGGT haplotype are associated with increased susceptibility to RHD, while GGAC haplotype is associated with moderate protection against RHD.

Advances in mitochondrial medicine and translational research

Current knowledge of mitochondrial biology and function has provided tools and technologies that helped a better understanding of the molecular etiology of complex mitochondrial disorders. Dual genetic control of this subcellular organelle function regulates various signaling mechanisms which are essential for metabolism, bioenergetics, fatty acid biosynthesis, and DNA replication & repair. Understanding nuclear mitochondrial crosstalk through advanced genomics as well as clinical perspectives is the overall basis of mitochondrial research and medicine, also the sole objective of Society for Mitochondrial Medicine and Research (SMRM) - India. The eighth virtual international conference on 'Advances in Mitochondrial Medicine and Translational Research' was organized at the Manipal School of Life Sciences, MAHE, Manipal, India, during 6 - 7 November 2020. The aim of the virtual conference was to highlight the recent advances and future perspectives that represent comprehensive clinical and fundamental research interests in the area of mitochondrial biology of human diseases. To systematically present the various findings in mitochondrial biology, the meeting was themed with specific aspects comprising (a) mitochondrial disorders: clinical & genomic perspectives, (b) mitochondria in cancer, (c) mitochondrial metabolism & disorders, and (d) mitochondrial diseases & therapy. This report provides an overview of the recent advancements in the area of mitochondrial biology and medicine that was discussed at the conference.

Dissecting the genetic history of the Roman Catholic populations of West Coast India

Migration and admixture history of populations have always been curious and an interesting theme. The West Coast of India harbours a rich diversity, bestowing various ethno-linguistic groups, with many of them having well-documented history of migrations. The Roman Catholic is one such distinct group, whose origin was much debated. While some historians and anthropologists relating them to ancient group of Gaud Saraswat Brahmins, others relating them for being members of the Jews Lost Tribes in the first Century migration to India. Historical records suggests that this community was later forcibly converted to Christianity by the Portuguese in Goa during the Sixteenth Century. Till date, no genetic study was done on this group to infer their origin and genetic affinity. Hence, we analysed 110 Roman Catholics from three different locations of West Coast of India including Goa, Kumta and Mangalore using both uniparental and autosomal markers to understand their genetic history. We found that the Roman Catholics have close affinity with the Indo-European linguistic groups, particularly Brahmins. Additionally, we detected genetic signal of Jews in the linkage disequilibrium-based admixture analysis, which was absent in other Indo-European populations, who are inhabited in the same geographical regions. Haplotype-based analysis suggests that the Roman Catholics consist of South Asian-specific ancestry and showed high drift. Ancestry-specific historical population size estimation points to a possible bottleneck around the time of Goan inquisition (fifteenth century). Analysis of the Roman Catholics data along with ancient DNA data of Neolithic and bronze age revealed that the Roman Catholics fits well in a basic model of ancient ancestral composition, typical of most of the Indo-European caste groups of India. Mitochondrial DNA (mtDNA) analysis suggests that most of the Roman Catholics have aboriginal Indian maternal genetic ancestry; while the Y chromosomal DNA analysis indicates high frequency of R1a lineage, which is predominant in groups with higher ancestral North Indian (ANI) component. Therefore, we conclude that the Roman Catholics of Goa, Kumta and Mangalore regions are the remnants of very early lineages of Brahmin community of India, having Indo-Europeans genetic affinity along with cryptic Jewish admixture, which needs to be explored further.

Advances in Mitochondrial Medicine and Translational Research

Current knowledge of mitochondrial biology and function has provided with tools and technologies that helped a better understanding of the molecular etiology of complex mitochondrial disorders. Dual genetic control of this subcellular organelle function regulates various signaling mechanisms which are essential for metabolism, bioenergetics, fatty acid biosynthesis, and DNA replication & repair. Understanding nuclear mitochondrial crosstalk through advanced genomics as well as clinical perspectives is the overall basis of mitochondrial research and medicine, also the sole objective of Society for Mitochondrial Medicine and Research (SMRM) - India. The eighth virtual international conference on 'Advances in Mitochondrial Medicine and Translational Research' was organized at the Manipal School of Life Sciences, MAHE, Manipal, India, during 6 - 7 November 2020. The aim of the virtual conference was to highlight the recent advances and future perspectives that represent comprehensive clinical and fundamental research interests in the area of mitochondrial biology of human diseases. To systematically present the various findings in mitochondrial biology, the meeting was themed with specific aspects comprising (a) mitochondrial disorders: clinical & genomic perspectives, (b) mitochondria in cancer, (c) mitochondrial metabolism & disorders, and (d) mitochondrial diseases & therapy. This report provides an overview of the recent advancements in the area of mitochondrial biology and medicine that was discussed at the conference.

Conformational scanning of individual EF-hand motifs of calcium sensor protein centrin-1

Centrin-1, a Ca²⁺ sensor protein of the centrin family is a crucial player for cell division in eukaryotes and plays a key role in the microtubule organising centre. Despite being regarded as a calcium sensor with a matched structure to calmodulin/troponin C, the protein undergoes mild changes in conformation and binds Ca²⁺ with moderate affinity. We present an in-depth analysis of the Ca²⁺ sensing by individual EF-hand motifs of centrin-1 and address unsolved questions of the rationales for moderate affinity and conformational transitions of the protein. Employing the more sensitive approach of Trp scanning of individual EF-hand motif, we have undertaken an exhaustive investigation of Ca²⁺ binding to individual EF-hand motifs, named EF1 to EF4. All four EF-hand motifs of centrin-1 are structural as all of them bind both Ca²⁺ and Mg²⁺. EF1 and EF4 are the most flexible sites as they undergo drastic conformational changes following Ca²⁺ binding, whereas EF3 responds to Ca²⁺ minimally. On the other hand, EF2 moves towards the protein surface upon binding Ca²⁺. The independent filling mode of Ca²⁺ to EF-hand motifs and lack of intermotif communication explain the lack of cooperativity of binding, thus constraining centrin-1 to a moderate affinity binding protein. Thus, centrin-1 is distinct from other calcium sensors such as calmodulin.

Heterogeneity in Hematological Parameters of High and Low Altitude Tibetan Populations

Introduction

High altitude hypoxia is believed to be experienced at elevations of more than 2500 meters above sea level. Several studies have shed light on the biochemical aspects of high altitude acclimatization, where participants were sojourners to the high altitude from low altitude areas. However, information regarding the difference between the high altitude adapted Tibetans living at high altitude and their counterparts who reside at low altitude are lacking. To understand this, we have measured various hematological parameters in the Tibetan populations, who are residing in both high and low altitudes in India.

Methods

A total of 168 individuals (79 from high altitude (≥4500 meters) and 89 from low altitude (~850 meters) were recruited for this study. Hematological parameters such as red blood cells (RBC) count, hematocrit (HCT), hemoglobin concentration (Hb), mean corpuscular volume (MCV), mean corpuscular hemoglobin (MCH) and mean corpuscular hemoglobin concentration (MCHC) were measured from the individuals from high and low altitudes. Serum erythropoietin (EPO) was measured by ELISA. Statistical analyses were performed to compare data from both of the altitudes. Gender-wise comparison of data was reported. Correlation analysis was performed within relevant parameters.

Results

Highly significant differences (p <0.0001) between high and low altitude Tibetans were detected in RBC count, HCT, Hb, MCHC in both males and females and in MCV in females. In the case of MCHC, however, age and BMI were potential confounders. Nominally significant differences (p <0.05) were detected in MCV and MCH within males. No significant difference in serum EPO level was found between altitude groups, in any gender. No significant correlation was found between serum EPO with Hb as well as serum EPO with HCT.

Discussion

Our study explores significantly lower RBC count, HCT, Hb, MCH, MCHC and higher MCV in long-term Tibetan residents living at low altitude compared to their high altitude counterparts, which is likely due to the outcome of hematological adaptation to a relatively hyperoxic environment in low altitude areas.

Mutations in the desert hedgehog (DHH) gene in the disorders of sexual differentiation and male infertility

PURPOSE

To identify the contribution of mutations in the Desert Hedgehog (DHH) gene to the disorders of sexual differentiation (DSD) and male infertility.

METHODS

The study included a total 430 subjects, including 47 gonadal dysgenesis cases, 6 patients with undescended testis and infertility characterized by azoospermia, 125 infertile male patients characterized by oligoasthenozoospermia, 24 patients with oligoasthenoteratozoospermia, and 200 ethnically matched normozoospermic fertile men who had fathered a child in the last two years. Sequencing of the complete coding region of the DHH gene was undertaken to find its contribution to the DSD and male infertility.

RESULTS

We observed four novel mutations in the DHH gene in the cases with different reproductive anomalies. A synonymous substitution, c. 543C>T (p.His181His) was observed in 6.6% oligoasthenozoospermic infertile males and 1.5% normozoospermic fertile control samples (RR = 4.4077, 95%CI 1.19-16.29). Another synonymous substitution, c.990G>A (p.Ala330Ala) was observed in an infertile patient with unilateral undescended testis (case #12). Insertion of G at c.1156insG (p.Arg385fs) was observed in a case with bilateral undescended testis and azoospermia (case #23). In gonadal dysgenesis category, two mutations, insertion of G at c.1156insG (p.Arg385fs) and c.997A>G (p.Thr333Ala) substitution were observed in one case (case #34). These mutations were completely absent in control samples.

CONCLUSION

Mutations in the DHH gene impact reproduction with mild mutations affecting fertility, and severe or multiple mutations resulting in gonadal dysgenesis.

CYP2C9 Variations and Their Pharmacogenetic Implications Among Diverse South Asian Populations

Introduction

Allelic frequency distribution of drug metabolizing enzyme genes among populations is important to identify risk groups for adverse drug reaction and to select representative populations for clinical trials. Although India emerged as an important hub for clinical trials, information about the pharmacogenetic diversity for this region is still lacking. Here, we investigated genetic diversity of cytochrome-P450-2C9 (CYP2C9) gene which metabolizes wide range of drugs and is highly expressed in the human liver.

Methods

In total, 1278 individuals from 36 diverse Indian populations, 210 individuals from in-house data-repository and 489 other South Asian samples from the 1000 Genomes Project were selected. Variants observed in CYP2C9 gene were subjected to various statistical analyses.

Results

High frequency of CYP2C9*3 (~13%) and CYP2C9*3/*3 (~1%) was observed among South Asians, compared to 21 populations living outside the Indian subcontinent. The allelic/genotypic frequency does not correlate with geographical location or linguistic affiliation, except populations speaking Tibeto-Burmans language, who have lower frequency of CYP2C9*3 and CYP2C9*3/*3. Since, South Asians practice strict endogamy, presence of unique mutation and high frequency of homozygous genotypes not surprising. CYP2C9*3 has been associated with therapeutic response.The effect of CYP2C9*3/*3 is more pronounced compared to heterozygous and wild type homozygous genotypes as evident in many in vitro studies. As South Asians have high frequency, it would be interesting to explore potential of CYP2C9*3 as a marker for personalized therapy. Our study revealed several rare functional variants, which form eight novel and rare haplotypes of CYP2C9 (CYP2C9*63–*70). Of which, CYP2C9*64, *65, *66, *68, *69 and *70 haplotypes are South Asian-specific.

Conclusion

Overall, we find high genetic heterogeneity within South Asians and identified South Asian-specific putative functional CYP2C9 haplotypes. High frequency of CYP2C9*3 and CYP2C9*3/*3 was observed in South Asian populations. Taken together, current study greatly enriches the knowledge of naturally occurring CYP2C9 variants and its diversity in South Asia, which are relevant to further CYP2C9-related functional research and for personalized medicine.

Contribution of nuclear and mitochondrial gene mutations in mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) syndrome

Background

Mitochondrial disorders are clinically complex and have highly variable phenotypes among all inherited disorders. Mutations in mitochon drial DNA (mtDNA) and nuclear genome or both have been reported in mitochondrial diseases suggesting common pathophysiological pathways. Considering the clinical heterogeneity of mitochondrial encephalopathy, lactic acidosis and stroke-like episodes (MELAS) phenotype including focal neurological deficits, it is important to look beyond mitochondrial gene mutation.

Methods

The clinical, histopathological, biochemical analysis for OXPHOS enzyme activity, and electron microscopic, and neuroimaging analysis was performed to diagnose 11 patients with MELAS syndrome with a multisystem presentation. In addition, whole exome sequencing (WES) and whole mitochondrial genome sequencing were performed to identify nuclear and mitochondrial mutations.

Results

Analysis of whole mtDNA sequence identified classical pathogenic mutation m.3243A > G in seven out of 11 patients. Exome sequencing identified pathogenic mutation in several nuclear genes associated with mitochondrial encephalopathy, sensorineural hearing loss, diabetes, epilepsy, seizure and cardiomyopathy (POLG, DGUOK, SUCLG2, TRNT1, LOXHD1, KCNQ1, KCNQ2, NEUROD1, MYH7) that may contribute to classical mitochondrial disease phenotype alone or in combination with m.3243A > G mutation.

Conclusion

Individuals with MELAS exhibit clinical phenotypes with varying degree of severity affecting multiple systems including auditory, visual, cardiovascular, endocrine, and nervous system. This is the first report to show that nuclear genetic factors influence the clinical outcomes/manifestations of MELAS subjects alone or in combination with m.3243A > G mutation.

Electronic supplementary material

The online version of this article (10.1007/s00415-020-10390-9) contains supplementary material, which is available to authorized users.

Adiponectin receptor 1 variants contribute to hypertrophic cardiomyopathy that can be reversed by rapamycin

Hypertrophic cardiomyopathy (HCM) is a heterogeneous genetic heart muscle disease characterized by hypertrophy with preserved or increased ejection fraction in the absence of secondary causes. However, recent studies have demonstrated that a substantial proportion of individuals with HCM also have comorbid diabetes mellitus (~10%). Whether genetic variants may contribute a combined phenotype of HCM and diabetes mellitus is not known. Here, using next-generation sequencing methods, we identified novel and ultrarare variants in adiponectin receptor 1 (ADIPOR1) as risk factors for HCM. Biochemical studies showed that ADIPOR1 variants dysregulate glucose and lipid metabolism and cause cardiac hypertrophy through the p38/mammalian target of rapamycin and/or extracellular signal-regulated kinase pathways. A transgenic mouse model expressing an ADIPOR1 variant displayed cardiomyopathy that recapitulated the cellular findings, and these features were rescued by rapamycin. Our results provide the first evidence that ADIPOR1 variants can cause HCM and provide new insights into ADIPOR1 regulation.

High-altitude adaptation: Role of genetic and epigenetic factors

After the completion of the Human Genome Project in 2003, the field of genetics has witnessed massive progress that spanned research in high-altitude biology also. Especially the decade of 2010s witnessed the most of it and revealed various genetic signatures of high-altitude adaptation in Tibetans, Andeans and Ethiopians. High-altitude area, with its extreme environment, harbors a tremendous potential for gene-environment interaction, an aspect that could be explored by epigenetic studies. There are only four original articles till now which explore the epigenetic aspect of high-altitude adaptation or acclimatization. However, there is no comprehensive review to provide complete information on the genetic and epigenetic aspects of high-altitude adaptations. Hence, we have prepared this mini-review to summarize the genetic and epigenetic studies that have correlated the high-altitude adaptation or acclimatization, until recently.

A Distinct Phylogenetic Cluster of Indian Severe Acute Respiratory Syndrome Coronavirus 2 Isolates

BACKGROUND

From an isolated epidemic, coronavirus disease 2019 has now emerged as a global pandemic. The availability of genomes in the public domain after the epidemic provides a unique opportunity to understand the evolution and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus across the globe.

METHODS

We performed whole-genome sequencing of 303 Indian isolates, and we analyzed them in the context of publicly available data from India.

RESULTS

We describe a distinct phylogenetic cluster (Clade I/A3i) of SARS-CoV-2 genomes from India, which encompasses 22% of all genomes deposited in the public domain from India. Globally, approximately 2% of genomes, which to date could not be mapped to any distinct known cluster, fall within this clade.

CONCLUSIONS

The cluster is characterized by a core set of 4 genetic variants and has a nucleotide substitution rate of 1.1 × 10-3 variants per site per year, which is lower than the prevalent A2a cluster. Epidemiological assessments suggest that the common ancestor emerged at the end of January 2020 and possibly resulted in an outbreak followed by countrywide spread. To the best of our knowledge, this is the first comprehensive study characterizing this cluster of SARS-CoV-2 in India.

Novel Biallelic NSUN3 Variants Cause Early-Onset Mitochondrial Encephalomyopathy and Seizures

Epitranscriptomic systems enable post-transcriptional modifications of cellular RNA that are essential for regulating gene expression. Of the ~ 170 known RNA chemical modifications, methylation is among the most common. Loss of function mutations in NSUN3, encoding the 5-methylcytosine (m5C) methyltransferase NSun3, have been linked to multisystem mitochondrial disease associated with combined oxidative phosphorylation deficiency. Here, we report a patient with early-onset mitochondrial encephalomyopathy and seizures in whom the novel biallelic NSUN3 missense variants c.421G>C (p.A141P) and c.454T>A (p.C152S) were detected. Segregation studies and in silico functional analysis confirmed the likely pathogenic effects of both variants. These findings expand the molecular and phenotypic spectrum of NSUN3-related mitochondrial disease.

Mitochondrial diversity of Yoruba and Fulani chickens: A biodiversity reservoir in Nigeria

Poultry are the most widely distributed type of livestock in Nigeria. Indigenous chickens are extremely common throughout the country. Indeed, approximately 83 million chickens are raised in extensive systems and 60 million in semi-intensive systems. To provide the first comprehensive overview of the maternal lineages in Southwest Nigeria, we analyzed 96 mitochondrial DNA control region sequences from 2 indigenous chicken ecotypes: Fulani and Yoruba. All samples belonged to the most frequent haplogroup (E) in Africa and Europe and showed noticeably low haplotype diversity. Although only 11 different haplotypes were detected, with 2 of them never found before in Nigeria, the presence of unique sequences among our indigenous samples testified to their status as an important genetic resource to be preserved. Furthermore, a total of 7,868 published sequences were included in the comparative analysis, which revealed an east-west geographic pattern of haplogroup distribution and led to the conclusion that the gene flow from Southeastern Asia mainly involved one mitochondrial clade. Moreover, owing to the extensive genetic intermixing among Nigerian chickens, conservation efforts are required to safeguard the extant mitochondrial variability in these indigenous ecotypes and establish future improvement and selection programs.

Ancient mtDNA from the extinct Indian cheetah supports unexpectedly deep divergence from African cheetahs

The Indian cheetah was hunted to extinction by the mid-20th century. While analysis of 139 bp of mitochondrial DNA (mtDNA) has confirmed that the Indian cheetah was part of the Asiatic subspecies (Acinonyx jubatus venaticus), the detailed relationships between cheetah populations remains unclear due to limited genetic data. We clarify these relationships by studying larger fragments of cheetah mtDNA, both from an Indian cheetah museum specimen and two African cheetah, one modern and one historic, imported into India at different times. Our results suggest that the most recent common ancestor of cheetah mtDNA is approximately twice as ancient as currently recognised. The Indian and Southeast African (Acinonyx jubatus jubatus) cheetah mtDNA diverged approximately 72 kya, while the Southeast and Northeast African (Acinonyx jubatus soemmeringii) cheetah mtDNA diverged around 139 kya. Additionally, the historic African cheetah sampled from India proved to have an A. j. jubatus haplotype, suggesting a hitherto unrecognised South African route of cheetah importation into India in the 19^th century. Together, our results provide a deeper understanding of the relationships between cheetah subspecies, and have important implications for the conservation of A. j. venaticus and potential reintroduction of cheetahs into India.

Novel Homozygous FAN1 Mutation in a Familial Case of Karyomegalic Interstitial Nephritis

Karyomegalic interstitial nephritis (KIN) is a rare genetic kidney disease associated with a mutation in FAN1 gene and is often underdiagnosed. The histomorphology demonstrates chronic interstitial nephritis with tubular epithelial cells showing bizarre enlarged nuclei. We present a case report of a 47-year-old multiparous South-Indian woman presenting with bilateral pitting pedal oedema and mild hypertension. At the time of presentation, her serum creatinine was 1.52 mg/dL and urine analysis showed mild proteinuria. Kidney biopsy showed features of tubular injury with bizarre enlarged nuclei and focal mild chronic tubulointerstitial nephritis. Immunohistochemistry was negative for cytomegalovirus (CMV) Ag and SV40 Ag. Real-time polymerase chain reaction (PCR) done for CMV and BK virus genomes was negative. Relevant family history was that her older brother was also diagnosed with kidney failure and is on renal replacement therapy. Genetic analysis for FAN1 gene of the proband and her sibling showed two rare mutations of the FAN1 gene in the exon 4, of which, one is non-synonymous mutation and the other is a stop-gain mutation in the proband. This case illustrates a rare presentation of karyomegalic interstitial nephritis in siblings with previous unknown FAN1 gene mutations.

Novel NR5A1 Pathogenic Variants Cause Phenotypic Heterogeneity in 46,XY Disorders of Sex Development

Steroidogenic factor 1 (NR5A1/SF1) is a key transcription factor that is known to regulate the development of adrenal glands and gonads and is also involved in steroidogenesis. Several pathogenic NR5A1 variants have been reported to cause 46,XY disorders of sex development (DSD), with varying clinical phenotypes ranging from hypospadias to complete gonadal dysgenesis. Most often, the primary cause of DSD is due to variants in gene(s) related to gonadal development or the steroidogenic pathway. In the present study, we have analyzed 64 cases of 46,XY DSD for pathogenic NR5A1 variants. We report a total of 3 pathogenic variants of which 2 were novel (p.Gly22Ser and p.Ser143Asn) and 1 was already known (p.Ser32Asn). Functional studies have revealed that the 2 mutations p.Gly22Ser and p.Ser32Asn could significantly affect DNA binding and transactivation abilities. Further, these mutant proteins showed nuclear localization with aggregate formation. The third mutation, p.Ser143Asn, showed unspeckled nuclear localization and normal DNA binding, but the ability of transcriptional activation was significantly reduced. In conclusion, we recommend screening for NR5A1 pathogenic variants in individuals with features of 46,XY DSD for better diagnosis and management.

The formation of human populations in South and Central Asia

By sequencing 523 ancient humans, we show that the primary source of ancestry in modern South Asians is a prehistoric genetic gradient between people related to early hunter-gatherers of Iran and Southeast Asia. After the Indus Valley Civilization's decline, its people mixed with individuals in the southeast to form one of the two main ancestral populations of South Asia, whose direct descendants live in southern India. Simultaneously, they mixed with descendants of Steppe pastoralists who, starting around 4000 years ago, spread via Central Asia to form the other main ancestral population. The Steppe ancestry in South Asia has the same profile as that in Bronze Age Eastern Europe, tracking a movement of people that affected both regions and that likely spread the distinctive features shared between Indo-Iranian and Balto-Slavic languages.

An Ancient Harappan Genome Lacks Ancestry from Steppe Pastoralists or Iranian Farmers

We report an ancient genome from the Indus Valley Civilization (IVC). The individual we sequenced fits as a mixture of people related to ancient Iranians (the largest component) and Southeast Asian hunter-gatherers, a unique profile that matches ancient DNA from 11 genetic outliers from sites in Iran and Turkmenistan in cultural communication with the IVC. These individuals had little if any Steppe pastoralist-derived ancestry, showing that it was not ubiquitous in northwest South Asia during the IVC as it is today. The Iranian-related ancestry in the IVC derives from a lineage leading to early Iranian farmers, herders, and hunter-gatherers before their ancestors separated, contradicting the hypothesis that the shared ancestry between early Iranians and South Asians reflects a large-scale spread of western Iranian farmers east. Instead, sampled ancient genomes from the Iranian plateau and IVC descend from different groups of hunter-gatherers who began farming without being connected by substantial movement of people.

Ancient DNA from the skeletons of Roopkund Lake reveals Mediterranean migrants in India

Situated at over 5,000 meters above sea level in the Himalayan Mountains, Roopkund Lake is home to the scattered skeletal remains of several hundred individuals of unknown origin. We report genome-wide ancient DNA for 38 skeletons from Roopkund Lake, and find that they cluster into three distinct groups. A group of 23 individuals have ancestry that falls within the range of variation of present-day South Asians. A further 14 have ancestry typical of the eastern Mediterranean. We also identify one individual with Southeast Asian-related ancestry. Radiocarbon dating indicates that these remains were not deposited simultaneously. Instead, all of the individuals with South Asian-related ancestry date to ~800 CE (but with evidence of being deposited in more than one event), while all other individuals date to ~1800 CE. These differences are also reflected in stable isotope measurements, which reveal a distinct dietary profile for the two main groups.

Remains of several hundred humans are scattered around Roopkund Lake, situated over 5,000 meters above sea level in the Himalayan Mountains. Here the authors analyze genome-wide data from 38 skeletons and find 3 clusters with different ancestries and dates, showing the people were desposited in multiple catastrophic events.

Peopling of India: Ancient DNA perspectives

To reconstruct and explain patterns of genetic diversity of modern humans, understanding their past and present genetic profile is crucial. While genomes of contemporary people can provide information about present day population structure, analysis of ancient genomes may provide unprecedented insights about the past demographic events that have shaped the contemporary gene pool. Population genetics has recently witnessed an explosion in studies on ancient human population histories, primarily from Europe and America. South Asia has no representation in the ancient genomics literature, despite the wealth of archaeological richness in the form of human skeletal remains that exist in collections all over the country. Representing one-fifth of present day humanity calls for understanding the demographic history of south Asia not merely as a prerequisite but as an urgent need to understand its genetic variations on a global scale. Although the overall picture is taking form, new archaeological and genetic information from the region has started to reveal a more complex scenario of ancient human migrations and admixtures than was ever known before. In this article, we discuss a meaningful insight on the current status of ancient DNA (aDNA) research in India. We have also summarized a few but important aDNA studies, which have been successfully carried out in India. Furthermore, we have highlighted the potential opportunity of aDNA research in the Indian subcontinent.

Retinoblastoma discordance in families with twins

Retinoblastoma has an increased inheritance risk of germline RB1 mutations in offspring and siblings, especially twins. Three families, each having one retinoblastoma-affected twin, were selected for genetic analysis and DNA profiling. Germline RB1 mutations were found in all probands. DNA profiling carried on similar-looking twins of families I and II, proved them to be fraternal. This study demonstrates the importance of genetic analysis of RB1 gene for risk prediction in retinoblastoma families. It also emphasizes that DNA profiling is a mandate for genetic screening of families with twins, thus adding a new dimension in counseling of retinoblastoma.

The peopling of Lakshadweep Archipelago

The archipelago of Lakshadweep is considered as a stopover to the maritime route since ancient time. It is not very clear when the human first occupied these islands, however in the long history of the islands, the local legends suggest that Lakshadweep has been ruled by different kingdoms. To have a better understanding of peopling of Lakshadweep, we have analysed 557 individuals from eight major islands for mitochondrial DNA and 166 individuals for Y chromosome markers. We found a strong founder effect for both paternal and maternal lineages. Moreover, we report a close genetic link of Lakshadweep islanders with the Maldives, Sri Lanka and India. Most of the Lakshadweep islands share the haplogroups specific to South Asia and West Eurasia, except Minicoy Island that also shares haplogroups of East Eurasia. The paternal and maternal ancestries of the majority of island populations suggest their arrival from distinct sources. We found that the maternal ancestry was closer to South Indian populations, whereas the paternal ancestry was overwhelmed with the haplogroups, more common in the Maldives and North of India. In conclusion, our first genetic data suggest that the majority of human ancestry in Lakshadweep is largely derived from South Asia with minor influences from East and West Eurasia.

High frequencies of Non Allelic Homologous Recombination (NAHR) events at the AZF loci and male infertility risk in Indian men

Deletions in the AZoospermia Factor (AZF) regions (spermatogenesis loci) on the human Y chromosome are reported as one of the most common causes of severe testiculopathy and spermatogenic defects leading to male infertility, yet not much data is available for Indian infertile men. Therefore, we screened for AZF region deletions in 973 infertile men consisting of 771 azoospermia, 105 oligozoospermia and 97 oligoteratozoospermia cases, along with 587 fertile normozoospermic men. The deletion screening was carried out using AZF-specific markers: STSs (Sequence Tagged Sites), SNVs (Single Nucleotide Variations), PCR-RFLP (Polymerase Chain Reaction - Restriction Fragment Length Polymorphism) analysis of STS amplicons, DNA sequencing and Southern hybridization techniques. Our study revealed deletion events in a total of 29.4% of infertile Indian men. Of these, non-allelic homologous recombination (NAHR) events accounted for 25.8%, which included 3.5% AZFb deletions, 2.3% AZFbc deletions, 6.9% complete AZFc deletions, and 13.1% partial AZFc deletions. We observed 3.2% AZFa deletions and a rare long AZFabc region deletion in 0.5% azoospermic men. This study illustrates how the ethnicity, endogamy and long-time geographical isolation of Indian populations might have played a major role in the high frequencies of deletion events.

Mitochondrial genome variations in idiopathic dilated cardiomyopathy

Idiopathic dilated cardiomyopathy (DCM) is a structural heart disease with strong genetic background. The aim of this study was to assess the role of mitochondrial DNA (mtDNA) variations and haplogroups in Indian DCM patients. Whole mtDNA analysis of 221 DCM patients revealed 48 novel, 42 disease-associated and 97 private variations. The frequency of reported variations associated with hearing impairment, DEAF, SNHL and LHON are significantly high in DCM patients than controls. Haplogroups H and HV were over represented in DCM than controls. Functional analysis of two private variations (m.8812A>G & m.10320G>A) showed decrease in mitochondrial functions, suggesting the role of mtDNA variations in DCM.

A Complete Absence of Missense Mutation in Myosin Regulatory and Essential Light Chain Genes of South Indian Hypertrophic and Dilated Cardiomyopathies

BACKGROUND

Myosin is a hexameric contractile protein composed of 2 heavy chains associated with 4 light chains of 2 distinct classes - 2 regulatory light chains (MYL2) and 2 essential light chains (MYL3). The myosin light chains stabilize the long alpha helical neck of the myosin head and regulate the myosin ATPase activities.

OBJECTIVES

Mutations in MYL2 and MYL3 are reported to be associated with cardiomyopathies. However, there is no study available on these genes in Indian cardiomyopathies, and therefore we planned to study them.

METHOD

For the first time we sequenced MYL2 and MYL3 genes in a total of 248 clinically well-characterized cardiomyopathies consisting of 101 hypertrophic and 147 dilated cases along with 207 healthy controls from south India.

RESULTS

Our study revealed a total of 10 variations - 7 in MYL2 and 3 in MYL3, of which 3 are novel variations observed exclusively in cases. However, the 15 causative missense mutations previously reported are totally absent in our study, which showed that the sequences of MYL2 and MYL3 are highly conserved in Indian cases/controls.

CONCLUSIONS

MYL2 and MYL3 mutations are rare and the least cause of cardiomyopathies in Indians.

Variations in macrophage migration inhibitory factor gene are not associated with visceral leishmaniasis in India

BACKGROUND

The host genetic factors play important role in determining the outcome of visceral leishmaniasis (VL). Macrophage migration inhibitory factor (MIF) is an important host cytokine, which is a key regulator of innate immune system. Genetic variants in MIF gene have been found to be associated with several inflammatory and infectious diseases. Role of MIF is well documented in leishmaniasis diseases, including Indian visceral leishmaniasis, where elevated level of serum MIF has been associated with VL phenotypes. However, there was no genetic study to correlate MIF variants in VL, therefore, we aimed to study the possible association of three reported MIF gene variants -794 CATT, -173G > C and non-coding RNA gene LOC284889 in Indian VL phenotype.

METHODS

Study subjects comprised of 214 VL patients along with ethnically and demographically matched 220 controls from VL endemic regions of Bihar state in India.

RESULTS

We found no significant difference between cases and controls in allelic, genotypic and haplotype frequency of the markers analysed [-794 CATT repeats (χ2=0.86; p=0.35; OR=0.85; 95% CI=0.61-1.19); -173 G>C polymorphism (χ2=1.11; p=0.29; OR=0.83; 95% CI=0.59-1.16); and LOC284889 (χ2=0.78; p=0.37; OR=0.86; 95% CI=0.61-1.20)].

CONCLUSION

Since we did not find any significant differences between case and control groups, we conclude that sequencing of complete MIF gene and extensive study on innate and adaptive immunity genes may help in identifying genetic variations that are associated with VL susceptibility/resistance among Indians.

SNPs in ERCC1, ERCC2, and XRCC1 genes of the DNA repair pathway and risk of male infertility in the Asian populations: association study, meta-analysis, and trial sequential analysis

PURPOSE

We investigated if substitutions in the ERCC1, ERCC2, and XRCC1 genes of the DNA repair pathway correlate with non-obstructive azoospermia and male infertility.

METHODS

A total of 548 azoospermic infertile males and 410 fertile controls were genotyped for XRCC1 399A > G, 280G > A, and ERCC1 C > A 3' UTR and 541 azoospermic infertile males and 416 fertile controls were genotyped for ERCC2 751A > C using iPLEX Gold Assay. Meta-analyses were performed on XRCC1 399A > G (1022 cases and 1004 controls), ERCC1 C > A 3' UTR (879 cases and 1059 controls), and ERCC2 751A > C (914 cases and 850 controls) polymorphisms to quantitatively estimate the significance of the association between these polymorphisms and the risk of infertility.

RESULTS

Statistically significant association between ERCC2 751A > C SNP and male infertility was found using the codominant model (p = 0.03). Results of meta-analysis suggested a lack of correlation with male infertility risk, which could be due to pooling of studies from different ethnic populations. Due to limited the number of studies, a stratified analysis for different ethnic groups could not be performed.

CONCLUSION (S)

In conclusion, AA genotype of 751A > C SNP in ERCC2 correlated with a higher risk of male infertility and may contribute to an increased risk of azoospermia and male infertility in Indian men.