Metagenome of Mehsani buffalo rumen microbiota: an assessment of variation in feed-dependent phylogenetic and functional classification

AIM

To gain a greater understanding of the ecology and metabolic potential of the rumen microbiome with the changes in the animal diet.

METHODS

Diet composed of varying proportion of green and dry roughages along with grains was given to 8 Mehsani buffaloes, and rumen metagenome was sketched using shotgun semiconductor sequencing.

RESULTS

In the present study, the Bacteroidetes were found to be dominant at the phyla level and Prevotella at the genus level. The ratio of Firmicutes to Bacteroidetes was found to be higher in the solid fraction as compared to the liquid fraction. In the solid fraction of the dry roughage group, the significant increment (p < 0.05) in Bacteroidetes abundance was observed with increment of roughage concentration. At the genus level, Clostridium significantly increased with the increment in roughage concentration. A comparison of glycoside hydrolase and cellulosome functional genes revealed more glycoside hydrolase 3 encoding genes with higher fiber diet and significant difference in carbohydrate-active enzymes family composition between green and dry roughage groups of the liquid fraction.

CONCLUSION

The present study provides a base to understand the modulating behavior of microbiota which can be manipulated to improve livestock nutrient utilization efficiency and for targeting the efficient catabolism of complex carbohydrate molecules as well.

Evaluation of genetic diversity and population structure of West-Central Indian cattle breeds

Evaluations of genetic diversity in domestic livestock populations are necessary to implement region-specific conservation measures. We determined the genetic diversity and evolutionary relationships among eight geographically and phenotypically diverse cattle breeds indigenous to west-central India by genotyping these animals for 22 microsatellite loci. A total of 326 alleles were detected, and the expected heterozygosity ranged from 0.614 (Kenkatha) to 0.701 (Dangi). The mean number of alleles among the cattle breeds ranged from 7.182 (Khillar) to 9.409 (Gaolao). There were abundant genetic variations displayed within breeds, and the genetic differentiation was also high between the Indian cattle breeds, which displayed 15.9% of the total genetic differentiation among the different breeds. The genetic differentiation (pairwise FST ) among the eight Indian breeds varied from 0.0126 for the Kankrej-Malvi pair to 0.2667 for Khillar-Kenkatha pair. The phylogeny, principal components analysis, and structure analysis further supported close grouping of Kankrej, Malvi, Nimari and Gir; Gaolao and Kenkatha, whereas Dangi and Khillar remained at distance from other breeds.

Parentage verification in field progeny testing program of Mehsana buffalo

The present study was undertaken to construct a multiplex microsatellite panel for parentage testing in Mehsana buffalo (Bubalus bubalis). The study was based on a total of 212 Mehsana buffalos (100 dams, 100 daughters, and 12 sires). Genomic DNA was extracted from blood and semen samples. A panel of 10 microsatellite markers (CSSM61, ILSTS29, ILSTS17, ILSTS28, CSSM57, CSSM22, ILSTS61, CSSM8, ETH152, and ILSTS11) was amplified in a single multiplex reaction and analyzed by capillary electrophoresis on an automated DNA sequencer. The expected heterozygosity ranged from 0.642 to 0.833 (mean 0.762). The total exclusion probability using 10 microsatellite loci with 1 known parent was 0.993. Seven out of 10 microsatellite loci revealed relatively high polymorphic information content (>0.7). Eighty-one daughters out of 100 daughters qualified by compatibility according to Mendelism. The results suggest that multiplex microsatellite panel is a fast, robust, reliable, and economic tool to verify the parentage as well as to assign the putative sire to daughters under progeny testing with very high accuracy and hence can be used in routine parentage testing.

Split-hand/split-foot malformation 3 (SHFM3) at 10q24, development of rapid diagnostic methods and gene expression from the region

Split-hand/split-foot malformation (SHFM, also called ectrodactyly) is a clinically variable and genetically heterogeneous group of limb malformations. Several SHFM loci have been mapped, including SHFM1 (7q21), SHFM2 (Xq26), SHFM3 (10q24), SHFM4 (3q27) and SHFM5 (2q31). To date, mutations in a gene (TP63) have only been identified for SHFM4. SHFM3 has been shown by pulsed-field gel electrophoresis to be caused by an approximately 500 kb DNA rearrangement at 10q24. This region contains a number of candidate genes for SHFM3, though which gene(s) is (are) involved in the pathogenesis of SHFM3 is not known. Our aim in this study was to improve the diagnosis of SHFM3, and to begin to understand which genes are involved in SHFM3. Here we show, using two different techniques, FISH and quantitative PCR that SHFM3 is caused by a minimal 325 kb duplication containing only two genes (BTRC and POLL). The data presented provide improved methods for diagnosis and begin to elucidate the pathogenic mechanism of SHFM3. Expression analysis of 13 candidate genes within and flanking the duplicated region shows that BTRC (present in three copies) and SUFU (present in two copies) are overexpressed in SHFM3 patients compared to controls. Our data suggest that SHFM3 may be caused by overexpression of BTRC and SUFU, both of which are involved in beta-catenin signalling.

Genomewide scan for nonsyndromic cleft lip and palate in multigenerational Indian families reveals significant evidence of linkage at 13q33.1-34

Nonsyndromic cleft lip with or without cleft palate (CL-P) is a common congenital anomaly with incidence ranging from 1 in 300 to 1 in 2,500 live births. We analyzed two Indian pedigrees (UR017 and UR019) with isolated, nonsyndromic CL-P, in which the anomaly segregates as an autosomal dominant trait. The phenotype was variable, ranging from unilateral to bilateral CL-P. A genomewide linkage scan that used approximately 10,000 SNPs was performed. Nonparametric linkage (NPL) analysis identified 11 genomic regions (NPL>3.5; P<.005) that could potentially harbor CL-P susceptibility variations. Among those, the most significant evidence was for chromosome 13q33.1-34 at marker rs1830756 (NPL=5.57; P=.00024). This was also supported by parametric linkage; MOD score (LOD scores maximized over genetic model parameters) analysis favored an autosomal dominant model. The maximum LOD score was 4.45, and heterogeneity LOD was 4.45 (alpha =100%). Haplotype analysis with informative crossovers enabled the mapping of the CL-P locus to a region of approximately 20.17 cM (7.42 Mb) between SNPs rs951095 and rs726455. Thus, we have identified a novel genomic region on 13q33.1-34 that harbors a high-risk variant for CL-P in these Indian families.

The phenotypic spectrum of GLI3 morphopathies includes autosomal dominant preaxial polydactyly type-IV and postaxial polydactyly type-A/B; No phenotype prediction from the position of GLI3 mutations

Functional characterization of a gene often requires the discovery of the full spectrum of its associated phenotypes. Mutations in the human GLI3 gene have been identified in Greig cepalopolysyndactyly, Pallister-Hall syndrome (PHS), and postaxial polydactyly type-A (PAP-A). We studied the involvement of GLI3 in additional phenotypes of digital abnormalities in one family (UR003) with preaxial polydactyly type-IV (PPD-IV), three families (UR014, UR015, and UR016) with dominant PAP-A/B (with PPD-A and -B in the same family), and one family with PHS. Linkage analysis showed no recombination with GLI3-linked polymorphisms. Family UR003 had a 1-nt frameshift insertion, resulting in a truncated protein of 1,245 amino acids. A frameshift mutation due to a 1-nt deletion was found in family UR014, resulting in a truncated protein of 1,280 amino acids. Family UR015 had a nonsense mutation, R643X, and family UR016 had a missense mutation, G727R, in a highly conserved amino acid of domain 3. The patient with PHS had a nonsense mutation, E1147X. These results add two phenotypes to the phenotypic spectrum caused by GLI3 mutations: the combined PAP-A/B and PPD-IV. These mutations do not support the suggested association between the mutations in GLI3 and the resulting phenotypes. We propose that all phenotypes associated with GLI3 mutations be called "GLI3 morphopathies," since the phenotypic borders of the resulting syndromes are not well defined and there is no apparent genotype-phenotype correlation.

The gene for autosomal dominant hidrotic ectodermal dysplasia (Clouston syndrome) in a large Indian family maps to the 13q11-q12.1 pericentromeric region

Hidrotic ectodermal dysplasia (HED), Clouston syndrome (MIM No. 129500), is an autosomal dominant disorder affecting the skin and its derivatives. It is characterized by alopecia, dysplastic nails in hands and feet, and hyperkeratosis of the palms and soles. We have studied a large Indian pedigree (UR005), from Gujarat region, consisting of a total 127 individuals including 41 affected (12 males and 29 females). The phenotype in this family ranged from atrichosis to hypotrichosis, sparsity or absence of eyebrows, and thickening of palms and soles. In order to map the disease locus by linkage analysis, DNA polymorphisms were used in DNAs from 23 affected and 8 normal individuals. While genotyping was in progress, Kibar et al. [1996] reported mapping of the locus of a similar disease in French-Canadian families to 13q around marker D13S141. We then utilized markers on 13q to genotype the members of the Indian family. Linkage with 13q11-12.1 markers was confirmed with a maximum lod score of 3.27 (theta=0.00) with locus D13S1316. Multipoint linkage analysis yielded a lod score of 5.04 at theta=0.00 with D13S1316; haplotype analysis indicated that the gene for the Clouston syndrome in this family is localized proximal to D13S292. These data suggest that the gene for the Clouston syndrome in this Indian pedigree is probably the same as that described in the French Canadian families. The combination of data from all available families linked to 13q11-12.1 will make it possible to narrow the critical region and facilitate the positional cloning of the elusive gene.

Mapping one form of autosomal dominant postaxial polydactyly type A to chromosome 7p15-q11.23 by linkage analysis

Postaxial polydactyly type-A (PAP-A) in humans is an autosomal dominant trait characterized by an extra digit in the ulnar and/or fibular side of the upper and/or lower extremities. The extra digit is well formed and articulates with the fifth, or extra, metacarpal/metatarsal, and thus it is usually functional. In order to map the gene responsible for PAP-A, we studied a five-generation Indian family of 37 individuals (15 of whom were affected). A genomewide search with highly informative polymorphic markers on part of the pedigree showed linkage between the PAP-A phenotype and markers on chromosome 7p15-q11.23 (no crossovers were found with D7S526, D7S795, D7S528, D7S521, D7S691, D7S667, D7S478, D7S1830, D7S803, D7S801, or ELN). The highest LOD score was obtained with marker D7S801 (zeta max = 4.21; theta = 0). Haplotype analysis enabled the mapping of the PAP-A phenotype in this family between markers D7S2848 and D7S669. Analysis of additional families with PAP-A will narrow down the critical genomic region, facilitate positional cloning of the PAP-A gene, and/or uncover potential genetic heterogeneity.

An autosomal dominant triphalangeal thumb: polysyndactyly syndrome with variable expression in a large Indian family maps to 7q36

Hereditary developmental abnormalities of the upper or lower limbs in humans are easily recognizable phenotypes that can be used in the mapping and cloning of genes involved in normal human development. We studied a large Indian pedigree (UR002) with an autosomal dominant triphalangeal thumb (TPT) and polysyndactyly (PSD). The abnormalities were present only in the upper limbs, and the phenotype was fully penetrant. The expression of the phenotype was variable and ranged from unilateral TPT to bilateral TPT, preaxial du-, tri-, or quadruplication of the thumb, or syndactyly of multiple thumbs. There were 112 affected individuals in the pedigree. Previous linkage analyses on apparently similar phenotypes have identified a locus at 7q36 [Heutink et al., 1994, Nature Genet 6:287-291; Tsukurov et al., 1994]. To map the gene responsible for the TPT-PSD in family UR002, we performed linkage analysis in DNA from 47 affected and 7 normal individuals. Marker D7S550, located at 7q36, yielded a maximum LOD score of 11.31 at theta = 0.00. Additional markers in the region also showed no recombination. These data indicate that the gene responsible for the hand abnormality in pedigree UR002 maps to the same region as that in previous pedigrees with similar phenotype. Further analyses of recombinants among all the linked families by using new polymorphic markers will narrow the critical genomic region and facilitate positional cloning of the elusive gene.

Polydactyly: a study of a five generation Indian family

Preaxial polydactyly was observed in up to five generations of an Indian family living in a village in the Rajkot district (Gujarat). Among the 71 affected members, 45 were males and 26 were females. All these affected members showed preaxial polydactyly manifesting as a well formed, articulated extra digit of the hand or foot. Twenty other cases were also identified with polydactyly involving triphalangeal digits replacing the thumbs or duplication of the big toe(s). To the best of our knowledge, the present family is the largest in which several members have preaxial polydactyly of different types. No other abnormalities were apparent. The present study strongly suggests that preaxial polydactyly with a well formed extra digit, triphalangeal thumbs, and duplication of the big toe can be manifestations of the same autosomal dominant gene. It is likely that other factors are modifying the expression of this gene.